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Preface 


Renewable energy is a key challenging problem that is increasingly gaining atten- 
tion worldwide. Renewable energy sources, such as solar and wind, are promising 
alternatives to conventional fossil fuels because they are clean, sustainable, safe, and 
environmentally friendly with zero CO2 emissions. For instance, it has been shown 
that 100 GW of photovoltaic (PV)-generated power in Europe in 2012 kept more 
than 53 million tons of CO2 from being emitted into the atmosphere. On the other 
hand, the global capacity of wind power has reached 539 GW in 2017, with an 
increase of 52.5 GW compared to 2016. According to the World Wind Energy 
Association, the overall capacity of all wind turbines installed by the end of 2018 
attained 600 GW. Furthermore, solar PV and wind energy, which are sustainable 
and economically competitive renewable energy sources, increase a countries' 
energy security by reducing their dependence on fossil fuels. Hence, the require- 
ment to operate these sources becomes a mandatory obligation to reduce environ- 
ment pollution and participate in sustainable energy development. 


However, the main crucial and challenging issue in solar energy production is the 
volatility of intermittent power generation due mainly to weather conditions. In 
particular, a variation of the temperature and irradiance can have a profound 
impact on the quality of electric power production. A drop of more than 20% of 
power PV production can be observed in real PV energy plants. This fact usually 
limits the integration of PV systems into the power grid. Hence, accurately fore- 
casting the power output of PV modules in the short-term is of great importance for 
daily/hourly efficient management of power grid production, delivery, and storage, 
as well as for decision-making on the energy market. Such management aims to take 
suitable control actions to balance the electricity supply and the large electricity 
demand, as well as to meet an electricity generation performance complying with 
power quality requirements and standards. 


Despite the progress in renewable energy systems to ensure their optimal operation, 
in practice several factors can significantly affect their performance by decreasing 
their efficiency. Indeed, PV systems and wind turbines are frequently exposed to 
different sources of faults and failures that affect the power generated by these 
systems. For example, ina PV system, many potential external interferences or 
faults can cause significant losses in power production, such as dust accumulation 
on the PV modules and shading. These faults could considerably reduce both their 
production efficiencies and lifespan. Therefore, it is crucial to detect and identify 
possible faults or failures in renewable energy systems (i.e., PV systems and wind 
turbines) as early as possible. Hence, accurate and prompt fault detection and 
diagnosis are required to improve renewable energy systems efficiencies, avoid the 
high cost of maintenance, and reduce risks of fire hazards, which could affect both 
personnel and installed equipment. 


The main objective of this book is to study and develop efficient statistical model- 
ing, forecasting, and fault detection techniques for renewable energy systems. To 
this end, different methods are presented to model and forecast power production 
using a model-based or data-based framework. Moreover, innovative techniques 


are proposed for fault isolation and control in renewable energy systems. The book 
is divided into ten chapters. These chapters are regrouped into three parts: fore- 
casting, modeling, and fault detection and control. 


We would like to thank Marijana Francetic, Author Service Manager, for her con- 
tinuous assistance during the preparation of this book. We would like to gratefully 
acknowledge the financial support by funding from King Abdullah University of 
Science and Technology (KAUST), Office of Sponsored Research (OSR) under 
Award No: OSR-2019-CRG7-3800. 


Fouzi Harrou and Ying Sun 

King Abdullah University of Science and Technology (KAUST) Computer, 
Electrical and Mathematical Sciences and Engineering (CEMSE) Division, 
Thuwal, Saudi Arabia 
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Forecasting in Renewable 
Energy Systems 


Chapter 1 


Forecasting of Photovoltaic Solar 
Power Production Using LSTM 
Approach 


Fouzi Harrou, Farid Kadri and Ying Sun 


Abstract 


Solar-based energy is becoming one of the most promising sources for producing 
power for residential, commercial, and industrial applications. Energy production 
based on solar photovoltaic (PV) systems has gained much attention from 
researchers and practitioners recently due to its desirable characteristics. However, 
the main difficulty in solar energy production is the volatility intermittent of pho- 
tovoltaic system power generation, which is mainly due to weather conditions. For 
the large-scale solar farms, the power imbalance of the photovoltaic system may 
cause a significant loss in their economical profit. Accurate forecasting of the power 
output of PV systems in a short term is of great importance for daily/hourly effi- 
cient management of power grid production, delivery, and storage, as well as for 
decision-making on the energy market. The aim of this chapter is to provide reliable 
short-term forecasting of power generation of PV solar systems. Specifically, this 
chapter presents a long short-term memory (LSTM)-based deep learning approach 
for forecasting power generation of a PV system. This is motivated by the desirable 
features of LSTM to describe dependencies in time series data. The performance of 
the algorithm is evaluated using data from a 9 MWp grid-connected plant. Results 
show promising power forecasting results of LSTM. 


Keywords: forecasting, deep learning, LSTM, solar power production 


1. Introduction 


Solar energy becomes one of the most promising sources for generating power 
for residential, commercial, and industrial applications [1, 2]. Solar photovoltaic 
(PV) systems use PV cells that convert solar irradiation into electric power. 
Renewable energy sources, in particular photovoltaic (PV) energy, has been pro- 
gressively increased in recent years because of its advantages of being plentiful, 
inexhaustible, clean energy and environmentally friendly [3-5]. As one of the most 
popular renewable energy sources, solar energy has the advantages of abundant 
resources, no pollution, free use, and no transportation [6-8]. This greatly acceler- 
ated the installation of solar photovoltaic (PV) systems around the world. 

Reliable and precise forecasting plays an important role in enhancing power 
plant generation based on renewable energy sources such as water, wind, and 
sun [9]. One of the most sustainable and competitive renewable energy sources is 
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solar photovoltaic (PV) energy which is becoming nowadays more attracting 
than ever before [3]. The main crucial and challenging issue in solar energy 
production is the volatility intermittent of PV system power generation due to 
mainly to weather conditions. In particular, a variation of the temperature and 
irradiance can have a profound impact on the quality of electric power production. 
A drop of more than 20% of power PV production can be observed in real PV 
energy plants. This fact usually limits the integration of PV systems into the 
power grid. Hence, accurately forecasting the power output of PV modules in 
a short-term is of great importance for daily/hourly efficient management of 
power grid production, delivery, and storage, as well as for decision-making on 
the energy market [10]. 

Precise forecasting of solar energy is important for photovoltaic (PV) based 
energy plants to facilitate early participation in energy auction markets and efficient 
resource planning [11]. Numerous methods have been reported in the literature 
for PV solar power forecasting. These methods can be classified into four classes: 
(i) statistical approaches based on data-driven formulation to forecast solar time 
series by using historical measured data, (ii) machine learning techniques, in par- 
ticular, deep learning approaches based artificial neural network, (iii) physical 
models based on numerical weather prediction and satellite images, and (iv) hybrid 
approaches which are the combination of the above methods. In [12], a combined 
approach merging seasonal autoregressive integrated moving average (SARIMA), 
random vector functional link neural network hybrid model and discrete wavelet 
transform has been introduced for forecasting short-term solar PV power produc- 
tion. It has been shown that the combined models provide improved forecasting 
results compared to individuals ones. In [13], Gradient boosted regression trees 
approach has been used to predict solar power generation for 1-6 h ahead. It has 
been that this approach outperforms the simpler autoregressive models. In [14], a 
model combining seasonal decomposition and least-square support vector regres- 
sion has been designed to forecast power output. This approach demonstrated good 
forecasting capacity compared to the autoregressive integrated moving average 
(ARIMA), SARIMA, and generalized regression neural network. In [15], a multi- 
variate ensemble forecast framework integrating ensemble framework with neural 
predictors and Bayesian adaptive combination is proposed for forecasting PV 
output power. 

Most conventional solar power forecasting approaches are limited in uncovering 
the correlation of the limited data but are not able to deep correlation and uncover 
implicit and relevant information. With the huge data from the modern power 
system, the use of conventional approaches is not suited for guaranteeing precise 
forecasting. Recently, deep learning (DL) approaches have emerged as powerful 
machine learning tools that enable complicated pattern recognition and regression 
analysis and prediction applications [16-18]. DL approaches are becoming increas- 
ingly popular due to their good capacity in describing dependencies in time series 
data. Deep Learning is the result of the concatenation of more layers into the neural 
network framework. Over the past few decades, many deep learning models have 
been proposed including Boltzmann machines, Deep Belief Networks (DBN) and 
Recurrent Neural Networks (RNNs) [19]. RNN is a type of neural networks that 
exploits the sequential nature of input data. RNNs are used to model time- 
dependent data, and they give good results in the time series data, which have 
proven successful in several applications domains [3, 20, 21]. Long Short-Term 
Memory Networks (LSTM) is a type of RNNs that is able to deal with remembering 
information for much longer periods of time [22]. It is also considered as one of the 
most used RNN models for time series data predictions, which is perfectly suited 
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to PV solar power production forecasting problems. In this chapter, we applied the 
LSTM model to accurately forecast short-term photovoltaic solar power. The 
effectiveness of this approach is tested based on power output data collected from 
a 9 MWp grid-connected plant. 

The next section introduces the core idea behind the LSTM model and how it 
can be designed and implemented. Then, Section 3 presents the results of solar 
photovoltaic power forecasting using the LSTM model. Lastly, conclusions are 
offered in Section 4. 


2. Deep learning and forecasting of PV power production 


Over the last decades, many studies have been dedicated to forecasting problems 
in several application domains. Recurrent Neural Networks (RNNs) have been 
successfully used in machine learning problems [23]. These models have been 
proposed to address time-dependent learning problems [22]. Figure 1 shows the 
basic concept of RNNs; a chunk of a neural network, A, looks at some input x, and 
outputs a value h,. It should be noted that RNNs are suited to learn and extract 
temporal information [24]. A general formula for RNN hidden state h t given an 
input sequence x = (X1, X2, ..., Xt): 


n= |? t=0 a) 


(Wx, Xt), otherwise 


where ọ is a non-linear function. The update of recurrent hidden state is 
realized as: 


hi =8(Wx, + uhy-1) (2) 


where g is a hyperbolic tangent function (tanh). 

Generally, it is not easy to capture long term time dependencies in time series 
when using recurrent neural networks. To bypass this limitation, Long Short-Term 
Memory Networks (LSTM) models were designed. LSTM is an extended version 
of RNN that are effectively capable to handle time dependency in data [22]. These 
models are flexible and efficient to describe time-dependent data, and they dem- 
onstrated success in several applications. LSTM is one of the most used RNN models 
for time series data predictions, which is perfectly suited to the PV forecasting 
problems [22]. Next, we present a basic overview of LSTM and how it can be 
designed and implemented. 


RNN hidden state 
he 


Input sequence x X: 
X=( X1, X2,... Xt) 


Figure 1. 
Basic illustration of RNN. 
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Figure 2. 
Illustration of LSTM unit. 


2.1 Long short-term memory (LSTM) models 


The Long Short-Term Memory (LSTM) is a variant of the Recurrent Neural 
Networks (RNN) that is capable of learning long term dependencies. LSTM models 
were initially proposed by Hochreiter and Schmidhuber [4] and were improved and 
popularized by many other researchers [4—6, 9]. LSTM models have an excellent 
ability to memorize long-term dependencies, are developed to deal with the 
exploding and vanishing gradient problems that can be encountered when training 
traditional RNNs. Relative insensitivity to gap length is an advantage of LSTM 
models over ANNs models, hidden Markov models and other sequence learning 
methods in several application domains. 

A common LSTM model is composed of cell blocks in place of standard neural 
network layers. These cells have various components called the input gate, the 
forget gate and the output gate. The cell remembers values over arbitrary time 
intervals and the three gates regulate the flow of information into and out of the cell 
[5]. Figure 1 shows the basic structure of RNN-LSTM. 

From Figure 2, the RNN-LSTM has two input features at each time, which 
include the current time step input X, (input vector) and the hidden state of the 
previous time step H,_, (previous input vector). The output is computed by the 
fully connected layer with its activation function (e.g., tanh, sigmoid, Softmax, and 
Adam). Therefore, the output of each gate can be obtained through logical opera- 
tion and nonlinear transformation of input. 

Let us denote the input time series as X,, the number of hidden units as h, the 
hidden state of the last time step as H,_1, and the output time series as Ht. 

The mathematical relationship between inputs and outputs of the RNN-LSTM can 
be described as follows. 


I, = o(XW și + Hp1W); + bi) (3) 
F; = 0(X:Wop + Hi-1Wyf + bp) (4) 
O; = o(XtW xo + Hi-1W no + bo) (5) 

C, = tanh (X;Wee + Hi-1Wpc + be) (6) 
C, = F, 0 C10 Č, (7) 
H, = O; o tanh (C;) (8) 


where 


° I,, F;, O; are input gate, forget gate, and output gate respectively, W xi, Wag, 
Wxo and Whi, Wig, Who are weight parameters and b;, bp, bo are bias 


Forecasting of Photovoltaic Solar Power Production Using LSTM Approach 
DOI: http://dx.doi.org/10.5772/intechopen.91248 


parameters. All these gates have the same dimensions and the same equations 
just with different parameters. They are called gates because the activation 
function transforms the element values between ranges ([0, 1], [—1, 1]). The 
input gate defines how much of the newly computed state for the current input 
you want to let through. The forget gate defines how much of the previous 
state you want to let through. The output gate defines how much of the 
internal state you want to expose to the external network (higher layers and 
the next time step). 


e Č; is the candidate memory cells, Wye, Wy, are weight parameters and b c isa 
bias parameter. LSTM model needs to compute the candidate memory cell C;, 
its computation is similar to the three gates (input, forget and output gates), 
but using a tanh function as an activation function with a value range between 
[—1, 1]. 


e C; is the memory cells, o is an operator that expresses element-wise 
multiplication. The computation of the current time steps memory cell C; 
combines the information of the previous time step memory cells (C;_1) and 
the current time step candidate memory cells (C,), and controls the flow of 
information through forgetting gate and input. 


° Hi is the hidden states, we can control the flow of information from memory 
cells to the hidden state H; through the output gate. The tanh function ensures 
that the hidden state element value is between [—1, 1]. It should be noted 
that when the output gate is approximately 1, the memory cell information 
will be passed to the hidden state for use by the output layer; and when the 
output gate is approximately 0, the memory cell information is only retained 
by itself. 


2.2 Proposed approach 


The proposed approach in this chapter aims to forecast solar power production. 
This methodology is based on the LSTM deep-learning model. Figure 3 summarizes 
the main steps of the proposed methodology. The proposed approach includes four 
key steps (Figure 3): 


i. Collect the SCADA data from the PV system. 


ii. Pre-process and clean data by removing outliers and imputing missing 
values. 


iii. Normalize the original data. 


iv. Train, validate and test the LSTM model. Various statistical indicators are 
used to quantify the accuracy of the developed model. Lastly, the designed 
LSTM model can be used for power production forecasting. 


2.3 Metrics for evaluating the forecasting models 


To assess the forecasting performance, numerous statistical indicators have been 
proposed in the literature including root mean square error (RMSE), mean absolute 
error (MAE), coefficient of determination (R°), and mean absolute percentage 
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Photovoltaic System 


Data collection 


Data preprocessing 


Data standardization 


Training and 


validation data 
Improve the 
~ performance of 
Test data Build and train ; 
È LSTM model 
i 


LSTM model 


Test and validate 
the model 


Figure 3. 
Schematic block of the proposed forecasting method. 


error (MAPE). In this study, we used R? and MAPE, which are frequently to 
evaluate the forecasting accuracy: 


MAE = Yf -x| (9) 
a 2 

RMSE = 2x) (10) 
24 pa: — #4) 

wes Goa? a 


where x are the measured values, £ are the corresponding forecasted values by 
the LSTM model and x is the number of measurements. 


2.4 Implementation steps 


Essentially, the LSTM model can be designed and implemented in four main 
steps. At first, define the LSTM model and train it, then fit the LSTM model, and 
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lastly, the trained LSTM model is used for forecasting. Table 1 summarizes the 
main steps (partial codes) performed in designing the LSTM model. 


2.5 Enhance LSTM models performance 


The key factors impacting the accuracy of the LSTM model are not only the 


amount of training data but also the architecture of the network, hyper-parameters 


and the utilized optimizers. Accordingly, the performance of LSTMs can be 
enhanced by acting on the following elements. 


e Activation functions: activation functions an important role in determining the 


final response of the neural network. Two families of functions are 
distinguished: linear and nonlinear functions. The output of the linear 
activation functions is linearly proportional to the inputs and is not limited 


between any ranges. They are more suited than a step function because they 
permit obtaining multiple outputs, not just binary output (i.e., yes and no). On 
the other hand, nonlinear activation functions are the most frequently utilized 
because they are flexible and permit obtaining nonlinear output and they are 
confined within a range. For instance, Sigmoid, Softmax, and Rectified Linear 
Unit (ReLU) activation functions permit rescaling the data to values in the 
interval [0, 1], while Hyperbolic Tangent (tanh) activation functions rescale 
the data within [—1, 1]. 


e Optimizer: in the training phase of the LSTM model, optimization algorithms 
are used for minimizing its error rate. The performance of an optimizer is 
generally characterized by convergence speed and generalization (the 
efficiency of the model on new datasets). The commonly used optimizers 
include Adaptive Moment Estimation (Adam) or Stochastic Gradient Descent 
(SGD) [25, 26]. 


e Dropout: it is a well-known stochastic regularization procedure applied to avoid 
overfitting and further enhance the prediction capacity of RNN models [27]. 
More details about dropout techniques can be found in [26-28]. 


¢ Epochs and batches: the number of epochs and batch are two important 
parameters when constructing deep learning models. It has been shown in the 


Steps Action 
Step 1: Define LSTM from keras.layers.recurrent import LSTM 
network from keras.models import Sequential 


from keras.layers.core import Activation, Dense, Dropout 

model = Sequential () 

model.add(LSTM (units=nb_neural, return_sequences=True, input_shape= 
(Xtrain.shape[1], 1))) 


Step 2: Compile the LSTM — model. compile(loss="mse", optimizer="adam", metrics=[rmse, ‘mae’, 
network Rsquare] ) 


Step 3: Fit the LSTM history = model.fit(Xtrain, ytrain, batch_size=batch_size, 
network epochs=num_epochs, 
validation_data=(Xval, yval), verbose=2) 


Step 4: Forecasting ypred = model.predict(Xtest) 


Table 1. 
Partial codes used for building the LSTM network. 
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literature that good results can be achieved when using large epochs and small 
batch sizes. 


© Weight regularization: another way to avoid overfitting and improves model 
performance is called weight regularization. This approach imposes constraints 
on the RNN weights within nodes to allow the network to maintain the weights 
small. Several penalizing or regularization approaches are commonly used in 
the literature based on L1 or L2 vector norm penalty. 


3. Results and discussion 


This study is based on real data collected from January 2018 to December 2018 
every 15 min from a 9 MWp grid-connected plant. Figure 4 shows the hourly 
distribution of PV power production day from January 2018 to December 2018. 


DC PV power production (kW) 


309 


10 ~ 200 
=- 100 
Hours 5 0 Days 


Figure 4. 
PV power production per hour for each day from January 2018 to December 2018. 
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Figure 5. 
Distribution of DC power output in the daytime. 
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From Figure 4, the solar PV power production reaches, every day, its maximum at 
the mid-day and falls to zero over the night. 

For the data in Figure 4, the box plots showing the distribution of DC power 
generation in the daytime are displayed in Figure 5. One can see that the maximum 
power production is achieved around mid-day. 

The monthly cumulative DC power generated by the inspected PV system from 
January 2018 to December 2018 is displayed in Figure 6. The highest and lowest 
monthly cumulative power are respectively achieved in March (6450.056 MW) and 
in October (4655.524 MW). 

Figure 7 shows the monthly distribution of DC power production during the 
monitored period. Figure 7 shows that the produced DC power is relatively high in 
January, February, and March. Also, it can be noticed that the production was 
relatively low from June to September (Figure 7). 

To investigate the interactions between the DC power and meteorological fac- 
tors (i.e., inclined irradiance 27, ambient temperature, and wind velocity) a Pearson 
correlation heatmap is displayed in Figure 8. From Figure 8, one can see that there 


7000 


6000 


5000 
4000 
3000 
2000 
1000 

0 


Jan Feb Mar Apr May Jim Jul Aug Sep Oct Nov Dec 


PV Power Production (MW) 


Figure 6. 
Monthly total DC power produced from January 2018 to December 2018. 
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Figure 7. 
Monthly distribution of DC power output. 
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Figure 8. 
Heatmap of the correlation matrix of data: inclined irradiance 27, ambient temperature, wind velocity, and 
power. 


Figure 9. 
ACF DC power measurements. 


Models | Parameter Value 
Number of LSTM layers 1 
learning: rate 0.001 
Loss function Cross-entropy 


Nuniber of hidden units 200 
Optimizer Adam 


Number of training epochs 100 


Table 2. 
Parameters in LSTM model. 
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is a high correlation between solar irradiance and power production. It should be 
noted that DC Power has a low correlation with wind velocity. 

Figure 9 shows the autocorrelation function (ACF) plot of the data shown in 
power generation data. A seasonality of 24 h can be seen from the ACF plot of PV 
power data, the time difference between two maximum in the ACF (Figure 9). In 
particular, this seasonality is mainly due to the variation of solar irradiance. 

The LSTM model has been constructed and then used for forecasting. Data were 
split into training and testing datasets (90% and 10% respectively). Parameters of 
the constructed LSTM are presented in Table 2. 

The evolution of the loss function and RMSE in the function of the number of 
iterations is displayed respectively in Figures 10 and 11. Figures 10 and 11 indicate 
the convergence of the loss function and RMSE when the number of epochs is 
around 60. 

Once the LSTM model has been constructed based on training data, it will be 
employed to forecast future values of power production. We attempt now to test 
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Figure 10. 
Evolution of LSTM loss function during training stage. 
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Figure 11. 
Evolution RMSE of LSTM model during training stage. 
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the capability of the above LSTM model to forecast future values of the PV power 
generation. Figure 12 shows the forecasting results of the PV power generation 
compared with the real data over a time horizon. Figure 13 shows the scatter plot of 
the measured and forecasted power production via the LSTM model. It can be seen 
from Figures 12 and 13 that the computed LSTM model has the ability of short-term 
forecasting of PV power generation. In addition, the forecasting result in Figure 12 
illustrates the efficiency of the LSTM model to forecast PV power production even 
under a cloudy day (i.e., the second day in Figure 12) where the power data is very 
dynamic. 

In summary, the LSTM model showed good forecasting capacity with the coef- 
ficient of determination R° = 0.98 close to 1 and relatively small mean absolute 
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Figure 12. 
Plot of collected solar power and forecasted one using LSTM model. 
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Figure 13. 
Scatter graph of measured and LSTM forecast solar power output. 
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percentage error (MAPE), MAPE = 8.93. It should be pointed out that the 
forecasting accuracy in cloudy days could be improved by including meteorological 
variables, such as solar irradiance, ambient temperature, and wind velocity, as 
input variables. 


4. Conclusion 


The major challenge in solar energy generation is the volatility intermittent of 
photovoltaic system power generation due mainly to weather conditions. Thus, 
accurate forecasting of photovoltaic power generation is becoming indispensable 
for reducing the effect of uncertainty and energy costs and enable suitable integra- 
tion of photovoltaic systems in a smart grid. This chapter employed a Long Short- 
Term Memory (LSTM) model to accurately forecast short-term photovoltaic solar 
power. This approach exploits the desirable properties of LSTM, which is a power- 
ful tool for modeling dependency in data. The forecasting quality of this approach 
has been verified using data from January 2018 to December 2018 collected from a 9 
MWp grid-connected plant. Promising results have been achieved by the proposed 
LSTM-based approach to short-term forecasting of photovoltaic solar power pro- 
duction. As future work, to further enhance the forecasting quality we plan to 
implement and test the performance of other RNN models like Gated recurrent unit 
(GRU) model and to incorporate other information such as meteorological data. 
Also, as most data from real plants are multiscale in nature and noisy, we plan in 
future work to merge the desirable LSTM model with the wavelet-based multiscale 
presentation [29]. This permits to get a multiscale LSTM model able to capture 
feature in both time and frequency and possess good ability to handle noisy data. 
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Optimal Design of Energy System 
Based on the Forecasting Data 
with Particle Swarm Optimization 


Yamin Yan, Haoran Zhang, Jianqin Zheng and Yongtu Liang 


Abstract 


Renewable energy source has developed rapidly and attracted considerable 
attention. The integration of renewable energy into the energy supply chain 
requires precise forecast of the output of energy supply chain, thereby reducing 
energy resource waste and greenhouse gas emissions. In this study, a coupled model 
system is developed to forecast energy supply chain for the design optimization of 
distributed energy system, which can be divided into two parts. In the first part, 
long short-term memory (LSTM) and particle swarm optimization algorithm (PSO) 
contribute to energy supply chain forecast considering time series, and particle 
swarm optimization is used to optimize the parameters of the long short-term 
memory model to improve the forecast accuracy. Results show that the mean 
absolute error and root mean squared error are 8.7 and 16.3 for the PSO-LSTM 
model, respectively. In the second part, the forecast results are used as input of the 
distributed energy system to further optimize the design and operation schemes, so 
as to achieve the coupling optimization of forecast and design. Finally, a case study 
is carried out to verify the effectiveness of the proposed method. 


Keywords: coupled model system, forecast, design optimization, renewable 
energy system 


1. Introduction 


With the development of social economy and the acceleration of urbanization 
and industrialization, the global energy consumption and greenhouse gas (GHG) 
emissions are expected to increase in the next few decades [1, 2]. At present, more 
than 80% of the world’s primary energy comes from fossil fuels, and its carbon 
emissions are the biggest cause of global greenhouse effect. Therefore, the develop- 
ment of low-carbon economy renewable energy technologies is an emerging trend 
and goal of the world. In the past decade, renewable energy power generation has 
grown at an annual rate of 16%, and China accounted for 45% of global growth [3]. 
Wind power, one of the most potential renewable energy sources, has entered the 
rapid development stage. In 2018, the wind power has a largest share (around 50%) 
in the renewable energy power generation and reached 142 TW h [4], and more 
than 10% of electricity consumption or 5 GW of electricity will be generated by 
solar energy by 2020 [5]. As reported in “Clean Energy Trends” published in March 
2011, the global production of biofuels, wind power, and solar energy is expected to 
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reach $112.8 billion, $122.9 billion, and $113.6 billion by 2020, with growth rates of 
116, 103, and 60%, respectively. Whether from the point of view of energy saving or 
environmental protection, renewable energy is certain to play a vital role in power 
generation [6, 7]. However, renewable energy has the characteristics of strong 
uncertainty and volatility, which not only increase the operation cost of power system 
but also reduce the cost-effectiveness of renewable energy resources [8]. Therefore, it 
is necessary to establish an accurate and stable prediction model for the renewable 
energy power to ensure the energy supply chain, so as to lay the foundation for 
optimizing the design and operation strategies of energy system [9, 10]. 

Many scholars have conducted the research on forecasting and optimizing the 
renewable energy supply system. Regarding the energy supply system forecast, 
many studies have developed effective forecasting methods for specific energy 
sources, which can be divided into three types of models, namely, time series 
models, data-driven models, and wavelet-based multiscale models [11]. Generally, 
the time series models and data-driven models are widely accepted. For instance, Lu 
et al. [12] compared the accuracy of four artificial intelligence methods in forecast- 
ing Taiwan’s renewable energy sources based on historical data from 2000 to 2015, 
and results showed that only grey forecasting model coupled with heuristic fuzzy 
time series method is suitable for small dataset forecast. Long et al. [13] used four 
famous data-driven method to forecast daily solar power after dimensionality 
reduction of data, and they inferred that different algorithms can outperform others 
in different considered scenarios. There are few literatures on the application of the 
wavelet-based multiscale models. Reikard [14] proposed a Kalman filter and time- 
varying regression for the time series prediction of wave energy. Aasim et al. [15] 
proposed a new repeated wavelet transform based ARIMA model for very short- 
term wind speed forecasting. Based on data-driven, Li et al. [16] developed a new 
deep machine learning algorithm to predict short-term wave energy. In summary, 
different forecast models may vary greatly in different scenarios. 

However, most of the previous work has been devoted to the forecast of renew- 
able energy in a certain region considering the limited range of measuring instru- 
ments, and few studies pay attention to the power generation in multiple regions. In 
fact, it is impossible to install the measuring instruments in each region from an 
economic point of view, thus bringing a great challenge for forecasting renewable 
energy in multiple regions. In addition, the forecast of renewable energy is not only 
related to the energy intensity but also depends on climate conditions, which have a 
strong time series relationship [17]. Meanwhile, the short stochastic characteristics 
and the dependence of the observation time series of renewable energy must be 
taken into account. Therefore, this paper is the first work to develop a forecast 
model based on the data-driven method by collecting the measurement equipment 
data and climate conditions in different regions and at different times. Also, the 
LSTM model is adopted to capture the dependence of time series forecasting for 
renewable energy output, and PSO is used to optimize the parameters of the LSTM 
model to improve the forecast accuracy. 

Regarding the optimization of renewable energy supply system, many scholars 
are devoted to integrating various renewable energy sources and have made a range 
of achievements [18-20]. Acevedo-Arenas et al. [21] described the methodology of 
hybridizing photovoltaic, wind, and forest biomass energy sources and developed a 
model to optimize the design and operation schemes of hybrid renewable energy 
system in consideration of economic and environmental impacts. Tajeddin and 
Roohi [22] put forward the use of biomass energy to improve the responsibility and 
efficiency of wind farms. Sakaguchi and Tabata [23] forecasted the power genera- 
tion potential and proportion of biomass, wind power, and PV in Awaji Island 
Japan, and results showed that it is possible to achieve self-sufficiency by the sole 
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use of renewable energy. Moreover, the renewable energy plants are usually 
equipped with the energy storage, which could improve the controllability of the 
output [24]. Based on the above research, the hybridization of renewable energy 
sources is able to reduce the uncertainty and volatility, thereby promoting the 
rational use of renewable energy and increasing environmental benefits. 

In this paper, we first focus on the accurate forecast of renewable energy supply 
chain, which is divided into three steps. The first step is to develop a forecast model 
based on the data-driven method by collecting the measurement equipment data 
and climate conditions in different regions and at different times. In the second 
step, the LSTM model is adopted to capture the dependence of time series forecast- 
ing for renewable energy output. In the third step, PSO is used to optimize the 
parameters of the LSTM model to improve the forecast accuracy. After that, the 
forecast results are used as input of the distributed energy system to further opti- 
mize the design and operation schemes, thereby achieving the coupling optimiza- 
tion of forecast and design. Overall, there are four contributions in this paper: 


1.LSTM model is established to forecast the output of renewable energy in 
multiple regions, considering the correlation of time series and the change of 
climate conditions. 


2.PSO is adopted to optimize the parameters of LSTM model, thereby improving 
the forecast accuracy. 


3. Based on the forecast results of LSTM model, a MILP model is established to 
optimize the design and operation schemes of renewable energy supply system. 


4. Real multi-region renewable energy plants are used as an example to verify the 
effectiveness and practicality of the proposed coupled model system. 


The rest of this paper is organized as follows: Section 2 describes the basic 
knowledge of LSTM model and PSO algorithm, and the optimization model of 
distributed renewable energy supply system is explained. Subsequently, the optimal 
results solved by the coupled model system and corresponding discussions are 
presented in Section 3. Finally, Section 4 concludes this paper. 


2. Problem formulation and solution methodology 
2.1 LSTM model coupled with PSO algorithm 


Figure 1 shows the framework of the coupled model system. Firstly, long short- 
term memory (LSTM) and particle swarm optimization algorithm (PSO) contribute 


Solar power plants 


c — 
| =A ar aed i a 
j ee, 


í LSTM E sod sess RES comfguration 
Solar meters Outahase L maa a J predic 


Figure 1. 
The framework of the coupled model system. 
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to energy supply chain forecast considering time series, and particle swarm optimi- 
zation is used to optimize the parameters of the long short-term memory model to 
improve the forecast accuracy. Then, the forecast results of solar power plants and 
the known wind power are used as input of the distributed energy system to further 
optimize the design and operation schemes, so as to achieve the coupling optimiza- 
tion of forecast and design. 


2.2 LSTM model coupled with PSO algorithm 
2.2.1 LSTM model 


As a kind of recurrent neural network, LSTM model combines short-term with 
long-term memory through subtle gate control and overcomes the shortcoming of 
gradient vanishing. Currently, LSTM model has great applications in time series, 
such as machine translation and speech recognition. In this section, the basic LSTM 
cell is introduced firstly. 

The basic LSTM cell is shown in Figure 2. At time t, the cell state at t—1, the 
hidden state at t—1, and the new information are denoted as ¢;, h,1, and xs, 
respectively, which are the input of the cell. Forget gate, input gate, and output gate 
are used to select information and depended on x; and /;.1. The output value of each 
gate is limited to between 0 and 1 by the activation function, which is the sigmoid 
function defined in Eq. (1). 


= 1 
~ +e 


o(x) (1) 

The LSTM cell included four steps as follows. Firstly, the cell decides what 
stored information is going to be thrown away from the previous cell state c;.1. 
Based on x; and h;.1, forget gate f , outputs a value between 0 and 1 for c,1. And 1 
means that the stored information is completely retained, and 0 means that the 
information is completely eliminated. f , is calculated by Eq. (2). 


f,= o(Wug Xp + Woe hea + by) (2) 


Secondly, the cell decides what new information is going to be stored. At first, 
tanh function defined in Eq. (3) is to ensure the normalization of the new informa- 
tion and converts x; and h;., to the new format between —1 and 1 denoted as in 
Eq. (4). Then, the input gate i; defined in Eq. (5) outputs a value between 0 and 1 to 


Figure 2. 
The basic LSTM cell. 
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select the new information. After that, these two parts are going to be combined to 
update the cell state in the next step. 


e—e* 
h = 
tanh (x) are (3) 
Ct ST tanh (Wee “Xe + Whe . hin T be) (4) 
i = o(Wxi -Xt + Whi -hia + bi) (5) 


Thirdly, it is going to update the cell state based on the previous steps. The old 
cell state is multiplied by f,, forgetting the information decided to be forgotten in 
step 1. The process information ¢; decided by tanh function is multiplied by i. Then 
these two parts are added to determine the new cell state c; in Eq. (6). 


Ct =f, * Cpa F ip * Cp (6) 


Finally, the cell decides what is going to be the output. Similar to other gates, 
output gate o, outputs a value to decide what parts of the cell state are going to be 
output in Eq. (7). Based on cz, the hidden state h, defined in Eq. (8) is multiplied by 
0; to decide the important information to be stored. At last, the output of the LSTM 
cell is shown in Eq. (9). 


o: = o(W xo -Xt + Who hea + bo) (7) 
h, = 0, * tanh (c+) (8) 
Yı = 0(Why «hy + by) (9) 


In this work, X is the measurement equipment data and climate conditions, and 
y is the output of the renewable energy. Wyf, Wii, Wxo, and W,, are the input 
weight matrices; Wpf, Whi, Who, and Wj,are the recurrent weight matrices; and 
Ww is the hidden output weight matrix. Vectors br, b;, bo, be, and by are the 
corresponding bias vectors. 


2.2.2 PSO algorithm 


Proposed in 1995, PSO has a great use in the field of optimization, especially 
coupling with machine learning, such as ANN [25]. In this study, PSO is used to 
select the optimal parameters for the LSTM model. The process of the algorithm is 
described as follows. 

Firstly, initializing the particle swarm, including the size of population and the 
dimension of each one. In this work, the number of units, learning rate, and time 
step are optimized, so the dimension of each individual is 3. 

Secondly, evaluating the fitness of each individual. In the process of training the 
LSTM model, the training loss is the fitness of each particle. 

Thirdly, finding the two “extreme value.” In each iteration, by calculating the 
training loss, the individual optimal value and optimal population value are found, 
which are used for updating the particles. 

Fourthly, updating the velocity and position of each particle. Based on the two 
“extreme value,” the particles are updated to form the next generation according to 
Eqs. (10) and (11): 


vilt +1) =w x vi ter x r1 x (p,(t) — xi(t)) +e2 X r2 x (p, (0) — x;(t)) (10) 


23 


Advanced Statistical Modeling, Forecasting, and Fault Detection in Renewable Energy Systems 


xilt +1) =x; (t) + vit + 1) (11) 


where v; and x; are the velocity and position of particle i, respectively. p; is the 
best position for a particle i and P, denotes the best position in the group at the tth 
iteration. 

Fifthly, reaching convergence condition. With the help of iterations, the best 
particle will be found within the maximum iteration. So the suitable parameters of 
the LSTM model are determined. 

The methodology of LSTM model coupled with PSO algorithm is composed of 
three parts. Firstly, the data of measurement equipment data and climate conditions 
are collected into the database with preprocessing. Secondly, the data is divided into 
two groups, and small part of the data is for the PSO algorithm to train the LSTM 
parameters. It contains the process of PSO algorithm mentioned in 3.1.2 section. 
Finally, after finding the optimal LSTM model, train the prediction model with the 
rest of the data. So far, a prediction model with high accuracy for the output of the 
renewable energy is established. 


2.3 The optimization model of energy supply system 
2.3.1 Objective function 


The objective of this model is to minimize the total annual cost of the energy 
supply system, Cror, defined as the sum of the total fixed cost, capital cost, and 
energy consumption cost. It should be noted that the input energy of the energy 
supply system is renewable energy sources; thus, the energy consumption cost is 
equal to zero: 


Cror = Ccac + Cr (12) 


where Ccac and Cpr are the total capital cost and fixed cost of the energy supply 
system, including the PV panels, wind turbines, and batteries [CNY]. 

The annual capital cost of technology i is obtained by the capital recovery factor 
CRF; which is defined as Eq. (15): 


Ccac = Cap;CuciCRF;i (13) 
Cr = B;CrciCRF; (14) 
1+r)™ 
A AE (15) 
(1+ry*-—1 


where Cap, is the rated power of technology i [kW, kWh]; B; is the binary 
variable, which is equal to 1 if the technology i is selected for installation [—]; Cuci 
represents the unit cost of technology i [CNY/kW for energy generation technology 
and CNY/kWh for energy storage technology]; Cc; is the fixed cost of technology i 
[CNY]; N; represents the lifetime of technology i [year]; and r defines the discount 
rate [—]. 


2.3.2 Constraints 
In this section, three types of constraints are considered: energy balances, 


design, and operation constraints. Regarding the energy balance constraints, the 
electricity demand is satisfied by the PV panels and wind turbines. When the 
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electric load changes, the batteries are used to effectively adjust the electricity 
balance and keep the voltage and frequency constant, and Eqs. (17) and (18) are 
required to describe the operation state of batteries. Notably, the batteries are only 
able to deal with the daily fluctuations, and Eq. (16) can be applied to any hour 
except for the first hour of the day, while Eq. (18) can only be applied to the first 
hour of the day: 


Q.pv + Qrwr + St — S} = D; (16) 
St = St-1M oss + Sony gg Ya /Ndis (17) 
S: = St+23Mloss + Sony a SŽ In gis (18) 


where D; is the electric load at time-window t [kWh]; Q, py and Q, wr indicate 
the energy output from the PV panels and wind turbines at time-window t [kWh]; 
S” and S% are the charging and discharging rate of the battery at time-window t 
[kWh]; S, represents the electricity stored in the battery at time-window t [kWh]; 
and Noss) Nch and ngis represent the battery’s loss, charging, and discharging 
efficiency, respectively, [—]. 

As for the design and operation constraints, the operating power of technology i 
should be less than the rated power, as defined in Eqs. (19) and (20). In addition, 
Eq. (21) requires that the installation capacity must be within the maximum limit to 
ensure safe and stable operation of technology i: 


Qr; < Cap; (19) 
Se < CAP patery (20) 
0 < Cap, < BiMax; (21) 


where Max; is the upper limit of the capacity of technology i [kW, kWh]. 


2.3.3 Model solving 


Based on the computer with 8 GB of RAM memory and solution environment 
with 1.6 GHz of CPU, the model was programmed with the software MATLAB R2015 
and solved by the business solver GUROBI 8.1.0 to obtain the global optimum. 


3. Case study 
3.1 The forecast results of LSTM model 


Firstly, this part is to predict the solar energy of the solar power plants. In area 
A, there are eight solar plants located in different regions. Historical data from 
January to May 2017 is used for experiment with a time interval of 30 minutes. The 
measurement equipment data including measured temperature and global solar 
radiation is from the power system, and the weather conditions are from the 
weather bureau. In this study, in order to evaluate the quality of the LSTM model, 
the evaluation indexes are root mean squared error (RMSE) and mean absolute 
error (MAE). The formulas of these two are defined as follows: 


1 n 
MAE = — 1; 22 
PD a (22) 
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Is 
RMSE = y DAAA (23) 


where n means the total number of the data, p, is the predicted value, and 7; is 
the real value. 

In this part, the data of January is used for training PSO-LSTM. The training sets, 
validation sets, and testing sets account for 80, 10, and 10%. And RMSE of the 
testing data is used as the fitness value of PSO. The parameters of the LSTM model 
are hidden units, learning rate, and time step. The corresponding range is 1-30, 
0.0001-0.1, and 1-50 respectively, as shown in Table 1. After executing the 
PSO-LSTM, the suitable parameters are optimized and shown in Table 1. The 
convergence of RMSE is figured in Figure 3. 

After determining the suitable LSTM model, the prediction model for solar 
energy is established. Data from February to mid-May is used to train this predic- 
tive model, and other data is used to test its accuracy. After the data is tested, the 
MAE and RMSE are 8.7 and 16.3, respectively. In order to illustrate the superiority 
of PSO-LSTM, the prediction results of basic LSTM and artificial neural network 
(ANN) are compared. And the predicted results and real value of the last 6 days of 
May are shown in Figure 4. 


3.2 The optimal results of the renewable energy system 


From the aforementioned analysis, the prediction model for solar energy is 
accurate. With the input of the measurement equipment data and weather condi- 
tions, the output of the solar energy can be predicted precisely. The hourly output 
forecast curve of the solar power plant and the known wind power energy and the 
electricity demand are shown in Figure 5. The technical and cost information of 


Neural units Learning rate Time step 
Range 1-30 0.0001-0.1 1-50 
Result 20 0.01 5 


Table 1. 
The parameters of the LSTM network. 
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Figure 3. 
RMSE of the PSO-LSTM. 
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Figure 4. 
The predicted results and the real value. 
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Figure 5. 
Forecast of solar power plant and wind farm and electricity demand. 


technologies are summarized in Table 2 [26-28], and the loss, charging, and 
discharging efficiency of battery are set as 99.9, 90, and 90%, respectively [29]. The 
discount rate is assumed to be 8% to annualize the total investment cost. 

By solving the MILP model, the optimal results of the renewable energy system 
can be obtained. The annual total cost is equal to 109.91 million. Specifically, the 


Technology Size range Efficiency Capital cost Fixed cost Lifetime 
(MW, MWh) (CNY/kW, CNY/kWh) (CNY) (a) 
PV panels 0-500 17% 1500 34,500 20 
Wind 0-500 35% 3500 17,900 25 
turbine 
Battery — — 85 11,600 20 
Table 2. 


The technical and cost information of technologies. 
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total capacities of PV panels, wind turbines, and batteries are 140, 270, and 185 
MWh. Obviously, the PV panels are selected due to the low investment cost, and 
the capacities reach the optimal to match with the limited solar power. The wind 
turbines are installed to work when the solar power is equal to 0. The batteries are 
installed to store the energy generated by the energy generation technologies to 
reduce the investment cost. 


4. Conclusion 


In this work, LSTM model is established to forecast the output of renewable 
energy in multiple regions, considering the correlation of time series and the change 
of climate conditions. Then, PSO algorithm is applied to optimize the parameters of 
the LSTM model. Results show that PSO-LSTM model is highly accurate with small 
error. Based on the forecast results, a MILP model is established to obtain the 
configuration of the renewable energy system with minimal total annual cost. 

In future work, other intelligent algorithms, such as differential evolution, the 
simulated annealing algorithm, and ant colony optimization, can also be applied to 
select suitable parameters for long short-term memory. Besides, it is also important 
to study solar and wind energy access to power systems. 
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Appendices and nomenclature 


Sets and indices 


ie] all technologies: PV panels, wind turbine, and battery 
teT number of time windows 


Continuous parameters 


Cac total capital cost [CNY] 

Cr total fixed cost [CNY] 

Cri fixed cost of technology i [CNY] 

Cuci unit cost of technology i [CNY/kW, CNY/kWh] 

CRF; capital recovery factor of technology i [-] 

Ct the cell state at time t 

Cr the new candidate information at time t 

D; electric load at time-window t [kWh] 

f: the forget gate at time t 

hi the hidden state at time t 

it the input gate at time t 

Max; the upper limit of the capacity of technology i [kW, 
kWh] 
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Noss the loss efficiency of the battery [-] 

Nch the charging efficiency of the battery [-] 

Ndis the discharging efficiency of the battery [-] 

N; lifetime of technology i [year] 

Or the output gate at time t 

Quri energy output from the technology i, at time-window t 
[kWh] 

X the input information at time t 


Positive variables 


Bi if the technology i is selected for installation, B; = 1; 
otherwise, B; = 0 

Cap; rated power of technology i [kW, kWh] 

S” charging rate of the battery at time-window t [kWh] 

Ya discharging rate of the battery at time-window t 
[kWh] 

S; electricity stored in the battery at time-window t 
[kWh] 
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Chapter 3 


Neutral Network Adaptive Filter 
with Application to Ocean Current 
Energy Estimation 


Hong Son Hoang and Remy Baraille 


Abstract 


This chapter proposes a new approach for the design of an adaptive filter (AF), 
which is based on an artificial neural network (NN) structure for estimating the 
system state. The NNs are now widely used as a technology offering a way to solve 
complex and nonlinear problems such as time-series forecasting, process control, 
parameter state estimation, and fault diagnosis. The proposed NN-based adaptive 
filtering (NNAF) is designed by considering the filtering algorithm as an input- 
output system and two-stage optimization procedure. The first concerns a learning 
process where the weights of the NNAF are estimated to minimize the error 
between the filtered state and the state samples generated by a numerical model. 
The adaptation is carried out next to minimize the mean prediction error (MPE) of 
the system outputs (error between the observations and the system output forecast) 
subject to the coefficients associated with the estimated NN weights. Simulation 
results for different numerical models, especially for state estimation of the chaotic 
Lorenz system as well as for the ocean current at different deep layers which is 
important for renewable energy device placements, are presented to show the 
efficiency of the NNAF. 


Keywords: neural network, Kalman filter, adaptive filter, machine learning, 
minimum prediction error approach, ocean current energy 


1. Introduction 


In this chapter, we develop an NNAF for solving filtering problems, especially 
for dynamic systems of high dimension (an order of 107-108). 

Optimal filtering originated from the astronomical studies at the end of the 
eighteenth century. It is worth of mentioning the early works on the filtering [1, 2], 
where the theory of linear minimum-variance filters is developed independently. 
In Ref. [1], the filter is developed to estimate the desired value of a noisy signal, 
whereas the filter in Ref. [2] minimizes the mean squared error between the 
estimated random process and the desired process. These two works constitute a 
foundation for the celebrated Kalman filter (KF) [3], which processes the measure- 
ments streaming in continuously and updates the estimate of the signal recursively 
upon the arrival of each measurement. 

Under the condition that the system state and observation processes satisfy 
linear equations driven by white Gaussian noises, the KF produces the unbiased 
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estimate with minimum error variance. Since in most applications, the dynamic and 
observation processes are nonlinear, and a linearization technique is applied, which 
results in an extended Kalman filter (EKF). However, the EKF is sometimes far 
from optimal and even divergent. 

The KF has received a great interest for many industrial and engineering fields 
since the 1970s, such as the system control, signal processing, robotics, seismology, 
communication, economics, and finance. 

Direct application of the KF to the data assimilation (DA) in geophysical systems 
such as meteorological or oceanic systems is very difficult, even impossible due to 
high computational effort and memory storage required to deal with very high 
system dimensions, not to say on nonlinearities and uncertainties in statistics of the 
model and observation errors. 

This situation forces researchers to search new, simple but effective algorithms 
to overcome these difficulties. We find here the approaches more appropriate for 
solving DA problems in high dimensional systems (HdSs). For example, the 4D-Var 
algorithm [4] adjusts a forecast, in space and time, to bring it into closer agreement 
with observations. Based on the KF version, the EnKF [5] uses the covariance 
estimated from samples instead of solving the Algebraic Riccati matrix equation. 
The EnKF is related to the particle filter (a particle is the same thing as ensemble 
member). Particle filtering [6] uses a set of particles (also called samples) to repre- 
sent the posterior distribution of some stochastic process given noisy and/or partial 
observations. The state-space model can be nonlinear, and the initial state and noise 
distributions can be non-Gaussian. These approaches are proposed with strategies 
to better control the uncertainty and to make it less non-Gaussian: targeting of 
observations, specific treatment of highly nonlinear degrees of freedom, localiza- 
tion, and model reduction. A review of the application and usefulness of 4D-Var, 
ensemble KF, and particle filters in geosciences is given in Ref. [7]. Among the 
efficient approaches for solving DA in high dimensional setting, we note the adap- 
tive filtering (AF, see Section 2.2). The idea of the AF will be exploited in this 
chapter to improve the neural network filtering (NNF) algorithms. 

To make more clear on why the AF is capable of dealing with the difficulties to 
assimilate data into the HdSs, it is useful to list the main differences between the KF 
and the AF as follows: 


1.The KF is aimed at seeking an optimal estimate in the space of measurable 
Borel functions of the system state. As for the AF, it seeks the estimate, 
optimal in the class of stable filters of a given structure. 


2. Optimality: the KF optimizes the mean squared error between the estimate and 
true system state given observations up to and included the recent moment k. 
The AF minimizes the mean squared error between the observation and the 
forecast of system output. 


3. The parameters to be estimated in the KF are all components of the system 
state, which are time varying. The tuning parameters in the AF are the 
parameters in the filter gain, which are assumed to be time invariant. It is 
mentioned that the tuning parameters in the AF have no physical sense. 


4. Nonlinear matrix equations are to be solved in the KF. 
In this chapter, a new NNAF is proposed and tested for solving DA problems 


with moderate and high dimensional systems. In particular, the problem of esti- 
mating the ocean underwater current by the NNAF will be presented in detail due 
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to its importance in producing energy from the ocean. It is mentioned that there are 
some basic ways to tap the ocean for its energy such as ocean wave, tidal range, tidal 
current, underwater current, and temperature differences in water. All these vari- 
ables can be estimated by the NNAF. Ocean currents are suitable locations for 
deploying energy extraction devices such as turbines. 

At the present, the use of ocean energy remains mainly in the prototype stage. 
Last few years, the prototypes are no longer confined only to laboratories, and some 
companies have set up commercial scale like RITE Project. In 2006, Verdant Power 
installed its first grid-connected Free Flow turbine (see Ref. [8]). 

The modern DA techniques represent a computational challenge, even with the 
use of parallel computing with thousands of processors. Today, many operational 
weather prediction centers use high resolution models, with the number of obser- 
vations (provided by a new satellite generation) growing exponentially. There is an 
increasing need to address new DA algorithms for HdSs and to improve their 
performance. This makes the problem of developing new DA algorithms like NNAF 
a great computational challenge. 

In this chapter, we propose a new version of an artificial NN-based AF as a new 
DA technique, with the objective to construct a high performance forecasting sys- 
tem. Recently, the NN has become a popular and helpful model for classification, 
clustering, pattern recognition, and prediction in many disciplines. For a survey of 
NN applications in the real-world scenario, see Ref. [9]. In Ref. [10], the focus is on 
explaining the concept and evolution of machine learning, and some of the popular 
machine learning algorithms are described and compared. The work of Ref. [11] 
describes the important problems that can arise when applying NN models to 
meteorology and oceanography, especially in large-scale, low-frequency studies, 
and the attempts to resolve these problems. For data assimilation in chaotic 
dynamics using NNs, see Ref. [12]. 

In this chapter, the NNF based on AF developed in Ref. [13] will be constructed 
and implemented to produce the estimation results, which will serve as the reference 
for performance comparison with other NNFs. The need in developing the AF is to 
construct a filter, which is suitable in a new environment. The AF is a powerful device 
for signal processing, control applications, identification, inverse modeling, and pre- 
diction in an uncertain environment. The AF is formally defined as a self-designing 
device with time-varying parameters that are adjusted recursively in order to minimize an 
error function—a distance measurement between the reference or desired signal and 
the output of the adaptive filter. For a review of the AF, see Ref. [14]. 

The idea on the AF for data assimilation in HdS has been first presented in Ref. 
[13]. The originality of this AF approach lies in considering the AF as a filter, which 
minimizes the MPE for the system output, with tuning variables chosen as some 
pertinent parameters of the filter gain. This idea allows to optimize the filter per- 
formance iteratively, in real time, but not on the basis of a batch of observations as 
done in the 4D-Var approach. Moreover, the dimensionless tuning parameters are 
chosen in a way to ensure a stability of the filter. This makes the filter always stable 
during an optimization process and easily controls the values of tuning parameters. 
For a detail presentation of the main characteristics and advantages of the AF in the 
context of data assimilation for HdSs, see Section 2.2. 

Assimilation experiments will be carried out with synthetic conventional data. 
The application of NNF produces a significant reduction in the computational 
effort. The goal of using the NN approach is to obtain a similar quality for analyses 
with better computational performance for prediction process. 

The chapter is organized as follows. In the next section, a brief description of the 
NN and the AF for filtering problems is given. Section 3 presents the scheme of the 
NNF algorithm. The adaptation scheme for the NNFs is developed in Section 4. In 
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Section 5, the simulation experiments on filtering problems for moderate models 
are presented. Here first the one-dimensional system is presented to see clearly the 
idea and structure of the proposed NNAF. The assimilation experiment with the 
widely used Lorenz system, being of chaotic character, is also considered. The 
experiment with the ocean model MICOM (Miami Isopycnic Coordinate Ocean 
Model) for estimation of the oceanic current in different layers by the NNF and 
NNAF is described in Section 6, which demonstrates a possibility of the use of 
NNAF for state estimation in very high dimensional systems. It will be seen that the 
NNAF can improve considerably the performance of the NNF, which is beneficial 
for decision on future planning installation of stations exploiting renewable energy 
resources such as the tides and currents, as well as the waves and the thermal energy 
that resides in the seas. The conclusions are given in Section 7. 


2. Artificial neural network filtering 
2.1 Artificial neural network 


Artificial neural networks (ANNs or NNs) are computing systems that are inspired 
by biological neural networks. Such systems are designed to “learn” by considering 
samples to perform the specified tasks. NNs are based on a collection of connected 
units or nodes called artificial neurons, which model the neurons in a biological brain. 
A standard simple NN scheme is shown in Figure 1. Each connection, like the 
synapses in a biological brain, can transmit a signal to other neurons. An artificial 
neuron that receives a signal then processes it and can signal neurons connected to it. 

Neural networks have been extensively investigated in the context of adaptive 
control and system identification [15], but at the present, not so much in state 
estimation algorithms. 

There are following principal types of NN: (i) feedforward neural network; 

(ii) radial basis function neural network; (iii) recurrent neural network (RNN); 

(iv) convolutional neural network (CNN); and (v) modular neural network. For 
more details, see Ref. [16]. For example, the RNN and CNN are designed to perform 
different tasks: the RNN is useful for learning the structure and syntax of the language, 
whereas the CNN is complex feedforward NN used for image classification and 
recognition because of its high accuracy, and it works by extracting features of image 
(identifying faces, individuals, street signs, and many other aspects of visual data). 

NNs offer the useful properties such as linearity, nonlinearity, learning with 
teacher (modification of the NN weights), and adaptivity (capability to adapt their 
weights to changes in the surrounding environment). A neuron model scheme is 
given in Figure 2. One has 


De = 5 WriXi + br Yp = P(e) () 


i=1 


Figure 1. 
Standard simple NN scheme: input, hidden, and output. 
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Figure 2. 
Model of a neuron. 


where Wz; is the synaptic weight between neuron i (at some / layer) and 
neuron k at the (J + i)”, X1, -Xn are the inputs, ¢(.) is the activation function, bx is 
the bias, and y, is the output signal of the neuron k at the (J + 1)” layer. 

There are different activation functions (AcF) such as Step, Linear, and Sigmoid 
functions. In Figure 3, three AcFs, Step, Linear, and Sigmoid, are displayed. 

One of the most important results concerning a capacity of NN as a powerful 
algorithm for function approximation is expressed in terms of the universal 
approximation theorem, which states that a feedforward network with a single 
hidden layer containing a finite number of neurons can approximate continuous 
functions on compact subsets of R” under mild assumptions on the AcF. The 
theorem thus states that simple neural networks can represent a wide variety of 
functions when given appropriate parameters. One of the first versions of the 
theorem was proved by Cybenko in 1989 for sigmoid activation functions [17]. 

In this chapter, the NN is designed in order to solve assimilation problems. The 
design procedure consists in the following steps: 

Step 1. The choice of NN structure: What is the structure for the NN? This 
concerns the questions on (i) How many layers to be chosen? (ii) What are the 
connections between different nodes in different layers? (iii) What are the values to 
be initialized for the weights? 

Step 2. The choice of variables for the input and output signals. 

Step 3. Formulate an optimization problem for estimating the NN weights. 

Step 4. Solve the formulated optimization problem: What is the optimization 
algorithm to be used? 

All these questions will be addressed in details in Sections 3-6. 

The objective of this chapter is to show that there exists a simple NNF structure, 
which can work efficiently and produce a good estimation results for the system 
states. As will be seen, all the NNFs employed in the experiments will have only one 
hidden layer. As to the AcFs, three types of AcF will be used: Linear, Sigmoid, and 
Step AcFs. 


2.2 Adaptive filter 


The AF for state estimation problem is defined as a filter, minimizing the MPE of 
the system output. One of the main features of the AF is that it is optimal among 


Step AcF Linear AcF Sigmoid AcF 


Figure 3. 
Typical AcFs: step, linear, and sigmoid functions. 
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members of the set of parameterized state-space innovation representations, with 
the vector of pertinent parameters of the gain as a control vector. The assimilation 
problem considered in this chapter is formulated as a standard filtering problem 


x(k +1) = x(k) + w(k) (2) 
z(k +1) = Hx(k+1) +v(k +1),k = 0,1, ... (3) 


where x(k) is the n-dimensional system state at the moment k, and z(k) is the 
p-dimensional observation vector. Under standard conditions, related to the noise 
sequences v(k), w(k), the optimal in mean squared estimate for x(k) can be obtained 
from the system of recursive equations known as the KF. 


x(k +1) = Ox(k) + Ke(k + 1), clk +1) =2(k + 1) — &x(k) (4) 


where K = K(M) is the gain matrix, which is a function of M :=M(k + 1)—the 
error covariance matrix (ECM) for the state prediction error (PE). 


e(k +1) = x(k +1/k) — x(k +1), x(k +1/k) = Ox(k) (5) 


It is mentioned that in the KF, there exists a nonlinear matrix Riccati equation 
for computing M(k + 1). For a system with dimension of order10’-108, under most 
favorable conditions (perfect knowledge of all system parameters and noise statis- 
tics), it is impossible to implement the KF (4) in the most powerful computers due 
to computational cost and insufficient memory. 

In order to exploit the optimal structure of the filter (4), an AF is proposed in 
Ref. [13]. The filter is of the form (4) with the gain K = K(@), which is a function of 
some vector 6 of pertinent parameters. The choice of gain structure and its param- 
eterization (i.e., 0) from the point of view of filter stability is studied in Ref. [18]. In 
Ref. [19], a simple and efficient algorithm is proposed for estimating the prediction 
error covariance matrix (ECM) based on the hypothesis on the separability of 
vertical and horizontal structure of the ECM and its parameterization. This allows to 
estimate optimally the unknown parameters of the structured parameterized ECM 
as well as to approach numerically the solution of the traditional Nearest Kronecker 
Problem in a simple and efficient way. 

The optimality of the AF is understood in the Minimum Mean Squared (MMS) 
sense, which minimizes the objective function: 


J (0) = E[¥(6)], ¥(s) = Ile(ks 8) |? 
J =J(0) > min argo (6) 


with E(.) is the mathematical expectation, and ||.|| is the Euclidean norm. In Eq. (6), 
¢(k; 0) is the PE for the system outputs, which depends on the vector of parameters 0. It 
is seen that if the KF minimizes the MMS filtered error for the system state (in state 
space), the AF minimizes the MMS prediction error for the system outputs in the 
observation space by tuning some pertinent parameters of the filter gain. 


2.3 Simultaneous perturbation stochastic approximation (SPSA) 


The choice of Eqs. (4) and (6) is important in many aspects in order to obtain a 
simple and efficient data assimilation algorithm. The idea lying behind Eqs. (4) and 
(6) is to select the vector 8 € © consisting of some pertinent parameters as control 
variables for minimizing the cost function (6). Here © is the set of admissible values 
of 0, which ensure a filter stability. 
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The solution to the problem (6) can be found iteratively using a stochastic 
optimization (SA) algorithm. 


A(k +1) = O(k) — y(k + 1)Vo¥l[c(k + 1; 0(k))] (7) 


where y(k) is a sequence of positive scalars satisfying some conditions to guar- 
antee a convergence of the estimation procedure. The standard conditions are 


(k)>0, So (k) = 00, $k) <o (8) 


The algorithm of Eqs. (4), (8), (11) is very simple to implement: at the k” 
assimilation instant, we need to integrate the numerical model ® by x(k) to produce 
the forecast x(k + 1/k) and to update 0(k) according to Eq. (7). The algorithm of 
Eq. (7) is completely defined if the gradient of the sample cost function 
Vo¥[ç(k + 1; O(k))] is available. This gradient can be computed using the adjoint 
equation (AE) approach [4]. However, writing the AE code is very complicated and 
hard task for complex physical systems. Moreover, it requires a linearization of the 
system (2) and (4). 

A much less computational burden can be achieved by measuring the sample 
objective function (but not based on a gradient formula): one approximates the 
gradient using the values of the cost function [on the basis of finite difference 
scheme (FDSA)]. This can be done by applying a so-called component-wise 
derivative approximation algorithm. However, for high-dimensional systems, this 
algorithm is inapplicable: for approximation of each partial derivative, we need to 
make two integrations of the model (2). 

Let us look at the class of SPSA algorithms [20, 21]. The algorithm SPSA is of the 
same structure as that of FDSA, with the difference that it perturbs not each 
component of @ separably, but perturbs stochastically and simultaneously all its 
components. Concretely, let Ay := (Agi, «..5 Aen)” be a random vector, A; j,i = 
1, „n are Bernoulli independent identically distributed (iid). The gradient of the 
sample objective function is estimated as 
T 
£ = Veh = (By) By) > (9) 
& = (W(A(k) + cp Ag) — P(O(R) — ceAx)]/(2ce Agi), i = 1, ...,0. (10) 


One sees from Eqs. (9) and (10) that all the directions in the parameter space are 
perturbed at the same time (the numerator is identical in all components). Thus, 
SPSA uses only two (or three) times the direct integration of the model, indepen- 
dently on the dimension of 6, which makes it possible to apply to high dimensional 
optimization problems. No development of the AE code is needed. Generally, SPSA 
converges in the same number of iterations as the FDSA, and it follows approxi- 
mately the steepest descent direction, behaving like the gradient method [20, 21]. In 
fact, SPSA, with a random search direction, does not follow exactly the gradient 
path. In average, though, it tracks the gradient nearly because the gradient approx- 
imation is almost an unbiased estimator of the gradient, as shown in Refs. [20, 21]. 
For SPSA algorithm, the conditions for the sequence of positive cp are 


(ylk) /ck}? < œ (11) 


Me 


foe) 
Ck>0, ) Ck = œ, 
k=1 


> 
ll 


T 


where y(k) satisfies (8). 
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3. NN filtering (NNF) 


Let us return to the input-output system (2)—(3). Instead of the AF, we will 
construct a NNF for estimating the system state x(k) based on the set of observations 
from beginning up to the moment k. According to the NN scheme, we have to define 
the set of inputs and the set of output signals. These two sets are introduced as 


Sr:= {Rn(k/k — 1),2(k)}, So = {2m (k)} (2) 


where S; is the set of inputs, and So is the set of output signals, and the estimate 
Xnn(k/k — 1) = £n (k — 1) is the forecast obtained from the integration of 
Xnn(k — 1) by the numerical model ®. It is mentioned that £„„(k — 1) is the filtered 
estimate produced by the NNF, and symbolically, it is obtained by application of the 
NNF on the basis of the inputs S; subject to the weights W(k), 


Xanlk) = f nn Xnn(k/k — 2),2(k); W (k)] (13) 


The difference between Xp,(k) in Eq. (13) and x(k) in Eq. (4) is that £„„ (k) is the 
estimate produced by the NNF, which is of the NN structure. Let us introduce a 
new model 


Xn (k +1) = Oxm(k) + Wn (k) (14) 


with w,,(k) being a stochastic process, which may be different from w(k). The 
reason is that in practice we are given not exactly the information on the process 
w(k) (in the experiments for simplicity, we assume w,,(k) = 0). The vector x,,(k) 
will serve as samples for learning process to estimate the NN parameters W. For 
illustration, one scheme of NNF structure for the system with two-dimensional 
state and two-dimensional observation vector is given in Figure 4. 

In Figure 4, W/(i,j) denotes the weight between the i node at the /™ layer and 
the j node at the next layer. The algorithm for updating the weights in the NNF 
looks as follows: 

Algorithm 3.1 (Learning algorithm for estimating NNF weights). 

Suppose the set of observations z(k),k = 1, ..., N is available. At the moment 
k =1, let the a priori values W)(k — 1; i,j), Xnn(k/k — 1) be given. Suppose we are 
given the set of observations z(k),k = 1, ...,.N. 

Step 1. Pass the inputs £„„n (k/k — 1), z(k) through the NNF. The NNF produces 
the output (estimate) Xny(k). 

Step 2. Compute the state error enn (k) = Xm(k) — Èn (k). 

Step 3. Using the SPSA algorithm to update the NNF weights according to 


Figure 4. 
Scheme of NNF structure. 
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Figure 5. 
Scheme for Algorithm 3.1. 


W(k) = W(k — 1) —y(k) Vw [W(k — 1)] (15) 
by minimizing the cost function 
Tang (W) = EN (em (W))] = E| lenn(s W)IP| > minarg(W) (16) 


Step 4. If k <N, go to Step 1. Otherwise, go to Step 5. 

Step 5. Stop. 

The scheme for Algorithm 3.1 is shown in Figure 5. 

Comment 3.1. The sample gradient Vw¥[W(k — 1)] is evaluated as shown in 
Section 2.3 with y(k)>0 (jointly with c,>0) satisfying two conditions (8) and (11) 
ensuring a convergence of the algorithm (15) [20]. The sequence {cp} of positive 
scalar values participates in generating simultaneous perturbations to approximate 
Vw¥[W(k —1)]. 

Comment 3.2. The estimates for the system state in (2) are Xn, (k), which are 
produced during running of Algorithm 3.1. To obtain the better estimates on the 
basis of nonadaptive NNF, the NNF should be rerun subject to the weights resulting 
from applying Algorithm 3.1. 

Comment 3.3. If the number of observations is insufficient for well estimating the 
NNF coefficients, it is possible to simulate a longer sequence of pseudo-observations 


Zm(k) = Hxm(k) + vm(k) 


where v,,(k) may be different from v(k) (but chosen as closely as possible to v(k) 
in the statistical sense). The similar experiment will be carried out in this frame- 
work for the ocean model MICOM in Section 6. 


4, Adaptive NNF (NNAF) 


To improve the performance of the NNF, the NNAF is developed in this section. 
The NNAF is of the same structure as that of the NNF, but it is designed to 
minimize the MPE of the system output. The NNAF is based on the second optimi- 
zation procedure with the initial weights resulting from application of Algorithm 3.1. 

For the design of the NNAF, the following objective function is introduced 


Janap (0) =E | (n(W(8))] = E|llem(k: WOI] > minarg(@) a7) 
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The difference between two objective functions (16) and (17) is that if (16) is 
expressed in terms of the (filtered) state estimation error, the error function (17) 
measures the difference between the observation and the prediction of the system 
output. As to the parameters to be adjusted, the following parameterizations are 
possible to be used: 


W(0)=W, +0 (18) 
W(0) = 06W, (19) 
where W, is the estimate obtained from the learning procedure (Algorithm 3.1). 
The initial value for 8 is 6(0) = 0 for the parameterization (18) and 6(0) = 1 for 
(19). It is mentioned that when the structure of the NNAF is closed to the AF in Ref. 
[19], the parameterization (19) is more appropriate. We have now 
Son ( — 1) = f ypl@m( — 1/k — 2),2(k — 1); W(0)] 
More concretely, the NNAF is a result of the following. 
Algorithm 4.1 (Adaptive procedure for improving NNF performance). 
At the initial moment k = 1 suppose we are given Xn,(k/k — 1) and W,(i,j) 
(W, (i,j) are obtained at the end of applying Algorithm 3.1). Suppose the set of 
observations z(k),k = 1, ..., N is available. 
Step 1. Pass the inputs X,,(k/k — 1), z(k) through the NNF initialized by the 


weights W/(i,j). The NNF produces in the output the estimate £n (k). 
Step 2. Compute the error ç(k) = z(k) — Hên (k/k — 1) where. 


Èml(k/k — 1) = Bun (k — 1), Xun (k — 1) = f „n nn (k — 1/k — 2), 2(k — 1); Wil, j)]. 
Step 3. Using the SPSA algorithm to update the weights 
O(k) = O(k — 1) — (Rk) Vora) FlA(R — 1)] (21) 


for seeking the optimal parameters 0 to minimize the cost function 


Ing (0) =| (cm (ki W(0)))] = E| 


cml WOI] > minargo (22) 


Step 4. If k <N, go to Step 1. Otherwise, go to Step 5. 
Step 5. Stop. 


war 


Figure 6. 
Scheme for Algorithm 4.1. 
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Figure 6 displays the scheme of Algorithm 4.1. Compared to Figure 3, here the 
weights are assumed given (as resulting from application of Algorithm 3.1). States of 
the numerical model are no longer participating in the SPSA algorithm; instead, the 
innovation vector is involved in the optimization process. 


5. Numerical experiments for systems of moderate dimensions 
5.1 Experiment 1: One-dimensional system 


We present in this section a very simple numerical experiment with an 
one-dimensional system to illustrate clearly the idea on the NNAF. 

Let us consider the system of Eqs. (2)-(3) subject to ® = 0.99,H =1, Q = 
0.09, R = 0.01. The experiment is carried out on the basis of a NNF, and the 
performance of the NNF is compared with that of the KF and NNAF. The structure 
of the NNF is shown in Figure 7. 

Thus, the NNF has one hidden layer structure. The Linear and Sigmoid AcFs are 
used in this NNF. The weights W/(i,7) are estimated by applying Algorithm 3.1. Here 
the error e(k) = x(k) — Xm(k) is the difference the “wrong” model state governed by 
Xm(k + 1) = ®xn(k),k = 0,1, ..., N and the true system state governed by 
x(k +1) = Ox(k) + w(k),k = 0,1, ..., N. The performance of the NNF is expressed 
in the terms of root mean square (RMS) of the state prediction error. 

Figure 8 shows the performances of three filters—NNF, NNF-FIX, and NNAF. 
The red curve “NNF” is the RMS of state prediction error produced by the NNF 
during the optimization process (Algorithm 3.1). The green curve “NFF-FIX” is 
produced by the NNF with the fix weights obtained at the end of Algorithm 3.1. As 
to the blue curve “NNAF,” it displays the performance of the NNAF by applying 
Algorithm 4.1. The reason for which the NNAF has produced the best performance, 
compared to the NNF and NNF-FIX, is that the NNAF minimizes the MPE for the 
system outputs. It means that learning process does not extract all useful informa- 
tion on the system state contained in the observations. 

Time evolution of the weights, associated with connecting links between the 
nodes in the hidden layer and the node of the output layer, is displayed in Figure 9. 

Figure 10 shows the performances of four filters implemented for estimating the 
system state of the 1d model: NAF (nonadaptive filter), AF developed in Ref. [18], 
NNF (Algorithm 3.1), and NNAF (Algorithm 4.1). The NAF is the filter with the 
constant gain computed on the basis of the stationary solution of the Riccati equa- 
tion associated with the KF. One sees that the NNF behaves better than the NAF at 
the beginning of the assimilation process. Application of adaptation allows to 
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Figure 7. 
Structure of the NNF for the one-dimensional system. 
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Figure 8. 

RMSs of state prediction error resulting from three filters: NNF (Algorithm 3.1, red line), NNF-FIX (green 
line), and NNAF (Algorithm 4.1, blue line). The error produced by the NNF during the optimization process 
is too high compared to that of the NNF-FIX. It (NNF) reaches the same level (of NNF-FIX) at the end of the 


assimilation period. As to the NNAF, the error remains low during all assimilation periods (after about the first 
10 iterations). 
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Figure 9. 

Time evolution of several weights (between the nodes in the second (hidden) layer and the nodes in the output 
layer) in the NNAF during the adaptation process: the weight of the connection link between the first node 
(hidden layer) and the output node (red curve); the weight of the connection link between the second node 


(hidden layer) and the output node (green curve); and the weight of the connection link between the third node 
(hidden layer) and the output node (red curve). 


improve considerably the performances of the two nonadaptive NAF and NNF. In 
general, two adaptive filters AF and NNAF are of the same performance. 


5.2 Experiment 2: Chaotic Lorenz system 
5.2.1 Lorenz equations 
The Lorenz attractor is a chaotic map, noted for its butterfly shape. The map 


shows how the state of a dynamical system evolves over time in a complex, non- 
repeating pattern. 
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Figure 10. 

RMSs of state prediction error resulting from four filters: NAF (nonadaptive filter, red line), AF (green line), 
NNF (blue line), and NNAF (pink line). The NNF, thanks to the learning process, behaves better than NAF. 
The AF developed in Ref. [18] and NNAF (Algorithm 4.1) are of the same performance. 


The attractor itself and the equations from which it is derived were introduced 
by Lorenz [22], who derived it from the simplified equations of convection rolls 
arising in the equations of the atmosphere. 

The equations that govern the Lorenz attractor are: 


dy, /dt = —o(y, — y2)» 
dy,/dt = py, — Vy — Y Y3 4y3/dt = yy — Py (23) 


where o is called the Prandtl number, and p is called the Rayleigh number. All 
o, P, p>0, and usually o = 10, 2 = 8/3, p are varied. The system exhibits chaotic 
behavior for p = 28 but displays knotted periodic orbits for other values of p. 


5.2.2 Numerical model 


In the experiments, the parameters ø, p, p are chosen to have the values 10, 28, 
and 8/3, respectively, for which the “butterfly” attractor exists. 

The numerical model (NM) is obtained by applying the Euler method (first- 
order accurate method) to approximate Eq. (23). Symbolically, we have the NM 


V tera) = Fl y(tk)] Y (Cea) = [v4 (te +1), Vo (th+1)> Y3 (ths 1)] j (24) 


where 6t :=t,41 — t, is the model time step. The observations arrive at the assim- 
ilation instants T and AT; := Tk41 — Tk. 


5.2.3 Numerical model and observation 


The numerical model is derived subject to ôt := 0.005. The true states x(k), 
produced by the NM corrupted by additive noise at the moment Tx subject to 
AT, = 25ôt, are taken as observation z(k),k = 1, ..., No contaminated with the noise 
v, having zero mean and variance o?, with op = 0.1. The true dynamical system, 
corresponding to the transition of the states between two assimilation instants T, 
and T;,41, is expressed as 
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x(k +1) = Flx(k)] + w(k) (25) 


In Eq. (25), w(k) simulates the model error. The sequence {w(k)} is assumed to 
be a white noise having diagonal covariance matrix with variances 2, 12.13, and 
12.13, respectively. The observation system is given by 


z(k +1) = Hx(k +1) +0(k) (26) 


where the observation operator H = [ar nt) hy = (1,0, 0), h2 = (0, 0, 1), that 
is, the first and third components x(k), x3(k) are observed at each time instant k. 
The noise sequence {vx} is white with zero mean and variance R = 21, where I, is 
the unit matrix of dimension n. The initial estimate in all the filters is given by 
&(0) = (5, 10,27)’. 

The true model (TM) state x * (0) is modeled as the solution of Eq. (25) subject to 
the initial state x * (0) = (1.509, 1.531, 25.461)". 

The problem in this experiment is to apply the NNF and its adaptive version to 
estimate the true system state using the observations z(k), k = 1, ..., No, and to 
compare their performances with those produced by the AF. 


5.2.4 NNF structure and numerical results 


The NNF used in this section has the following structure. 

The NNF has only one hidden layer, and its AcFs are chosen as Linear and Step 
AcFs (Figure 3). The structure in Figure 11 is chosen to be as simple as possible. It is 
interesting to ask whether there are more complicated NNFs, which can yield a 
better performance? However, in this chapter, we restrict our self to considering 
only the NNF with the structure in Figure 11. 

The weights in the NNF are first estimated (learning process) using the sequence 
of the erroneous model (EM) states [which is a noise-free version of the TM (true 
model) (25)]. The initial state in the EM is the same as that in the TM. Typical 
evolution of weights during the learning process is shown in Figure 12. It is seen 
that the weights’ estimates are relatively stable, though we cannot say they are 
already stabilized due to an insufficient number of iterations. All the weights are 
initialized with the values 0.5. We have also tested the NNF for the situation with 


Figure 11. 
Structure of the NNF for the Lorenz data assimilation system. 
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Figure 12. 

Typical time evolution of NN weights (between the nodes in the second (hidden) layer and the node in the 
output layer) in the NNAF during the learning process: the weight of the connection link between the fourth 
node (hidden layer) and the first node of the output layer (red curve); the weight of the connection link between 
the fourth node (hidden layer) and the second node of the output layer (green curve); and the weight of the 
connection link between the fourth node (hidden layer) and the third node of the output layer (red curve). One 
notes that the weights do not yet attend a stationary regime. 


the initial weight values equal to 1. The results show that there is no significant 
difference in filter performance between these two situations. 

It is mentioned that the Lorenz system is very sensitive to the initial state 
condition. Different initial states will lead to a variety of solution behavior: some 
solutions are steady; others oscillate between two or more states; and still others 
vary in an irregular manner. 

We generate observations at 160 assimilation instants (4000 model time steps). 
The different filters will be implemented to assimilate the observations. 

Figure 13 shows the performance of three filters—NFF1, NNF2, and NNF3, 


which are initialized with the same initial condition (0) = (5,10, 27)’. It is seen 
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Figure 13. 

RMS of state filtered error resulting from three filters: NNF1 (red line)—the Lorenz model alone (without data 
assimilation), which is running subject to the another initial state (5,10,27); NNF2 (green line)—the NNF 
running with fix initial (connection link) weights (i.e., without learning process); and NNF3 (blue line)—the 
NNF (applying Algorithm 3.1) with the initial weights in NNF2. One sees that the estimation error is too high 
if no observations are assimilated into the model. If the observations are assimilated, with the learning process, 
the error is decreasing as the learning process progresses and becomes lower than that in NNF2 at the end of the 
learning process. 
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Figure 14. 

RMSs of state filtered error resulting from three filters: NNF3 (red line)—the same NNF3 in Figure 13 
(Algorithm 3.1); NNF-FIX (green line)—the NNF subject to the fix coefficients obtained at the end of the 
optimization process in NNF3; and NNAF (blue line)—the NNAF, which is obtained by tuning the weights in 
NNF (initialized by the values obtained at the end of NNF3) to minimize the prediction error for the system 
outputs (Algorithm 4.1) AF (pink line)—the AF (4) and (6) developed in Ref. [19]. 


that without data assimilation, the erroneous initial condition leads to big errors 
between the true states and the states of the erroneous model (red curve, NNF1). By 
simply initializing the weights in the NNF, it is possible to reduce an important 
portion of the errors (green curve, NNF2). At the beginning of the learning process 
(NNF3), the error is relatively high (even higher than that in NNF2), but it quickly 
decreases after about 20 assimilation instants. 

The SPSA algorithm is really capable of well searching optimal weights as seen 
from comparison of performances of the two filters NNF3 and NNF4 (Figure 14). 
Finally, tuning the weights in the NNF to minimize the prediction error (NNAF) of 
the system output yields the best results (among all the filters of the neural network 
structure) as seen from the curve “NNAF” (Algorithm 4.1). To verify the optimality 
of the NNAF, we have also implemented the AF (4) and (5) based on tuning the 
parameters of the filter gain [19]. The performance of the AF is shown by the pink 
curve “AF” in Figure 14, which is almost the same as that of the NNAF. As before, 
all the filters are starting from the same initial condition as used in the NNF1. It is 
interesting to note that the two adaptive filters NNAF and AF are of different 
structures, but they can yield the same performance. If for the NNAF, the choice of 
the NN structure (sets of inputs, outputs; number of layers, initial values of the 
weights to be selected, and learning process) is most important, in the AF the 
question about the choice of stabilizing gain with appropriate parameterization is 
primordial for the success of the AF. 


6. Numerical experiment for high dimensional ocean model MICOM 
The oceans represent a vast source of renewable energy. In general, ocean 


energy can be divided into six types of different origin and characteristics: ocean 
wave, tidal range, tidal current, ocean current, thermal energy, and salinity 
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gradient. Knowledge of these quantities is important for decision making on sizing, 
placement, and optimization of installation of system energy. 

These oceanic quantities are now can be well modeled through an ocean numer- 
ical model, which is constructed according to the ocean configuration of interest. 
The difficulties arising here are that the model solution gives us only the relative 
values of quantities of interest because the model is always far from the reality. 
Therefore, in order to more accurately assess the true values of oceanic quantities, 
the measurement data are required to be inserted into the model to correct the 
model solution. The role of DA algorithms is just to properly (in some sense, in the 
best way) insert data into a numerical model so that to ensure a stability of the 
algorithm and to minimize the estimation error. 

In this section, the experiment with the ocean model MICOM is carried out, 
which is aimed at estimating the ocean current in different depths (vertical layers). 


6.1 Numerical model 


The ocean model used for the DA experiment with the NNF and NNAF is the 
MICOM (Miami Isopycnic Coordinate Ocean Model), which is identical to that 
described in Ref. [18]. The model configuration is a domain situated in the North 
Atlantic from 30°N to 60°N and 80°W to 44°; for the exact model domain and some 
main features of the ocean current produced by the model, see Ref. [19]. 

The MICOM used here has the state x = (h, u,v), where h = h(i,j,k) is a layer 
thickness, and u = u(i,j,k) and v = v(i,j,k) are two velocity components. In the 
model, i = 1, ...,140;7 = 1, ..., 180; and N; = 4 vertical layers, resulting in the state 
vector of dimension 302,400. The model is integrated from the state of rest during 
20 years (ys). During the period of 2 years 19 and 20, every 10 days, the sea surface 
height (SSH) is computed as a linear function of the layer thickness h, which is 
considered as observation in assimilation experiments (totally 72 observations are 
available). To be closer to realistic situations with satellite observations, available 
only at along-track grid points, we assume that the observations are noise-free sea 
surface height (SSH), available only at the set S,of horizontal grid points S, = 
{(ij) i = 1,11, ..., 131, j = 1,11, ..., 171} (but not at all horizontal grid points). 


6.2 Structure of NNF 


We first construct a NNF block for estimating the layer thickness variable h 
using the SSH observations. The Linear and Step AcFs are used in the NNF. The two 
velocity components (u,v)are computed using the geostrophy hypothesis. The 
scheme of the NNF is shown in Figure 15. 


Figure 15. 
Scheme for NNF for data assimilation in MICOM. 
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It is mentioned that as the observations are available only at the set S,, one sub- 
NN block (with sigmoid AcF), named “optimal interpolation,” is designed to inter- 
polate the observations (i.e., to estimate the values of observations), available at S,, 
on all the horizontal grid points of the surface. This operation is very important to 
ensure smooth corrections to all layer thicknesses and to estimate all layer velocities. 

The NNF has one hidden layer as that used in the Lorenz experiment. One note 
from Figure 15 that in the MICOM experiment, the Algorithm 3.1 has been 
implemented not on the basis of the true observations z(k) in Eq. (3) but using the 
other set of “observations” (calculated from x,,(k), see Figure 15) for estimating the 
weights in the NNF (see Comment 3.3). To test the influence of “true samples” and 
“erroneous samples” on the NNF performance, we have implemented two learning 
procedures, one is based on the set of erroneous samples x,«X¢,,, and another on 
the basis of true samples x,,«x*. 


6.3 Numerical results 


In Figure 16, we show the evolution of some NNF weights during the learning 
process on the basis of the set of erroneous samples. 
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Figure 16. 


Typical evolution of estimated weights in different layers during learning process. Initial weights are taken from 
application of the Cooper and Haines’s water column lifting-lowering method [19]. 
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Figure 17. 
Variances of prediction error for the velocity (v component) at the first layer by NNFT (red curve) and NNFE 
(green curve). 
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Figure 18. 
Performance of different filters (in variance) of SSH prediction error (objective function). 


The weights, obtained at the end of learning process, are inserted into two NNFs 
—NNFT and NNFE, where the weights in the NNFT are obtained by minimizing 
the difference between the NNF state and the true system state, and NNFE—with 
the weights resulting from minimizing the difference between the NNF state and 
the state of the erroneous model. Figure 17 displays the variances of the prediction 
error for velocity (v component) at the first layer produced by NNFT (red curve) 
and NNFE (green curve). One sees that learning on the basis of samples of the true 
state yields more optimal weights, especially in avoiding the error pick, but gener- 
ally speaking, only a slight improvement is observed. 

We have tested different filters to see the effect of optimization in finding the 
weights in the NNF and those in the NNAF. The performances of these filters are 
shown in Figure 18. From Figure 18, the error in the NNFO on the basis of the 
initialized weights (red curve) is too high compared to those of the other filters. The 
improvement is observed for the NNF (Algorithm 3.1) during the optimization 
process (green curve). If we run the NNF — FIX with the fix weights, resulting at 
the end of the optimization process in theNNF, a much better performance has been 
produced (blue curve). Finally, the best results are obtained by running the NNAF 
(Algorithm 4.1) (pink curve). 


r 

2 í PONT 
Toor Cece Ten 
on 


4 4 4 
0 10 20 10 a so a 7 tò 
assiemal avian instants 


Figure 19. 
Spatial temporal average RMS of velocity at first two layers for u and v component (by NNF). 
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Time average EMS of filtered error (evelecity) at surface 


Figure 20. 
Spatial temporal average RMS of u velocity at surface produced by NNF, NNF-FIX, and NNAF. 


To have a more complete idea on how work the NNF and NNAF in high 
dimensional setting, in Figure 19, we show the spatial temporal average RMS of 
velocity (u,v) at the first two layers. As expected, estimating the velocity at the 
surface is more difficult than in deeper layers. One remark is that the u-velocity 
component is easier to estimate (i.e. with less error) than the v component. 
Figure 20 gives a picture on performance of different filters. It displays the error 
(for u velocity at the surface) produced by three filters: NNFO, NNF — FIX, and 
NNAF. If in Figure 18, the difference (in SSH forecast error) between NNF — FIX 
and NNAF is less visible for the velocity variable, and the superiority of the NNAF 
over the NNF — FIX is remarkable. This tendency is observed in the comparison of 
velocity estimation errors produced by the three mentioned filters, at different 
layers as well as for the v-velocity component. 


7. Conclusions 


In this chapter, a constructive algorithm offering an attractive approach for the 
design of a NNAF is proposed. In particular, a simple NN structure with one hidden 
layer has proved to be efficient: it can yield a high performance of all the NNAFs 
implemented in the experiments—either with moderate or with high dimensional 
systems. The superiority (in performance) of the NNAF is confirmed even for the 
Lorenz system with chaotic character, which is very popular due to its high sensi- 
tivity to the initial condition. It would emphasize here that in the NNF and NNAF, 
no initial condition is used as a control variable to formulate the optimization 
problems. The experiment with the HdS ocean model MICOM demonstrates a 
possibility to exploit the NNF for state estimation problems for realistic HdSs. 

It is found from the experiments that an appropriate choice of initial values for 
the weights is important for yielding a high performance of the NNF, with a 
moderate number of hidden layers and iteration steps. In the MICOM experiment, 
we have initialized the weights with the values corresponding to the lifting- 
lowering water column method. A more efficient way to initialize the weights (to 
achieve a better performance) is to follow the method described in Ref. [23] based 
on computing the dominant Schur vectors of the system dynamics. It is highly 
probable that the NNF with a larger number of hidden layers and nodes could 


52 


Neutral Network Adaptive Filter with Application to Ocean Current Energy Estimation 
DOI: http://dx.dot.org/10.5772/intechopen.90148 


produce a better performance. The question on how many numbers of layers and 
nodes to be chosen for the NNAF is open and is left for the future study. 

The experimental results confirm that including adaptation as a mechanism for 
minimizing the MPE of the system outputs is an efficient way to improve the NNF 
performance. 

We want to stress that the success of the NNAF is possible partly due to appli- 
cation of SPSA optimization algorithms. These algorithms are relatively ease of use 
for solving difficult multivariate optimization problems and are simple to imple- 
ment and possess a good convergence rate. They require only two objective func- 
tion measurements per iteration regardless of the dimension of the optimization 
problem. These features make SPSA as a powerful tool for solving optimization 
problems with nonlinearities and high dimensionality of complex dynamical 
systems. 

Finally, as the developed NNAF in this chapter is very effective and easy to 
implement, which makes it quite applicable for solving different real-world tasks, in 
particular for solving DA problems in more sophisticated ocean models like the 
HYCOM (HYdrodynamic Coordinates Ocean Modeling). With the HYCOM model, 
it is possible to estimate, for example, thermal energy and salinity gradients—the 
quantities, serving important sources of renewable energy to generate electricity for 
global electrical demand in future. It is mentioned that the model HYCOM is used 
by the French Navy Hydrographic and Oceanographic Service (SHOM) for the 
prediction of ocean currents, to predict variations on the climatic scale, with a 
model of waves and a metric resolution, with very short-lived physics at the littoral 
scale. Implementation of the NNAF represents real challenges for the future opera- 
tional forecasting systems [24]. 
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Abstract 


The accurate parameters extraction is an important step to obtain a robust PV 
outputs forecasting for static or dynamic modes. For these aims, several approaches 
have been proposed for photovoltaic (PV) cell modeling including electrical circuit- 
based model, empirical models, and non-parametrical models. Moreover, numerous 
parameter extraction methods have been introduced in the literature depending on 
the proposed model and the operating mode. These methods can be classified into 
two main approaches including automatic numerical and analytical approaches. 
These approaches are commonly applied in the static mode, whereas they can be 
employed for dynamic parameters extraction. In this chapter, as a first stage, the 
static parameters extraction for both single and double diodes models is exposed 
wherein Genetic Algorithm and outdoor measurements are considered for fixed 
irradiation and temperature. In the second stage, a dynamic parameters extraction 
is carried out using Levenberg-Marquardt algorithm, where 1 day profile outdoor 
measurement is considered. After that, the robustness of the proposed approaches is 
evaluated and the parameters obtained by the static method and that given by the 
dynamic technique are compared. The test is carried out using 3 days with different 
weather conditions profiles. The obtained results show that the parameters extrac- 
tion by dynamic techniques gives satisfactory performances in terms of agreement 
with the real data. 


Keywords: photovoltaic module, static parameters extraction, dynamic 
parameters extraction, empirical model, electrical model 


1. Introduction 


The increasing development of PV technologies brought out their potential to 
provide the energy abundance across the world. Hence, they have been interested 
by several research groups in the purpose to improve their behavior and extend 
their life-time. Meanwhile, giving an accurate forecasting of the PV outputs 
behavior has been always a real issue related to their nonlinearity. Two modes can 
characterize the PV module in terms of modeling, the first one is the static mode 
wherein the obtained model is characterized and validated for fixed weather 
conditions (irradiation and temperature), while the second is the dynamic mode 
where the validation is carried out using variable weather conditions. In this 
context, several models of the PV cell/module have been introduced in the 
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literature based mainly on several approaches including electrical, empirical, and 
non-parametrical modeling. For the non-parametrical models, two approaches 
are introduced namely, Artificial Neural Network (ANN) and Neuron-Fuzzy 
based models. The accuracy and the robustness of these approaches rely strongly 
on the richness of the training dataset in terms of scenarios. Neuron-Fuzzy tech- 
niques has been used to predict cell short-circuit current and open-circuit voltage 
for the static representation [1]. The ANN-based technique is introduced to 
model the PV array power for the embedded systems implementation [2]. This 
technique has been tested for dynamic mode dealing [3]. The ANN technique is 
also used for the prediction of the PV cell/module voltage directed to amorphous 
silicon PV technology wherein the obtained network has been tested by real 
dynamic data [4]. 

In the empirical or the analytical approaches, several models have been proposed 
to estimate accurately both static and dynamic modes of PV cell/modules. An 
analytical model is proposed by [5] based on manufacturer characteristic. This 
model provides acceptable results for both static and dynamic working. Another 
model has been introduced by Sandia National Laboratory [6], which is widely 
employed for PV cell/model forecasting especially for the large-scale arrays. Other 
empirical models have been proposed in order to estimate the PV array power 
under uniform shading [7, 8]. 

In the electrical approaches, two models widely prevail owing to their 
simplicity based on equivalent circuits namely: single (SDM) and double (DDM) 
diodes-based models. These models can deal with both static and dynamic 
modes with an acceptable accuracy. Besides, several enhancements have 
been introduced in these models to minimize parameters number and give 
more simplicity [9, 10]. The expressions of both photo-generation and diode 
saturation currents have been improved to give more accuracy in the 
dynamic working [11]. 

All models possess unknown parameters, which should be identified according 
to the module used in the practice. For this end, numerous approaches have been 
introduced in the literature including analytical and automatic numerical methods. 
The analytical methods usually rely on specific points on the I-V curve and on some 
value given by manufacturer. However, a significant error can be engendered if one 
of more of selected points is incorrect [12]. 

Owing to their adequate results, automatic numerical methods are prevail in PV 
models parameters extraction either through the use of deterministic algorithms 
such as: Newton model modified with Levenberg [13], Levenberg-Marquardt [14], 
Simulated Annealing algorithm (SA) [15], Pattern Search (PS) [16], Nelder-Mead 
Simplex algorithm (NMS) [17], and hybrid Nelder-Mead and modified particle 
swarm optimization [18] or by introducing metaheuristic algorithms such as: 
Genetic Algorithms (GA) [19], Particle Swarm Optimization (PSO) [20-26], 
Cuckoo Search (CS) [27], Artificial Bee Colony (ABC) [26, 28], and Artificial Bee 
Swarm (ABS) [29]. Moreover, other algorithms have been introduced like (FPA) 
[30, 31], hybrid Bee Pollinator Flower Pollination Algorithm (BPFPA) [31, 32], 
Harmony Search (HS) [33], Artificial Fish Swarm Algorithm (AFSA) [34], and 
other algorithms. The majority of the aforementioned algorithms have been applied 
for static parameters extraction. 

Numerical algorithms have been also applied for dynamic parameters 
extraction wherein the identification process is carried out using variables weather 
conditions [35, 36]. 

In this chapter, modeling and parameters extraction of PV cell/module are 
detailed. Where, comparison study among three models by applying static and 
dynamic identification using out-door measurement. 
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2. PV cell/module modeling 


The PV cell presents outputs variation, which depends on weather conditions 
namely, irradiation and temperature. As illustrated in Figure 1, for load variation 
from open circuit to short circuit, the PV cell shows nonlinear characteristic that 
possess a maximum point of power. For an optimal working, the load should be 
adapted at this point. In this section, three PV cell models will be employed and 
improved using automatic parameters extraction namely: the empirical Sandia 
model and both single and double diodes electrical models. 


2.1 Sandia model 


This empirical model given by SANDIA National Laboratories provides rela- 
tively accurate dynamic forecast for PV cell/module by describing the thermal, the 
electrical, and the optical characteristics. Also, this model can be destined for any 
technology and can be adapted with any scale of PV arrays. Furthermore, its sim- 
plicity can qualify it to be used for real-time online prediction. Expressions (1)—(4) 
describe the variation of Impp, Vmpp, and Pmpp, respectively. 


Imp = Imp_stc(CoEe + CE?) (1+ amp(T — Tsrc)) (1) 
Vinp = Vinp_stc + CoN, x 8(T) x In (Ee) + C3N,(5(T) In (Ee))” + Bp (T — T src) 
(2) 
A(T) = n x k x (T + 273.15) (3) 
q 
Pmp = Imp X Vmp (4) 


where, Co_3 are empirical parameters to be identified, Imp src, Vmp_stc are the 
current and the voltage in the maximum power point under standard test condition, 
E, is the effective irradiation, K is the Boltzmann constant, q is the electron charge, 
S(T) is the thermal voltage, Amp and Bmp are, respectively, the current and the 
voltage temperature coefficient [36]. 


“|——41000 Wim2 
_|—B00 Wim2 


z 
È 
Š 


Figure 1. 

PV characteristic for different irradiation and temperature. (a) Current versus voltage for different irradiation; 
(b) Power versus voltage for different irradiation; (c) Current versus voltage for different temperature; (d) 
Power versus voltage for different temperature. 
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2.2 Single diode based model (SDM) 


This physical model is based on the electrical approach illustrated in Figure 2 
wherein the PV cell is composed of: a photo-generation current source and a diode 
while joule losses are represented by two serial and parallel resistors. 

From this electrical representation, expression (5) and (6) can be obtained to 
describe the evolution of both current and voltage. The output current is expressed 
as a sum of the photo-generation current Iph the diode current Ig, and the shunt 
current In. 


I= Ipr — Ia — Lsn (5) 
V+RJI 
lu =1o| exp (C) = (6) 
V+RI 
Ln Ra (7) 


where, K is the constant of Boltzmann, q is the electron charge, T is the cell 
temperature, n is the diode ideality factor, and Ip is the current saturation due to 
diffusion and recombination. 

After the substitution of Eqs. (6) and (7) in (5), the following expression is 
obtained: 


= = 1 
I= Ip — Io [ep ( ORT Ry, (8) 


Rs Rsh and n are parameters to be identified in the static study and can be 
adjusted in the dynamic study. 

Diode saturation current Ip is expressed in Eq. (9) function of the cell tempera- 
ture and the energy band-gap [11]. 


KT 
4X Voc 3 q-Egr 
(exp (Se) = 1) x T,” x exp (- KT, 


Eq. (10) describes the evolution of the energy band-gap E, as function of the cell 


temperature. 
a x T? 
E; = Ego — 10 
g gO (5 + F) ( ) 


where Ego and E,, are the energy band-gap of the silicon at 0°C and at the 
reference temperature T,, respectively, a and £ are constants of the material. 


I x T? x exp (- Z) 
igs 


(9) 
) 


Figure 2. 
Single diode equivalent circuit. 
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The photo-generation current is given by Eq. (11) as a function of the irradiation 
and the cell temperature. 


G 
Lp = —— [Ie +u(T -T, 11 
where G is the input irradiation, T is the cell temperature, Isc is the module 
short-circuit current, and p is the coefficient temperature/short-circuit current 
(given by the manufacturer) [11]. 


Ego, &, B, p, Isc) and Voc are parameters to be identified in the dynamic study. 


2.3 Double diode-based model 


From the electrical representation illustrated in Figure 3, the PV cell can be 
represented by a source of current that represents the photo-generation, two diodes 
and both parallel and serial resistances representing the loss of energy inside the 
cell. 

After applying nodes law, the output current is expressed as sum of: photo- 
generation current lpn, shunt current Isp and the diodes currents Ig; and Iaz 
(Eqs. (12)-(15)). 


I =Ipn — Lar — Ia — Ian (12) 

In = Io È (E) — | (13) 
In = Toa exp (E — | (14) 
In = — (15) 


In which I91_2 are currents saturation of the two diodes that resulted from 
diffusion and recombination, n,_ are ideally factors. 

By substituting Eqs. (13)—(15) in (12), final description of the output current 
versus the voltage is obtained which is expressed in Eq. (16) [11]. 


1 = Ip ~ or exp (LEER) | Toal exp (2089) 1 VERI 


n,KT mKT Ron 
(16) 


ni2 and Rs, Rsa are parameters which will be identified in the static study and 
they can be adjusted in the dynamic study. 

Eqs. (17) and (18) express the evolution saturation currents of the diodes Io, _» 
versus energy band-gap E, and cell temperature [11]. 


Figure 3. 
Double diodes equivalent circuit. 
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3 qx Ei 
Is x T? x exp (- cr) 


Io1 = (17) 
Von E 
(exp (Se) 7 1) x T,? x exp ( See ) 
Is x T? x exp (- Te) 
KT (18) 


Io2 = GX Voc T 3 qX Egr 
( exp (Lee) — 1) x T, x exp (- =] 


The photo-generation current is represented by the same expression of the single 
diode model (Eq. (11)). 

Parameters E,o1-2, 1-2, B1-2, H, Isc and Voc will be identified in the dynamic 
study. 


3. Static parameters extraction of PV module 


For fixed irradiation and temperature, a static parameters extraction will be 
done to extract five parameters in SDM and seven parameters in DDM. A numerical 
stochastic optimization algorithm is used in this identification. This algorithm 
namely, Genetic Algorithm (GA), is employed to minimize the cost function given 
in Eq. (19) which expresses the root mean square error (RMSE) between the 
measured PV module I(v) characteristic and that given by the models. For this and, 
outdoor static measurements have been carried out using the peak measuring 
device tracer (PVPM 2540C), whose characteristics are illustrated in Table 1. This 
device has been programed to provide both I(V) and P(V) curves of 101 samples 
per 1 min. 


RMSE = fey Vm P) (19) 


where I is the simulated current, V,,, is the measured voltage, N is the number of 
sample in I(V) characteristics. The error between the measured and simulated I(V) 
characteristics for the aforementioned models are expressed in Eqs. (20) and (21). 


Vn +R) | Vin + RoI I, (20) 


Vp, eee ee 
FC ) = Iph of exp ( nKT Re 


f (Vins 1, P) = Inn — Tor| exp (C $ zD) 1| Ioa| exp ae = 1| 


nı KT mKRT 
Vin + RI 
=o Bel 
Rsn 
(21) 

Application DC DC Temperature Irradiance Measuring I-V 

voltage current period single curve 

measurement samples 

PV modules 25/50/ 2/5/10/ —40°C to 0-1300 (W/m?) 0.02-2 (s) 101 


and small strings 100/250V 40A  +120°C with Pt1000 (standard-sensor) 


Table 1. 
PVPM2540C characteristics. 


64 


Static and Dynamic Photovoltaic Cell/Module Parameters Identification 
DOI: http://dx.doi.org/10.5772/intechopen.89449 


3.1 Genetic Algorithm 


The Genetic Algorithm (GA) is a stochastic algorithm imitated from the biolog- 
ical genetic process used to find an approximate solution for optimization problems. 
Like in the natural concept, the chromosome is the holder of the genes that the child 
can probably get from his parents. By analogy, these genes represent the variables 
(parameters) of the function to be minimized. Five steps can characterize the GA 
namely, generation of initial population, evaluation of fitness, selection, crossover 
and mutation [37, 38]. 


3.1.1 Initial population 


The process starts by the generation the initial population of N chromosome 
coded in binary. Each vector chromosome is formed of group of parameters in 
which its length M is given in Eq. (22) wherein n is the number of parameters and 
Np is the length of the sub-string (number of bits) of each parameter as shown in 
Figure 4. The length of the integer part given by the vector Conv (Eq. (23)) is used 
to limit the research domain in which, P; (Eq. (24)) is the parameter value in 
decimal code [37, 38]. 


Conv = [2" 271,20 271,20 ga (23) 
P; = [aoay...an,-1] x Conv” (24) 


3.1.2 Fitness 


In this stage, the parameters values that have been randomly generated and 
decoded in decimal base will be substituted in the cost function to be optimized. 
The fitness is the solution of the parameters in the RMSE (x) function calculated in 
Eq. (19). Its value is mathematically expressed in Eq. (23) [37, 38]. 


1 


= TP RMSE) (25) 


Fitness (x) 


3.1.3 Selection 


The chromosomes that will participate as parents to generate a new child are 
chosen in this step. Any chromosome in the generated population can be chosen 
however, the individual that presents a good fitness have a high probability. The 
technique used for the chromosome choice is the roulette wheel illustrated in 
Figure 5, wherein the selection probability P, expressed in Eq. (24), is calculated, 
consists of a cumulative sum of the fitness of each chromosome orderly relative to 
the sum of all fitness. After that, the process generates a random drawing probabil- 
ity Pr. Hence, the first chromosome corresponds to P, < P, is chosen for the next 
steps (crossover and Mutation). 


P, : sescees 2 P, 


Figure 4. 
Chromosome string. 
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E Si Fitness(k) 


P,(i) = E? Fitness (j) 


(26) 


where, k is an integer counter that varies from 1 to the current chromosome, and 
j is an integer counter that varies from 1 to the population size N [37, 38]. 


3.1.4 Crossover 


After selecting the chromosomes, the algorithm gives birth to new children by 
performing a crossover between each two chromosomes. For this end, a drawing 
probability P, is generated and compared with the crossover probability P. (usually 
high probability). Hence, the parents chromosomes that corresponds to P, < P, will 
be chosen for child generation. If not, the same chromosomes are kept. As illus- 
trated in Figure 6, the crossover by point is used wherein the bits after the point 
randomly chosen are swapped [37, 38]. 


> 


u chromosome 1 m chromosome 2 


u chromosome 3 s chromosome 4 


Figure 5. 
Roulette wheel example. 


Parent 1 Child 1 
P/O LE /O OF] L Oo; 1) O;8 1/0/1/1/0j/0/0/1/1/0/0/1 
Parent 2 F Sa Child 2 
y + v 7 + 4 P - x T y = 7 = = 
[o[alolijoliafofalifofo]i] fofifolifejilijifosijfess 
Figure 6. 
Crossover process. 
Child 1 afler crossover Child 1 alicr mulation 
1 itap oafi ao ofa 1oo2 vo uli i olvi 
Child 2 aller erossower Chiki 2 alier mutation 
01i d1 nf @ 1 7 1 OFT on 1@tn1tt}t l 6 
Figure 7. 
Mutation process. 
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3.1.5 Mutation 


In this step, the algorithm introduces a change in some characters of the selected 
chromosomes in order to expand the search space if the initial population does not 


Description SANYO mono-crystalline 
Cell number 96 
Cell type Mono-crystalline 
Cell size 156 x 156 mm 
PV module dimension 1319 x 894 x 35 mm 
Nominal power 180 W 
Open circuit voltage Voc 66.4 V 
Short circuit current Isc 3.65 A 
Voltage Vinpp 54V 
Current Lmpp 3.33 A 
Nominal operating temperature NOCT 45 + 2°C 
Temperature coefficient (Pmax) —0.33%/°C 
Temperature coefficient (Isc) 1.10 mA/°C 
Temperature coefficient (Voc) —0.173 V/°C 
Table 2. 


PV modules characteristics. 


PV module software 


Figure 8. 
(a) Pictures for I-V experimental platform and (b) synoptic scheme for I-V measurement. 
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fall in the optimal solution. In the binary coding, the selected bit will change from 1 
to 0 and vice versa as described in Figure 7. The mutation has low probability Pm in 
which, it will be affected for characters that correspond to P, < P,, in which, P, is 
the drawing probability (randomly generated) [37, 38]. 

Our system is formed of mono-crystalline PV module SANYO technology with 
the characteristics listed in Table 2, peak measuring tracer and the necessary sen- 
sors. The experimental platform is illustrated in both Figure 8a and b. 

After running of GA for 1000 cycles with the parameters listed in Table 3 for 
both SDM and DDM using outdoor measurement of the systems, wherein extracted 
parameters are listed in Table 4. Figure 9a and b illustrate the agreement between 
the measured and simulated I(V) and P(V) characteristics for SDM model whose 


GA parameters Value 

Number of cycle 1000 

Population length 500 

Crossover probability 0.7 

Mutation probability 0.2 
Table 3. 


The GA parameters. 


The electrical parameter Iph [A] Ip [A] n Rs [Q] Rsh [KQ] 
The identified value 3.0195 49591e—005 1.874*96 0.3273 8.1514 
Table 4. 


The obtained parameters for SDM. 


(a) {b} 
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Figure 9. 
Matching between measured and simulated characteristics for SDM. (a) I(V) and (b) P(V). 


The electrical Iph Io; [A] Io2 [A] ny n Rs Rsh 
parameter [A] [Q] [K9] 


The identified value 3.0289 6.1035e—005 3.8147-006 1.3658*96 1.9179*96 0.1017 5.992 


Table 5. 
The obtained parameters for DDM. 
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Current [A] 
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35 35 
Voltage [V] Voltage [V] 
Figure 10. 

Matching between measured and simulated characteristics for DDM. (a) I(V) and (b) P(V). 


E AE 

NS NL ZE PL PS 
NS PS PS PL PS PS 
NL PS PL PL PL PS 
ZE NL NL ZE PL PL 
PL NS NL NL NL NS 
PS NS NS NL NS NS 

Table 6. 


The inference table. 


obtained parameters are summarized in Table 5. Figure 10a and b show the agree- 
ment between the measured and simulated I(V) and P(V) characteristics for DDM 
model whose extracted parameters are summarized in Table 6. 

Some parameters will be identified again in the dynamic study including the 
parameters involved in Io and I, equations, while the remaining will be adjusted to 
give more accuracy under variable weather conditions. 


4. Dynamic parameters extraction of PV module 


In this section, dynamic parameters identification will be described wherein the 
process is done by using 1 day profile of measurement. This allows to improve the 
nominal values given by the manufacturer, which can cause a significant error due 
to operating conditions and the consumed lifetime. Moreover, parameters obtained 
by static method can be adjusted by dynamic identification. For this end, automatic 
parameters adjustment using Levenberg-Marquardt optimization algorithm is 
employed. 

As illustrated in Figure 11, the main idea is to take both PV module model and 
the MPPT as a single system with three outputs namely, Impp, Vmpp, and Pmpp- These 
outputs will be compared with 1 day profile of outdoor measurements. The process 
consists in minimizing the error between the model outputs and the real data. The 
whole system has been implemented in Matlab/Simulink tool. 
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For Sandia model, the process is carried out without the use of MPPT consider- 
ing that this model has been established to the dynamic forecasting. 


4.1 The MPPT used 


An Accurate fuzzy logic MPPT algorithm is employed in our system (for SDM 
and DDM) in order to get satisfactory results in terms of precision and accuracy. 


Papp V myy Tape 


Lrradiation / Temperature 


Extracted parameters 


Figure 11. 
Dynamic parameters extraction process. 


\/ 
| AXA | 


Fuzzy logic algorithm 
Output Duty cycle 
Figure 12. 
Fuzzy logic algorithm steps. 
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The algorithm is used to control a DC/DC boost converter for the purpose to keep 
the PV module working at the maximum point of power. Mamdani inference model 
is used with two inputs namely, the error E and the variation of the error AE. The 
calculation of these attributes is expressed in Eqs. (27)—(30). 


APp (n) = Pov (n) — Pp (n — 1) (27) 
AVp(n) = Vpn) — Veo(n — 1) (28) 


Figure 13. 
(a) Error membership; (b) variation of error membership; and (c) duty cycle membership. 
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E(n) = an (29) 
AE(n) = E(n) — E(n — 1) (30) 


Three steps can characterize the fuzzy algorithm; the first one is the 
fuzzification process that consists on the conversion of the numerical inputs 
values (E and AE) into linguistic values by the substitution in the membership 
functions. The second step is the inference process, which is considered as the main 
stage in the fuzzy algorithm wherein the relation between the inputs and the 
output is done. The third step is the defuzzification where the process converts the 
linguistic decision into numerical output. Figure 12 describes briefly the fuzzy 
processing steps [39]. 


SDM parameters Boundaries DDM parameters Boundaries Sandia parameters Boundaries 


n [0, 2]*96 nı Ny [0, 2]*96 Co [0, 2] 
R, [Q] 0,1] R, [Q] [0, 1] Ch [-1, 1] 
Ry, [Q] [0, 104] Rg, [Q] [0, 10*] G [—10, 10] 
Ego 0,1 Eyo1-2 [0, 2] C3 [-10, 50] 
A 0,1 a2 [0, 1] Amp [°C] [0, 1] 
B [0, 104] Pra [0, 10*] Bmp [V/°C] [-1, 0] 
u 0,1 u [0, 1] n 
Isc [3, 3.7] Ise [3, 3.7] 
Voc [60, 66.8] Voc [60, 66.8] 

Table 7. 


Lower and upper boundaries selected for each model parameters. 


SDM Values DDM Values Sandia Values 
parameters parameters parameters 
n 105.73/96 n1, N2 90.73/96; 73.39/ Co 1.058 
96 

R, [Q] 0.82495 R, [Q] 0.3219 Cy 0.020 
Ren [Q] 8.371 x 10° Ron [Q] 4.9664 x 10° C2 —0.341 
Ego [ev] 1.4525 E,o1-2 [ev] 1.649; 1.31 C3 —9.997 
a 6.56 x 1074 a2 0.0018; 0.0132 imp [°C] 2.53 x 1071 
B 126.11 Pr2 694.84; 1020.76 Bvmp [VI°C] —0.203 
H 0.0121 H 0.0112 n 1.221 
T 3.671 Isc 3.629 
Vo 66.208 Vor 65.527 

Table 8. 


Dynamic extracted parameters. 
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For both inputs and output, five trapezoidal and triangular membership func- 
tions have been employed namely: NS (negative strong), NL (negative low), ZE 
(zero), PL (positive low) and PS (positive strong). The center of gravity based 
method is used for the defuzzification to provide the control duty cycle after 
applying the Mamdani inference model given in Table 6. Figure 13a, b and d 
describes the used membership functions [39]. 


4.2 Simulation study 


The Levenberg-Marquardt algorithm is implemented using 1 day profile of 
outdoor real measurement of dynamic PV outputs (Pmpp, Impp and Vmpp). The 
process consists in minimizing the error between simulated outputs of both 
SDM and DDM and 8 h of real data (09:00 am—05:00 pm). The peak measuring 
device tracer (PVPM 2540C) has been programmed to provide 1 sample per 
minute. Table 7 lists the lower and upper limits search of the extracted 
parameters. The extracted parameters using the dynamic method are summarized 
in Table 8. The inputs measurement of the irradiation and the temperature are 
illustrated in Figure 14a and b, respectively. Satisfactory results have been 
obtained in terms of matching between the real data and the simulated outputs 
Pinpp» Impp and Vinpp for SDM, DDM and Sandia as shown in Figure 15a-c, 
respectively. 
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Figure 14. 
Eight hours profile (09:00 am-05:00 pm). (a) Irradiation and (b) temperature. 
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Figure 15. 
Eight hours profile (09:00 am—05:00 pm) of: measured output VS obtained output for SANYO mono- 
crystalline. (a) Measured Ppp vs obtained Pmpp from: SDM, DDM and Sandia, (b) measured Impp vs 


obtained Impp from: SDM, DDM and Sandia and (c) measured Vinpp VS obtained V mpp from: SDM, DDM 
and Sandia. 


5. Experimental validation 


In this section, a validation with an unseen data is carried out to test and 
compare the effectiveness of the proposed enhancement. The three developed 
models will be compared with real measurement profile (09:00 am—05:00 pm) of 
irradiation and temperature for different weather conditions. Wherein, the SDM 
and DDM models using the developed parameters are compared with the former 
nominal parameters listed in Table 9, Sandia model and the real data. Three differ- 
ent skies of real measurement have been used for this validation namely, clear day, 


74 


Static and Dynamic Photovoltaic Cell/Module Parameters Identification 
DOI: http://dx.doi.org/10.5772/intechopen.89449 


Parameters Value 
Ego 1.16 eV 

a 4.73 x 1074 
B 1000 

H 1.10 mA/°C 
Tye 3.65 A 
Voc 66.4 V 

Table 9. 


Nominal parameters. 


odlom Wom nom iam Hoop Poop ~~ pC 000 pm 
Tene 


Figure 16. 

Eight hours profile (09:00 am-05:00 pm) of: measured Pippy vs obtained Pp, from: SDM, SDM with nominal 
parameters, DDM, DDM with nominal parameters and Sandia. (a) clear day, (b) semi-cloudy day and (c) 
cloudy day. 
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Figure 17. 

Eight hours profile (09:00 am-05:00 pm) of: measured Impp vs obtained Impp from: SDM, SDM with nominal 
parameters, DDM, DDM with nominal parameters and Sandia. (a) clear day, (b) semi-cloudy day and (c) 
cloudy day. 


semi-cloudy day and cloudy day. The matching in the power (P,p,) between and 
the real data and SDM and DDM with nominal parameters, SDM and DDM with the 
new parameters and Sandia model is illustrated in Figure 16a—c for clear day, semi- 
cloudy day and cloudy day respectively. Besides, the agreement in the voltage and 
the current (Vinpp and Impp) for these models with the real data is shown in 

Figure 17a-c and Figure 18a-c, respectively. 
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Figure 18. 

Eight hours profile (09:00 am-05:00 pm) of: measured V npp vs obtained V mpp from: SDM, SDM with nominal 
parameters, DDM, DDM with nominal parameters and Sandia. (a) clear day, (b) semi-cloudy day and (c) 
cloudy day. 


It is clearly found that an improved agreement has been shown by models with 
new parameters compared to that given by the nominal parameters and the static 
method. 

For more clarity, the hourly power efficiency given by the presented models and 
the real data has been calculated. It consists on the average of the power during 1h 
versus the optimal PV module power (Eq. (31)) [40]. 

— oeerPmpe(k)/N 


Hourly efficiency = “= (31) 
d 4 Ppv_optimal 
_optim 


in which, N is the number samples per hour. 
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Eight hours profile (09:00 am—05:00 pm) of: measured hourly efficiency vs obtained hourly efficiency from: 
SDM, SDM with nominal parameters, DDM, DDM with nominal parameters and Sandia. (a) Clear day, (b) 
semi-cloudy day, and (c) cloudy day. 


Figure 19a—c show the bar-graph of the hourly power efficiency of the proposed 
models for the three weather conditions, namely the clear day, semi-cloudy day and 
cloudy day, respectively. The enhanced models present higher hourly power effi- 
ciency versus models with the former parameters and those given by the static 
technique. Furthermore, root mean square error (RMS) and the mean absolute 
error (MAE) between the real data and the studied models are calculated by 
Eqs. (19) and (32) to show the enhancement of the proposed method. 


N . 
|Simulated output — Measured output| 
MAE = 
N 2 Measured output 


i=1 


(32) 
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Day Weather Error Current Voltage 
[%] 


SDM SDMnp DDM DDMnp Sandia SDM SDMnp DDM DDMnp Sandia 


22/01/ Clear RMS 0.72 1.12 0.69 0.93 0.68 1.48 2.21 1.63 2.34 1.49 
2018 day 


MAE 4.31 5.05 4.24 5.12 4.38 1.27 4.07 1.77 3.70 1.17 


29/01/ Semi- RMS 0.22 0.45 0.24 0.32 0.23 3.30 9.10 6.53 7.87 1.13 
2018 cloudy 


MAE 3.14 5.40 2.96 4.35 2.85 1.50 5.34 2.37 4.76 0.57 


28/01/ Cloudy RMS 0.43 0.48 0.47 0.31 0.52 9.05 8.77 12.2 12.3 1.83 


2018 day 
MAE 4.38 6.62 5.02 3.93 5.60 2.06 4.63 4.05 6.18 0.77 
Day Weather Error [%] Power 
SDM SDMnp DDM DDMnp Sandia 
22/01/2018 Clear day RMS 3.65 8.49 3.45 4.34 3.79 
MAE 4.24 6.79 3.74 4.44 4.51 
29/01/2018 Semi-cloudy RMS 2.48 4.30 2.48 4.30 0.82 
MAE 5.14 11.7 4.58 10.3 2.33 
28/01/2018 Cloudy day RMS 5.35 4.82 5.51 5.24 2.79 
MAE 1.37 1.56 1.55 1.57 5.63 
Table 10. 


Calculated RMS (%) and MAE (%). 


where N is the number of samples [41]. 

The aforementioned results show clearly that the extracted parameters of the PV 
module using dynamic techniques present more accuracy compared with the static 
method and the parameters given by the manufacturer. Indeed, the parameters 
obtained by the static method are clearly improved for variable weather conditions 
(irradiation and temperature), which is confirmed using different skies. Table 10 
summarizes the calculated RMS and MAE errors values which show obviously that 
the developed models present advantages comparing with real outdoor data of 
different weather conditions. 


6. Conclusion 


In this chapter, both dynamic and static parameters identification methods have 
been highlighted and compared with real measurement. The SDM and DDM nom- 
inal parameters involved in Io and Ipn equations have been developed by dynamic 
method. This improved result has been compared with that given by the static 
technique and Sandia model versus out-door real data for different skies (clear day, 
semi-cloudy day and cloudy day). It was found that SDM and DDM based on the 
parameters extracted by dynamic method give satisfactory accuracy, which is con- 
firmed by some calculated indicator such as: the hourly efficiency and both root 
mean square error (RMS) and the mean absolute error (MAE). This allows to solve 
modeling problems of PV module that apply for several applications such as fault 
detection. 
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Chapter 5 


Performance Analysis of a20 MW 
Grid-Connected Photovoltaic 
Installation in Adrar, South of 
Algeria 


Nouar Aoun 


Abstract 


This chapter presents the performance of a20 MWp grid-connected PV system 
installed in a harsh environment, Adrar in the South of Algeria. The results were 
monitored over a period of 1 year, from January 2018 to December 2018. The PV sys- 
tem assessment includes final and reference yields, system efficiency, performance 
ratio, capacity factor, and total system losses. The total yearly energy delivered to 
the grid was 35892.22 MWh, and the monthly average reference and final yields of 
the system are in the range 5.92-8.1 (h/day) and 4.39-5.56 (h/day), respectively. 
Furthermore, the annual average daily PV system efficiency, performance ratio, 
capacity factor, and total losses were 10.82%, 71.71%, 20.76%, and 2.04 h/day, 
respectively. 


Keywords: photovoltaic, grid-connected, final yield, PV system efficiency, 
performance ratio, total losses 


1. Introduction 


In recent years, Algeria has shown great interest in new and renewable energies. 
Therefore, on February 19, 2015, the first PV plant was commissioned in Djanet 
(stand-alone system) in the South-East of Algeria with a total power of 3 MW. In 
the same year, five PV plants with a total capacity of 48 MW were integrated into 
the national grid. Most of the photovoltaic plants were installed in the south of 
Algeria where the insulation time exceeds 3500 h/year. Evaluating the performance 
of the installed photovoltaic systems is indispensable; it gives researchers the 
required information regarding the technology (PV, inverter) used and its suitabil- 
ity for the installation site. 

In the literature, many studies have investigated, evaluated, and analyzed PV 
grid-connected systems in different locations, such as India [1], Italy [2], Ireland 
[3], Ghana [4], Iran [5], Brazil [6], Mauritania [7], Korea [8], Greece [9], and 
Morocco [10]. 

In India (South India), the performance of 5 MWp photovoltaic grid-connected 
systems was presented in [11], by using the RETScreen software for the valida- 
tion of the performance of the PV plant in real time. It was found that the annual 
average energy generated by the system is 24116.61 kWh/day, and observed that the 
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annual average daily array yield is of 5.46 h/day, reference yield of 5.128 h/day, final 
yield of 4.810 h/day, module efficiency of 6.08%, inverter efficiency of 88.20%, 
and system efficiency of 5.08%. Moreover, a performance study of a190 kWp 
grid-interactive photovoltaic plant was conducted by Sharma and Chandel [12], 
who found that the final yield, reference yield, and performance ratio vary in 

the range 1.45-2.84 h/day, 2.29-3.53 h/day, and 55-83%, respectively. In [13], the 
authors presented a grid-connected photovoltaic system installed on the rooftop of 
the Siksha ‘O’ Anusandhan University. The results showed the total energy gener- 
ated during the period test was 14.960 MWh. The array yield, final yield, module 
efficiency, inverter efficiency, and performance ratio of the system were found to 
be 4.09 h/day, 3.67 h/day, and 13.42, 89.83, and 78%, respectively. In a different 
study, Doolla and Banerjee [14] analyzed the variation in the capacity factor of PV 
grid-connected systems around India. They found that the capacity factor varies 
between 16 and 20%, and lower values were observed in northeastern states of the 
country due to poor solar irradiance throughout the year. 

In other studies, Yaacob et al. [15] evaluated three different PV systems in 
Malaysia, namely concentrating PV system (CPV), PV sun tracking flat (TF), and 
fixed flat PV system (FF). It was concluded that the most suitable system is the 
tracking flat PV system with an average daily generation, system efficiency, power 
efficiency, average daily yield, and capacity factor of 4.7 kWh, 11%, 85%, 2.3 kWh/ 
kWp, and 32%, respectively. Mpholo et al. [16] evaluated the performance of a 
281 kWp grid-connected photovoltaic system in Lesotho, Roma (Italy). The results 
show a satisfactory performance, with a weighted performance ratio of 70%; a 
monthly yield for the array (YA), reference (YR), and final (YF) in the ranges 
(71.7-168.0), (161.6-199.3), and (64.4-151.4) kWh/kWp, respectively. 

In Algeria, although different photovoltaic systems were installed around 
the country, they were poorly investigated. Only some efforts were found in the 
literature to evaluate the performance of mini-systems connected to the power 
grid [17, 18]. 

For the first time, this chapter attempts to predict the performance of a grid- 
connected system installed in a desert area in the Adrar region, South of Algeria, 
and a performance evaluation of a20 MW PV system over a 1-year period (January 
2018 to December 2018) is presented. The IEA 61724 standard (International 
Energy agency) guidelines were used to analyze the performance of the system. 
The PV system assessment includes final and reference yields, system efficiency, 
performance ratio, capacity factor, and total system losses. 


2. PV plant description 


The 20 MW grid-connected PV plant of Adrar comprises 20 subfields of 1 MWc 
each and divided into two subsystems. In total, the numbers of modules per subfield 
are 4092. The modules are mounted at an angle of 28° and oriented southward. The 
type of module is polycrystalline with a maximum power of 245 Wp and total module 
area of 132,538,775 m°. The manufacture characteristic of the PV module is shown 
in Table 1. The station was put into service on October 28, 2015, and was expected to 
produce annually 36,414 MWh, with an economy of about 10,080 m*/year of natural 
gas. However, the results presented in this study consist of the measured data over a 
period of 1 year (January 2018 to December 2018). A block scheme of the considered 
photovoltaic system is illustrated in Figure 1. 

The in-plane global solar radiation, wind speed, ambient temperature, and the 
output AC power were measured and recorded after every 15-min interval. 


86 


Performance Analysis of a 20 MW Grid-Connected Photovoltaic Installation in Adrar... 
DOI: http://dx.doi.org/10.5772/intechopen.89511 


Parameters Specification 
Manufacturer Yingli Green Energy Holding Co. Ltd. 
Cell type Polycristallin 
PV Model YL245P-29P 
Maximum power, P, (W) 245 
Maximum power voltage, V, (V) 29.6 
Maximum power current, I, (A) 8.28 
Open circuit voltage, V» (V) 375 
Short circuit current, Is (A) 8.83 
Number of cells 60 
Module dimensions (mm) 1640*990*35 
Module efficiency (%) 15.1 
Maximum system voltage (V) 1000 
Temperature coefficient of P,, app (%/°C) —0.42 
Temperature coefficient of Le, are (Y/°C) 0.05 
Temperature coefficient of Vye, avoc (Y%/°C) —0.32 
Table 1. 


Photovoltaic module specifications. 


Subsystem 1 


$ vvvy } 


Subfield 1 Subfield 10 


vvvy 


Subsystem 2 


Power transformer 


Power transformer 


Power grid 


Figure 1. 
Schematic representation of the considered PV system. 


We observe a few days missing in the months of: January (2-day data is miss- 


ing for both subfields), July (5-day data is missed for the second subsystem), and 
October (3 days are missed concerning the two subsystem). 
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3. Performance analysis 


Equations (1)—(9) are applied in accordance with the IEC 61724 standard to 
investigate the performance of the grid-connected PV system. 


3.1 Final yield (Yp) 


Final yield (Yp) is the total AC energy generated by the PV system over the 
monitored period divided by the rated output power of the PV system: 


E 
Yra = ps (1) 
PV,rated 
1 N 
Vem = EDV pa (2) 
1 


The total daily and monthly AC energy generated by the system was obtained as: 


24 

Eaca = T° 2 Pp (3) 
N, 

Eacym = dEaca (4) 


where N is the number of days per month, t is the recording interval (t = 15mn), 
P, is the AC power in kW, and Eac is expressed in kWh. The terms d and m indicate 
the daily and monthly period. 


3.2 Reference yield (YR) 


Reference yield is the ratio of total in-plane solar irradiation H, (kW/ m7’) and the 
reference irradiance Herc (=1 kW/m’). It can be calculated as follows: 


_T:2G, 
YRa = Gsrc (5) 
The monthly average daily reference yield is obtained as: 
Yam= 2 SY 6 
R,m = we R,d ( ) 


3.3 Performance ratio (PR) 


The performance ratio is a dimensionless quantity which provides important 
information concerning the system losses (modules, inverters, cables, weather con- 
ditions, loss due to non-STC temperature) in the PV systems. This quantity defined 
as the ratio of the final yield (Yp) and the reference yield (Ypg). It indicates the 
percentage of the real energy supplied by the system [19, 20] and can be expressed 
as follows: 


_ Yr 
PR= (7) 
3.4 Total energy losses (Lr) 


The total energy loss (T,) of the system is obtained from the difference between 
the reference yield Yg and the final yield Yp. The total losses represent the PV losses 
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due to irradiance level and array temperature and quality module, ohmic wiring, 
mismatch, and total inverter losses [21]: 


Lr=Yp-Yp (8) 
3.5 Capacity factor (CF) 
The capacity factor (CF) is the ratio of the actual energy output of the PV 
system and the PV system energy generates operating at full rated power. The yearly 


and monthly capacity factors were calculated by Eqs. (9) and (10), respectively: 


Eac 
Poen a 
G K Pry rated x 8760 (9) 


Eacun 
CF = Pir aa AEN (m 
where the term (y) indicate the yearly period. 
3.6 System efficiency 


The monthly system efficiency is given as: 


Eacum 
Nsys,m = Hem xA, (11) 


where 
N/ 24 
Aim = AE G) (12) 
T\ T 


where H, m monthly total in-plane insolation (kWh/m’). 


4. Results and discussion 


The monthly average daily ambient temperature and wind speed over the moni- 
tored period is shown in Figure 2. The monthly average daily ambient temperature 
varies from a minimum of 14.05°C in December to a maximum of 40.8°C in July. 
The months of May, June, July, August, and September (summer months) show 
higher ambient temperature compared to the other months. The monthly variation 
of daily average wind speed is from a minimum of 3.59 m/s in October to a maxi- 
mum of 5.14 m/s in April. Figure 2 shows that the wind speed over the monitored 
period (during all the months) is higher. 

The monthly total AC energy generated to the grid and in-plane irradiation 
over the monitored period is shown in Figure 3. The total in-plane irradiation 
varies from a minimum of 171.7 kWh/m? in January to a maximum of 251.2 kWh/ 
m’ in May. The highest values were recorded in the months of March, April, May, 
June, July, and August (Spring and Summer). The monthly total variation of the 
AC energy generated to the grid ranged from 2725.15 MWh in July to 34475 MWh 
in March. As shown in Figure 3, there is a correlation between the total AC energy 
and the total in-plane irradiation. As the total irradiation increases, the total energy 
likewise increases. In July, a 5-day loss of energy value data for one of the subsystem 
was observed. This explains the low energy value in the month of July despite high 
irradiation. From Figures 2 and 3, it is evidently clear that increases in ambient 
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Figure 2. 
Monthly average daily ambient temperature and wind speed over the monitored period. 
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Figure 3. 
Monthly total AC energy generated to the grid and in-plane irradiation over the monitored period. 


temperature have a negative effect on the energy produced by the system due to the 
increased PV modules temperature, although the wind flow on the photovoltaic 
modules contributes to reducing the PV cell temperature [7, 19]. Generally, the irra- 
diation had a greater impact on energy generation than ambient temperature [22]. 
Figure 4 shows the average monthly performance ratio and efficiency of the 
PV system over the monitored period. The results show that the performance ratio 
varies from a min. of 55.29% in July to a maximum of 81.96% in January, with an 
annual average of 71.71%. Similarly, the efficiency varies from 8.34% in July to 
12.37% in January. The performance ratio and efficiency refer to the level of losses 
in the photovoltaic solar system. There are three main reasons for losses of per- 
formance are dust, ambient temperature, and solar radiation. In [22], the authors 
indicated that ambient temperature had a greater impact on the performance ratio 
and efficiency than irradiation, and high values for both parameters were observed 
during days with low ambient temperature. 
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Figure 5 shows the monthly and annual capacity factor of the PV system. The 
monthly capacity factor can be seen to vary between 18.31% in July and 23.17% 
in March. The annual capacity factor is 20.77%, with an overall monthly average 
of 20.76%. The capacity factor of other PV plants based on existing literature: in 
India [1], it ranged from 15.4 to 20%, Oman from 13 to 20% [23] and 21% [24], and 
Morocco from 6.55 to 21.42% in [10] and for a polycrystalline PV technologies the 
CF ranged from 20.014 to 23.799% [25]. 

Figure 6 shows the monthly daily average variation of the final yield, reference 
yield, and total energy losses. As shown in Figures 3 and 6, the reference yield is 
directly proportional to the irradiation, respectively; they vary from a minimum 
of 5.92 (h/day) and 171.69 kWh/m? in January to a maximum of 8.1 (h/day) and 
251.19 kWh/m? in May. 

In the present study, the annual monthly average daily final yield of the PV plant 
is 4.98 (h/day). The reason for this high value is the high irradiation and the extended 
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Figure 4. 
Performance ratio and efficiency of the PV system over the monitored period. 


Ga PV Plant — — -Yearly Plant 


20 --m@ oe pe ua i 

| | | | | | 

0- T ; T ; r 

Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec 
Months 


Capacity factor, CF (%) 
S 


Un 


Figure 5. 
Monthly and annual capacity factor over the monitored period. 
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Monthly average energy yield and energy losses over the monitored period. 


sunshine duration in Adrar [10]. This value is close to other final yields of different 
PV systems found in the literature such as that of Oman (Sohar) 5.1 (h/day) [24], 
and higher than in Kuwait 4.5 (h/day) [26], in Oman (Muscat) 4.1 (h/day) [23] and 
in Morocco 4.45 (h/day) [10]. The monthly daily average yields (final and reference) 
increase during the months from May to August due to higher irradiation and more 
sun hours. 

The total energy losses varied between 1.06 h in January and 3.55 h in July. The 
higher total energy losses were observed in the hot months, i.e., May to August. The 
reasons for the losses are generally due to dust and ambient temperature (heating of 
the photovoltaic cells which reduce the power provided by the photovoltaic modules). 
And in these conditions, periods of downtime of the system were observed (2 days in 
January, 5 days in July and in October stopped for 3 days with 1-day data incomplete). 


5. Conclusion 


The present chapter investigates the performance analysis assessment of 
ground-mounted large-scale grid-connected PV plant in Adrar, South of Algeria. 
The PV system is made up of 20 MW of polycrystalline silicon modules produced 
by Yingli Green Energy Holding Co. Ltd. and monitored between January 2018 and 
December 2018. 

The PV system assessment includes final and reference yields, system efficiency, 
performance ratio, capacity factor, and total system losses. The main findings are 
the followings: 


e The monthly total variation of the AC energy generated to the grid ranged 
from 2725.15 MWh in July to 34475 MWh in March, and the total yearly energy 
delivered to the grid in the year 2018 was 35892.22 MWh, which is close to the 
expected energy (36,414 MWh). 


e The monthly average reference and final yields of the system are in the range 


5.92-8.1 (h/day) and 4.39-5.56 (h/day), respectively, and the annual monthly 
average daily final yield is 4.98 (h/day). The monthly daily average yields 
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(final and reference) increase from May to August due to higher irradiation and 
extended sunshine duration. 


The minimum and maximum efficiency and performance ratio were 8.34 and 
55.29% in July, and 12.37 and 81.96% in January, respectively. The minimum 
capacity factor was 18.31 in July and the maximum was 23.17% in March. The 
minimum total energy losses was 1.07 h in January and the maximum was 

3.55 h in July. Furthermore, the annual average daily PV system efficiency, per- 
formance ratio, capacity factor, and total losses were 10.82%, 71.71%, 20.76%, 
and 2.04 h/day, respectively. 


The reasons for the losses are generally due to dust, ambient temperature and 
solar radiation. 


The results obtained show that the installation of photovoltaic stations in Adrar in 
Southern Algeria gives good results encouraging investments in this field and region. 
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Appendices and nomenclature 


Pimp maximum power (W) 

Vimp maximum power voltage (V) 

Lap maximum power current (A) 

Voc open circuit voltage (V) 

Isc short circuit current (A) 

Opmp temperature coefficient of Pmp (%/°C) 
Osc temperature coefficient of Isc (%/°C) 
Oyoc temperature coefficient of Voc (%/°C) 
STC standard test condition (25°C and 1 kW/m’) 
Ta ambient temperature (°C) 

W, wind speed (m/s) 

Pp AC power (kW) 

Eac AC energy generated (kWh) 

Yr final yield (h) 

YR reference yield (h) 

TL total energy losses (h) 

CF capacity factor (%) 

PR performance ratio (%) 

H, total in-plane solar irradiation (kW/m?) 
Hote reference irradiation (=1 kW/m’) 

Nsys system efficiency (%) 

N number of days per month 


d, m, and y daily, monthly, and yearly period 
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Chapter 6 


Modeling of the Small Wind 
Energy in Saharan Region 
of Algeria 


Harrouz Abdelkader, Fadila Tahiri, Boussaid Brahim 
and Fatiha Bekraoui 


Abstract 


In the last century, research for the wind power industry started to gain impor- 
tance in the field of renewable energies. This research is mainly on the side of big 
wind power. This wind form, are grouped and connected to the electricity grid. In 
the other form, the small wind power for production in isolated areas. This wind are 
applicated for example: telecommunication stations, sailboats and isolated areas. 
This power is from 100 W to a few tens of kilowatts. They are applied to power 
installations, for example, telecommunication stations, sailboats, and isolated areas. 
The use of small wind power in the region of Adrar—the southwest of Algeria—is 
an economical and durable solution. This chapter will focus on the modeling of the 
main components of a small wind turbine adapted to Saharan regions. The wind 
chain consists of a wing coupled directly to a synchronous generator that delivers on 
a continuous bus via a rectifier; it is the overall structure of the chain that we retain 
for modeling in this chapter. In order to control the system, the modeling of this 
study touches all the parts of the system: the turbine, the generator PMSG, and the 
converter with the load. At the end, simulation results are presented to show the 
good performance of the choice of control type that is applied to the wind system. 


Keywords: wind, modeling, electricity, MPPT 


1. Introduction 


Wind energy is one of the fastest-growing renewable energy in the world. 

The wind resource comes from the movement of air masses that is directly related 
to the sunshine of the earth. By warming certain areas of the planet and cooling 
others a pressure difference is created and the air masses are in perpetual motion. 

Since the use of the windmill, wind sensor technology has been evolving. It was 
in the early 1940s that real prototype wind turbine blades were successfully used to 
generate electricity. Several technologies are used to capture wind energy (vertical 
axis or horizontal axis sensor) and the sensor components, the mechanical charac- 
teristics of the wind turbine, and the synchronous and asynchronous machines are 
efficient [1]. 

In the region of Adrar (located in the southwest of Algeria), there are isolated 
sites, where the use of this abundant source of energy, is an economical and durable 
solution. The permanent magnet synchronous machine is characterized by high 
volumetric torque, very low inertia and low inductances. All of these features 
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Wind Turbin PMSG AC/DC converter DC bus 


Figure 1. 
Schematic diagram of a wind chain. 


provide the generator with high performance, high efficiency and better control. 
This puts this machine in competition with the asynchronous machine [2]. 

To meet the demand for autonomous electric charge at remote sites, small wind 
turbines smaller than 5 kW have become a very important system for the produc- 
tion of electricity. Current research will study the functional analysis, control and 
modeling techniques appropriate for this small turbine system [3]. 

Modeling of the wind turbine recently published by many authors, Harrouz et al. 
[4] have modeling of small wind energy based on PMSG in the south of Algeria. 
Langlois [5] have modeling and study the performances of wind system out the link. 
Tran [6] has present work of optimal design integrate of wind energy “passive” 
robustness analysis, experimental validation. Jedli and Hidouri [7] present a power 
drive scheme for an Isolated pitched wind turbine water pumping system based on 
DC machine, this work can be a good application of win system in Sahara region. 
Modeling, simulation, and control of a variable speed wind turbine based on the 
asynchronous dual-feed generator [8]. 

The wind energy conversion chain is schematized as shown in Figure 1. 

This consists of a wing coupled directly to a synchronous generator that delivers 
on a DC bus via a diode rectifier; it is the structure that we retain for this work. 

The control strategies of this system and their possible load connection inter- 
faces must make it possible to capture as much energy as possible over the widest 
possible range of wind speeds, with the aim of improving the profitability of wind 
turbine installations. 

The section after will focus on the modeling of the main components of a small 
wind turbine adapted to saharan regions such as Adrar. 


2. Modeling of the small wind 
2.1 Modeling the wind 


The wind is highly variable, both geographically and over time. It varies from 
one place to another, from 1 day to another, and from 1s to the next. It is the energy 
input vector of a wind chain, determining for the calculation of electricity produc- 
tion, therefore for profitability [2]. 

The dynamic properties of the wind are crucial for the study of the whole system 
of energy conversion because the wind power, in the optimal conditions, is at the 
cube of the wind speed (see formula (2)) These fluctuations are influenced by the 
displacement of air masses at altitude, but also by the relief, the type of soil cover 
and the thermal stability of the atmosphere. The wind varies in direction and 
intensity, but for the purposes of the model, it is limited to the variation of the wind 
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in intensity, in a single direction in order to determine a valid and representative 
sequence of real winds [2-8]. 

Different approaches used in the literature for the generation of a synthetic 
series of wind, in our case, the wind speed will be modeled by a sum of several 
harmonics [2, 4, 9] 


Voent(t) =A gs sin ( w,.t ) (1) 
k=1 


where az: amplitude of the harmonics; œp: frequency of harmonics. 
The wind speed is represented by the function: 


V(t) = 10 + 0.2sin (0.1047t) + 2 sin (0.2665t) + sin (1.2930t) + 0.2 sin (3.66451) (2) 


It should be noted that this particular wind profile corresponds to [9-12]. 


2.2 Modeling of turbine 


A wind turbine is a machine that, by definition, transforms wind energy into 
mechanical energy. To begin, it is necessary to quantify the energy source available, 
i.e., the energy associated with the wind. If the wind has a certain speed “V” at a 
given moment and crosses a certain area “A,” the instantaneous wind power is 
given by the following relation 


Pg= ; „pA V? (3) 
where p is the density of the air, which is approximately 1.2 kg/m?; A is the area 
swept by the turbine in mĉ; and V is the wind speed in m/s. C, (power coefficient). 
The turbine used in our work is a “Savonius” wind turbine with a vertical axis 
(see Figure 2). 
Surface “A” given by the following formula: 


A=2R.H (4) 


where R represents the radius of the wing (m) and H its height (m). 
The output power is given by the following equation that we will normalize in pu: 


Pm = Cp DpHRV" (5) 


where C, is a coefficient that expresses the efficiency of the wing in the transfor- 
mation of the kinetic energy of the wind into mechanical energy, which is in fact often 


S$=2R.H 


Figure 2. 
Wind vertical axis of Savonius type [7]. 
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given as a function of the reduced speed (the specific speed A), which is the ratio of the 
linear speed at the end of turbine blades reduced to wind speed, defined by: 


1 = (R.Q)/V (6) 


where Q the angular rotation speed of the blades in rad/s. 

The evaluation of the power coefficient is a data specific to each wind turbine. 
From the readings taken on a wind turbine, the expression of the power coefficient 
has been approximated, for this turbine [13, 14], by the analytical equation as a 
function of A (resulting from the interpolation) according to: 


C,(A) = —0.12992 43 — 0.116810? + 0.45406 A (7) 


Figure 3 represents the characteristics of power coefficient as a function of A, it 
presents a maximum power factor Cymax = 0.19 for Amax = 0.55. 

The wind torque noted C, is obtained from Eq. (5) and using the expression (6) 
of the reduced speed (Figure 4): 


Pm _ Gp(d).p.R? HV? 
Q A 


C: = (8) 

The Model of the wind torque under Simulink presented in Figure 4. The 
Figure 5 shows the characteristic of the turbine through the wind power extraction 
of the turbine as a function of the wind speed. 


o 
za 
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o 
o 9 
© & 


Power Coefficient 


į Ae i 
0 0.2 0.4 0.6 0.8 1 
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Figure 3. 
Power coefficient C, as a function of A for Savonius wind turbine. 


= 
Torque coefficient + 


Or 
Figure 4. 


Model of the wind under Simulink. 
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Nominal power - - ~ 


Optimal speed 


Mechanical power (MW) 


10 15 20 25 30 
rotor speed (rpm) 


Figure 5. 
The output power for different values of wind speed (m/s). 


It can be seen from Figure 5 that the electric power varies with the variation of 
the mechanical angle of the turbine. The shape of the electric power as a function of 
the mechanical angle is seems to be that of the Cp as a function of à. We also note 
that the higher the wind speed, the greater the power delivered because the power is 
expressed by the cubic wind speed. 


2.3 Mechanical modeling of the turbine-generator coupling 


The dynamic (mechanical) behavior of the turbine and generator assembly can 
be represented by the following relation: 


J = Ce Cem Cr ( 9) 


where Q the speed on the generator shaft, Cy the friction torque, Cem the elec- 
tromagnetic torque developed by the generator, C, the mechanical torque applied to 
the alternator shaft and J is the total moment of inertia, calculated with: 


J = Jen -4 (10) 
We take note that: 
Cr =fp Q (11) 


where f, is the coefficient of viscous friction in N.m. 

The wind torque provided by the turbine, drives the generator. Taking into 
account the torque drops caused by the friction and the inertia of the turbine, we 
can establish the electromechanical model of Figure 6. 


2.4 Model of the GSAP 


According to Chapter, the Permanent Magnet Synchronous Generator (PMSG) 
is classically modeled in the Park coordinate system, giving rise to the following 
equation: 
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dI 
Va = —Ryla — Lage + Lgoy 
dl, | 
Va = Rs Tg Ly a Lgq@lq 4 pro (2) 
dQ 
= Cm — Com — F.Q 
J dt 


3 
Com = 5P| (Lg — Ld lay + 9 1, 


where: 0 is the angle between a reference axis of the stator and an axis of the 
north pole of the rotor, p-the number of pairs of poles, R, the resistance of a stator 
phase. 

Also, @, is the excitation flux produced by the permanent magnets, L4 and L, are 
the equivalent inductances on the d axes, respectively q. 

Cem is the electromagnetic torque, C,, is the motor torque applied to the gener- 
ator with F the coefficient of friction, J: moment of inertia and p the number of pairs 
of poles. 

In order to know the dynamic characteristics of the various state variables of the 
machine and to validate the mathematical model, we proceed to a simulation of the 
machine using the Runge-Kutta algorithm of order 4. 


cp).p.R? av? 


À 


Figure 6. 
Model of wind turbine torque. 


Figure 7. 
Stator current lap: of PMSG after simulation. 
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Uabe (V) 


Figure 8. 
Stator voltage of PMSG after simulation. 


The response to an empty speed step is very fast. The torque peaks at startup and 
is proportional to the stator current. At full speed, the speed decreases while the 
torque increases to drive the load (Figures 7 and 8). 

To see the efficiency of our model, the next part will be the subject of test on 
matlab/simulink with real parameters of the wind turbine and the wind speed. 


3. Simulation of model PMSG with wind system 


The wind will be modeled in deterministic form by a sum of several harmonics: 


V(t) = 10 + 0.2sin (0.1047t) + 2sin (0.2665t) + sin (1.2930t) + 0.2 sin (3.66451) 
(13) 


The simulation results of the permanent magnet synchronous generator associ- 
ated with the wind turbine and with the wind speed simulated by Eq. (13) are given 
by Figure 9: 

We notice that the results take the form of the wind, the coefficient of power is 
0.13 (it is the maximum value to have a maximum of power extracted from the 
wind). The active power is of the order of 400 W. 

But this system presents fluctuations due to the variation of the speed of the 
wind for that the system needs a command to maintain this power at its nominal 
value. The next partie present the proposed control of system. 


4. Command of wind system 


The wind system studied is a complex system composed of several subsystems 
belonging to several different physical domains, namely, traditional mechanics, 
electrical engineering and electrochemistry. Energy from the wind passes through 
the wind turbine that is coupled to the GSAP generator and allows the transforma- 
tion of mechanical energy into electrical energy. The electrical components such as 
the static converter disposed downstream of the generator, have a role of active 
adaptation of the characteristics of the electrical energy between the generator and 
the final load. 
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Wind speed (m/s) 


i 
| 
| 


Figure 9. 
The results of the wind turbine simulation. (a) The wind speed as a function of time. (b) The speed of the 
turbine as a function of time. (c) Power factor. (d) The electrical speed of the GSAP as a function of time. 


The PWM rectifier was studied by simulation under MATLAB/SIMULINK 
according to the diagram of Figure 10, the current was controlled in the reference 
frame abc by hysteresis regulators, the reference of the voltage at the output of the 
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MSC 


Wind profile 


Figure 10. 
Structure of control the wind system. 


wind speed (m/s) 


Figure 11. 
Wind speed (m/s). 


Voltage (V) 
25 
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Oo 
Time (s) 
Figure 12. 
Voltage rectified. 


rectifier is taken equal to 45.8 V, We simulated the model of wind profile is given by 
Figure 11. Figure 12 shows the DC voltage at the output of the rectifier with its 
proposed reference, it is observed that the DC bus voltage remains on average equal 
to 45.8 V. Figure 13 shows the zoom of the rectified voltage. 

According to the simulation results obtained, we conclude that the response of 
the voltage at the rectifier output (rectified) to a speed variation is relatively fast 
and does not exceed 2% of the reference value, during disturbances. 
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Figure 13. 
Zoom of voltage rectified. 


In spite of the increase or the decrease of the wind speed, the shape of the 
voltage of the continuous bus is established with 45.8 V with a time of response 
which depends on the control of the rectifier, of the order of 0.5 s in the example 
treated. 


5. Conclusion 


This chapter has established a global model of the wind energy conversion chain 
which consists of a GSAP permanent magnet synchronous machine associated with 
a wind turbine. The turbine used is the vertical axis Savonius type. They presented 
the simulation results of the complete wind energy production chain. The end of 
this chapter presents the application of the control of this wind energy production 
system. This control shows that the response of the voltage at the rectifier output 
(rectified) to a speed variation is relatively fast. In spite of the increase or decrease 
in wind speed, the speed of the continuous voltage is established with a response 
time which depends on the control of the rectifier. 
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Implement Using KY Converter 
for Hybrid Renewable Energy 
Applications: Design, Analysis, 
and Implementation 


Pushpavalli Murugan and Jothi Swaroopan Nesa Mony 


Abstract 


This chapter mainly focuses on meeting the energy demand and methodolo- 
gies of renewable energy. Nowadays, researchers are mainly focusing on renew- 
able energy from the sun, wind, biomass, etc. due to energy crises and the lack 
of non-renewable energy. The potential for solar energy is high and this demand 
can best be met with hybrid systems, which can provide an uninterruptible power 
supply. This chapter looks at the performance metrics of hybrid energy as well as 
the methodologies and various control techniques connected with power manage- 
ment. The chapter also defines the photovoltaic (PV)-based, novel, dual KY boost 
converter. Dual PV sources act as input for the dual KY boost converter to generate 
as much energy as possible from the dual PV system, using the inverter module to 
produce single-phase alternating current output. A dual KY boost converter can 
provide higher maximum power, a faster response, and smaller voltage ripple. KY 
boost converters are designed to generate stable output values according to various 
conditions because of various control techniques and the maximum power point 
tracking control algorithm. 


Keywords: hybrid energy, grid system, control techniques, dual KY converter 


1. Introduction 


In India 53% of energy is produced from thermal power stations. The major prob- 
lems with thermal energy are lack of coal, ash production, fuel emission, high running 
costs, and disruption to the ecology. This chapter mainly focuses on renewable energy 
sources like photovoltaic (PV), wind, and stable conversion using a hybrid system. In 
this introduction, the focus is on wind energy-related challenges because wind energy 
is not a linear energy. Wind is always fluctuating in nature so an appropriate controller 
is always required for optimal power. Here, an appropriate controller for an induction 
generator to achieve steady-state output will be described. An induction generator acts 
as a drive and a neural network controller is introduced with an ant colony optimiza- 
tion method. The controller is intended to drive a turbine to extricate extreme power 
from the wind. This strategy is dependable and reliable [1]. 

A wind energy system requires constant power at the conversion stage, 
and to avoid this variable, a frequency transformer technique is adopted. This 
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methodology proposed in this chapter is a standalone wind system. The suggested 
strategy is modest, and has basic control with no harmonics [2]. The system consists 
of solar power, wind power, a diesel engine, and an intelligent power controller. 

To attain an active and balanced response for active power control, the controller 
consists of a radial basis function network and an improved Elman neural network 
for maximum power point tracking (MPPT). The pitch angle of the wind turbine 
and the PV system uses a radial basis function network, where the output signal 

is used to control the DC/Elman neural network DC boost converters to achieve 
MPPT [3]. Different types of permanent synchronous generator with various power 
capacity wind farms are connected to the common grid. This setup has three col- 
lected models of variable wind speed fitted with permanent magnet synchronous 
generator (PMSG) wind turbines for dynamic investigation [4]. 

This chapter mainly targets a suitable converter for a renewable energy system. 
It discusses the KY converter, which is suitable for hybridizing the energy. The 
KY converter gets its name from a paper titled “KY converter and its derivatives” 
by K.I. Hwu and YT. Yau (KY converter). A dual KY boost converter is presented, 
which is a KY converter combined with a normal synchronously rectified boost 
converter. The input and output inductor currents are continuous, with a higher 
voltage conversion ratio suitable for low-ripple applications. By using soft switching 
with the surge current suppressed the device can be used in high-power applica- 
tions. The KY converter has better performance than a boost converter [5] and 
provides a fast response to reach the maximum power point compared to a boost 
converter. A novel voltage-bucking/boosting converter is called a KY buck-boost 
converter. It has fast transient responses, synchronous rectification, load regula- 
tions, and low-output voltage ripples. Cuk, Sepic, Zeta, and Luo converters could 
also be considered but they possess right half-plane zeroes, thereby causing system 
instability and slow load transient responses [6]. The operation of KY converters is 
also given, along with experimental results to verify ideal line and load regulations, 
and low output voltage ripples due to non-pulsating output currents. 

Most solar panel energy efficiency ratings range from 15 to 20% [7]; therefore, 
various MPPT algorithms have been introduced. Numerous controllers are utilized 
to find MPPT and the efficiency is increased [8]. Power from any renewable energy 
source can be delivered to the grid individually or a balanced MPPT feature can be 
realized for renewable energy, power converters, and control algorithms along with 
purposeful energy resources for efficient operation of the microgrid. The proposed 
microgrid, based on hybrid energy resources, operates in autonomous mode and 
has an open architecture platform for testing multiple different control configura- 
tions. AC/DC, DC/AC, and DC/DC converters are integrated with the distributed 
energy system due to different types of output voltages [9]. Several algorithms that 
are used to operate DC/DC converters around the MPPT are reported in the litera- 
ture. Comparing all algorithms, a fuzzy logic controller coupled to other controllers 
works well under partial shading conditions. There is also a newly designed fuzzy 
logic controller coupled to a Hopfield neural network maximum tracking technique 
[10]. The MPPT technique is very important because it increases the energy effi- 
ciency of a renewable energy system. 

Another major challenge is hybrid renewable energy. The configuration of the 
hybrid renewable energy system and interfacing power converters for connecting 
the energy sources to the AC bus is extensively discussed. An outline of the control 
process in a hybrid renewable energy system and the application of the relevant 
control methods for system stabilization, effectively inducing real power and 
proper load-sharing methods, are available. Different approaches for hybrid renew- 
able energy system design and control methods for power converters in the recent 
research literature are also briefly discussed. The AC bus-linked hybrid renewable 
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energy system configuration reduces the number of power conversion stages and 
losses in power transferred to the load/utility. The master/slave control with the 
droop concept does not require a communication link and provides good load 
sharing, such as the flexibility, expandability, and modularity of the hybrid renew- 
able energy system. Both single-master and multimaster approaches are used in the 
inverter control strategies [11]. In addition, a multi accumulated model utilizing 
comparable wind speed display has excellent similarity to the total model in transient 
and steady states and it also expends less time in reproduction when contrasted with 
the entire model [12]. Finally that the hybrid control system challenges not so distant 
future with its very own potential which clears another time with powerful and 
proficient energy balance. Globally, renewable energy has a green future and will 
undoubtedly become a credible alternative to fossil fuels [13, 14]. 

In conclusion, hybrid energy devices combining a variable velocity wind turbine 
and PV array-generating device could deliver a non-stop power to weight ratio to 
produce a viable hybrid controller. The hybrid controller manages the strength flow 
between the system’s additives such that the cost is minimized and load necessities 
might be met throughout the year [15]. Hybrid renewable energy’s main challenges 
are energy storage capacity and energy management. To utilize 100% renewable 
energy, the design requires intelligent charge control, battery state of charge, and 
estimation of the impact of various working routines on battery life. Different 
architectures of energy sources and power management on peak load sharing are 
analyzed with a reduction of varying losses [16]. 

The next important analysis concerns a multi-input system. A novel multi-input 
inverter is used for the grid-connected hybrid solar/wind power system to simplify 
the power system and reduce the cost. The multi-input inverter is designed with 
a buck/buck—boost fused to a multi-input DC/DC converter and a single-phase 
bridge-type DC/AC inverter [17]. A semi-isolated multi-input converter for a 
hybrid PV/wind power charger system that can simplify the power system, reduce 
cost, deliver continuous power, and overcome high-voltage-transfer ratio problems 
is proposed [18]. 

Various literatures on the use of hybrids with renewable energy systems for 
uninterrupted power supply to the load have been analyzed. Another requirement is 
the MPPT algorithm along with the multi-input converter, which is used to reduce 
the volume of the system. Hybrid renewable energy systems in combination with 
existing setups are presented to investigate the possibility of specialized aggressive- 
ness [15]. In this chapter, the operational principle of the proposed semi-isolated 
multi-input converter is explained: a novel, executed multi-input KY boost 
converter for hybrid systems associated with the grid. The control parameters are 
reliant on the wind, sun, storage, and grid conditions and checked under dynamic 
conditions. DC bus voltages are used to deal with the contribution of various work- 
ing mode controls amid different working conditions [19, 20]. 


2. Hybrid energy system methods 
2.1 Types of integration system 


There are numerous approaches to coordinate distinctive alternative energy 
to form the integration of hybrid systems. For the most part, strategies can be 
arranged into three classes: DC coupled, AC integrated, and hybrid integrated. The 
AC-integrated system can be grouped into power-frequency AC-integrated and 
high-frequency AC-integrated systems. The strategies of a DC-coordinated setup 
are shown in Figure 1. 
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Direct Current 


(—) Loads 


Figure 1. 
Block diagram of a DC-coupled integration system. 


2.1.1 DC-coupled integration system 


The DC sources may be associated with the DC bus directly if they are suitable. 
If there are any DC loads, they can be clearly related to the DC bus, or through the 
DC/DC converters, to fit DC voltages to DC loads [21]. 


2.1.2 AC-coupled integration system 


AC coupling can be split into two further categories: power frequency 
AC-coordinated and high-frequency AC-coordinated systems. An AC-coordinated 
system is shown in Figure 2, where diverse energy sources are coordinated through 
their control switching interfacing circuits to a power recurrence AC bus. Inductors 
may likewise be required between the control switching circuits and the AC bus to 
accomplish the desired control through the board. 


2.1.3 Hybrid-coupled integration system 


Rather than interfacing all the distributed generation sources to only one DC or 
AC bus, as already discussed, the distinctive distributed generation sources can be 
associated with the DC or AC bus of the mixture system. Figure 3 demonstrates a 
hybrid-coupled coordinated system, where distributed generation assets are associ- 
ated with the DC bus as well as the AC bus. In this arrangement, energy sources can 
be coordinated specifically without any additional circuits. Therefore, the frame- 
work can have higher energy effectiveness and diminished cost. Meanwhile, control 
and power flow management may become increasingly difficult compared to DC 
and AC bus-coordinated systems [22]. 


2.2 System arrangement 


The system comprises renewable energy and alternative energy. The renew- 
able energies are wind, solar, and biomass. Alternative energies are a diesel 
generator and a storage system. In a storage system, a battery, fuel cells, superca- 
pacitors, and aqua electrolysis are used [23]. The level of dependability of a PV 
and wind-coordinated framework to suit a specific load can be given by a mix of 
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Figure 2. 


AC Bus 


Figure 3. 
Block diagram of a hybrid-coupled integration system. 


appropriately estimated wind energy, PV, battery units, and alternative sources. 
In atmospheres, alternative sources are expected to be highly dependable also, 
keep away from gross over-plan of the sunlight based and wind framework. 

In other words, an assistant energy provider can be viewed as an optional 
technique for proposing an alternative energy source fused into the framework 
compared to expanding equipment sizes unnecessarily for those periods when 
energy generation is reduced [24]. The PMSG-based wind turbine works at the 
greatest power point following the MPPT mode to extract power from the wind, 
although wind power may need to be reduced under specific conditions. Wind 
reduction can be accomplished by utilizing the turbine’s capacity and pitch 
control systems [25]. 

The execution of crystalline silicon PV modules limits the physical elements 
of the PV cell material, the temperature of PV array cells and the sun situated 
irradiance revealed on the sun’s situated cells. One solved proper model for 
the most outrageous power yield of PV modules is used [4-7, 12, 19, 20, 22-24, 
26-30]. 
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Changes in wind speed and sunlight-based illumination are associated with the 
wind turbine generator, and PV can recreate a variety of intensities of AC and DC 
sources and test the MPPT control calculation [28]. 

Regarding the size of the battery bank to use for these sorts of uses requires 
examination of the battery’ charge and discharge capabilities, including burden, 
yield, and in the case of sunlight-based or elective critical sources, the working 
temperature and capability of the charger and other parts of the structure. By and 
large, problems occur while charging the battery bank and its capability drops when 
the battery ages or is misused. At the point when the aggregate yield of solar and 
wind energy is more than the energy requested, the storage unit is charged. Two 
properties of the battery are associated with the hybrid system’s execution, i.e. the 
state of charge and the charge voltage [29]. 

The hybrid system grid plays a vital role and can work in two modes. Grid Tied 
mode the power converters to produce constant voltage to a DC bus and exchange 
of power can be accomplished using DC and AC buses, while the grid working 
as an autonomous mode storage unit plays a vital role in balancing voltage and 
power. 


2.3 Expert control 
2.3.1 Intelligent controller 


The control framework for hybrid energy system designs should limit fuel 
utilization by expanding energy from inexhaustible sources. Nonetheless, there are 
control changes by which the sustainable power source can change, which causes 
unsettling influences that can influence the nature of the energy conveyed to the 
heap. To deal with the stream of energy effectively, with fewer harmonics, a con- 
trolling strategy should be developed. 

The energy sources are modeled from PV, wind, and battery. To coordinate all 
the converters used in the hybrid system, intelligent controllers are essential. A 
proportional-integral controller, fuzzy logic controller, neural network, and genetic 
algorithms are used as an intelligent controller [7]. 


2.3.2 Control methods 
e The hybrid system has three types of controller. The controller used in PV and 
wind is said to be a generator side converter to obtain maximum power point 


tracking from solar and wind. 


e A grid side converter is used to sustain constant voltage in a DC bus and fre- 
quency in an AC bus, to control the exchange of reactive power in the AC bus. 


e A storage side converter maintains constant DC voltage and state of charge of 
the battery. 


These converters will coordinate the entire hybrid system. 
2.3.3 Power management 

Power management is categorized as two types: rule based and optimization 
based. Rule based mainly depends on real-time applications. Depending on the 
grid requirements, power flow will be controlled from source to grid. Optimization 


based depends on the cost of energy. Its main focus is to minimize cost. 
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As the arrangement of hybrid renewable energy and alternative energy sys- 
tems as a microgrid expands, so does the necessity for continuous energy for such 
systems. Powerful correspondence between the separate energies of the microgrid 
is imperative and merits further consideration. In addition, efficient methodologies 
and institutionalization are required for the production and safe organization of 
such systems. 


2.3.4 Grid-connected mode 


At the point where the hybrid grid works in islanding mode, the boost converter 
and the consecutive AC/DC/AC converter of the doubly-fed induction generator 
may work in on-MPPT or off-MPPT mode depending on framework control equality 
and criticality goals. The essential converter is a voltage source that gives a consistent 
voltage that repeats for the AC framework and works either in inverter or converter 
mode for smooth power exchange between the AC and DC link. The battery con- 
verter works either in charging or discharging mode and is reliant on the power 
balance in the framework. DC interface voltage is maintained by either the battery 
or the boost converter and is reliant on structure working conditions. Forces under 
different load and supply conditions should be adjusted. 

When the hybrid system works in grid-connected mode, the goal of the boost 
converter is to pursue the MPPT of the PV display by coordinating its terminal volt- 
age. The sequential AC/DC/AC converter of the doubly-fed induction generator is 
controlled to guide the rotor side current to achieve MPPT and to synchronize with 
the AC framework. The energy overflow of the hybrid grid framework can be sent 
to the utility framework. The activity of the battery for energy stockpiling is less 
imperative since the power is adjusted by the utility framework. For this situation, 
the fundamental purpose of the battery is to eliminate customary control trade 
between the DC and AC link. The DC/DC converter of the battery can be controlled 
as an energy pad using this strategy. The principal of converter is intended to work 
bidirectionally to fuse corresponding normal for wind whats more, sun-powered 
sources. The goals of the primary converter are to maintain a steady DC-connected 
voltage for the variable DC stack and to synchronize with the AC connection and 
utility framework. 


3. Design of a dual KY boost converter 
3.1 Working of a dual KY boost converter 


Figure 4 demonstrates the proposed dual KY boost converter built using the 
KY converter joined with the boost converter. The dual KY converter is made 
using four switches M1, M2, M3, and M4, one diode D, one output inductor Lo, 
two input inductors Lil, Li2, a support capacitor Cs, and one output capacitor 
Co. The depository capacitor Cs is a support between the KY converter and the 
conventional synchronously rectified boost converter. The multi-input DC voltage 
is 12 V and 110 V are obtained at the yield. Because of the conduction status of the 
switches the dual KY converter has two modes of operation. Figure 5 shows mode 
1 operation of the dual KY boost converter. M1 and M4 are turned off and M2 and 
M3 are turned on. The diode D is forward biased and begins to conduct, and Cs 
is released. In this manner, the voltage across Lil is V1, making Lil charged, while 
the voltage across Lo is Vo minus VCs, subsequently making Lo demagnetized. 
Additionally, the current flowing through Co is equivalent to ILo minus the current 
moving through RL. 
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Figure 4. 
Dual KY boost converter connected to DC input. 
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Figure 5. 
Mode 1 operation of the dual KY boost converter. 


Figure 6 shows mode 2 operation of a dual KY boost converter. M2 and M3 are 
turned on and M1 and M4 are turned off. Consequently, the voltage across Lil is VCs 
subtracted from voltage source V1, subsequently making Lil demagnetized, while the 
voltage across Lo is Vo subtracted from 2VCs making Lo polarized. Additionally, the 
current flowing through Co is equivalent to ILo minus the current moving through 
RL, while the current moving through Cs is equivalent to the total of ILi1 and —ILo. 
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The input voltage V1 = V2 = 12 V and the input power of 12 W is made up from each 
source of 6 W. Figure 8 shows the output power of a dual KY boost converter of 22 W 
with output voltage of 110 V. The results are shown in Figures 7 and 8. 


Figure 6. 
Mode 2 operation of a dual KY boost converter. 
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Figure 7. 
Dual KY boost converter output voltage is 110 V. 
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Figure 8. 
Total input power = 12 W and dual KY converter output power = 22 W. 
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3.2 Disadvantages of the existing method 


Regular converters for interfacing numerous sustainable sources uti- 
lize a typical DC connection or normal AC recurrence interface, which 
has disadvantages like the necessity for increased bundling, and expanded 
switches over the ordinary structure becoming delineated [31]. The multiport 
converter has numerous points of interest compared to the traditional struc- 
ture regarding the number of energy gadgets and transformation steps that 
can enhance system cost. Certain non-segregated voltage-bucking/boosting 
converters make the system cost effective. The interleaved converter 
can likewise be utilized. It comprises two single-stage support converters in 
parallel. 


3.3 Advantages of the proposed method 


This dual KY support converter topology highlights higher proficiency with 
less powerful hardware gadgets and fewer power change processes. The proposed 
dual KY support converter has the benefits of basic topology and the least number 
of energy switches. The dual KY support converter has persistent data and yields 
inductor streams and high-voltage proportions. The dual KY support converter 
structure is promising from the perspectives of minimal effort, concentrated 
control, and small size [27]. 


4. Dual KY boost converter connected to inverter 


The dual KY converter is connected to a single-phase inverter. Analysis is done 
for a PV array input of 18 V given to V1 and V2 and the KY converter voltage is 230 
V with a current of 3.5 A. The dual KY boost converter output is connected toa 
single-phase inverter. The single-phase inverter output voltage is 230 V and the total 
harmonic distortion (THD) value is nearly 8.22%. 


4.1 Dual KY boost converter connected to PV array input 


Table 1 listed the design value of PV array. 


Parameters Values 
Isc 5.96 A 
Voc 6.2V 
Series resistance 0.18 ohm 
Parallel resistance 2 ohm 
No. of modules connected in series 4 
No. of modules connected in parallel 4 
Irradiance 1000 W/m? 
Temperature 25°C 
Table 1. 
Design value of PV array. 
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4.2 Simulation results of dual KY boost converter connected to inverter 


The above design values of PV arrays are connected as input and the waveforms are 
obtained. The analysis is done for a PV array input of 18 V shown in Figures 9 and 10 
and an output of 230 V obtained across the load resistance shown in Figure 11. The 
component determinations are given as follows: (1) evaluated PV array input voltage 
is 18 V; (2) DC voltage is 110 V; (3) exchanging recurrence is 2 kHz; (4) estimation of 
Lil and Li2 is 400 uH and the estimation of Lo is the same as that of Li; (5) one 1 pF 
capacitor is decided for Cs; (6) one 470 pF capacitor is decided for Co; and (7) param- 
eters of the PI controller, kp and ki, are set to 1 and 10 individually. These outputs are 
replicated using the simulation shown in Figures 12 and 13. 

Figure 14 shows the perturbation and observation algorithm used to gen- 
erate the duty ratio for switches. The simulation also shows a dual KY boost 
converter connected to battery terminals and the state of charge is checked. The 
dual KY boost converter output is connected to a single-phase inverter. The 
single-phase inverter output voltage is 230 V and the THD value is nearly 8.22%. 


Figure 9. 
Dual KY boost converter connected to PV array input. 


Figure 10. 
Dual KY boost converter connected to single-phase inverter output. 
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The corresponding waveforms are shown in Figures 15-18. Table 2 shows the 
output values of KY boost converter. Table 3 shows the comparison of KY boost 
and boost converter. 
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Figure 11. 
Dual KY boost converter DC output voltage of 230 V. 
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Figure 12. 
Dual KY boost converter DC output current of 3.5 A. 
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Figure 13. 
Ppv = 450 W 
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Figure 14. 


Perturbation and observation algorithm for generating pulses. 


Parameters Values 
Isc 26.3 A 
Voc 19.5 V 
PV panel voltage 18V 
PV panel current 25A 
PV panel power 450 W 
KY boost converter voltage 230 V 
KY boost converter current 3.5 A 
Inverter AC output voltage 230 V 
Inverter AC output current 5A 
Power factor 0.93 
Table 2. 
Output of KY boost converter. 
Converter Vin Vo THD (%) 
KY boost converter 18V 230 V 6.58 
Boost converter 18 V 110 V 10 
Table 3. 
Comparison of output voltage and THD. 
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Figure 15. 
Inverter AC output voltage is 230 V. 
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Figure 16. 
THD = 8.22% for single-phase inverter output voltage. 
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Figure 17. 
Single-phase inverter AC output current is 5 A. 
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Figure 18. 
THD = 6.58% for single-phase inverter output current. 
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5. Conclusion 


Because of fast improvements in the field of sustainable power sources and 
expanding costs of ordinary items like oil and gas, the use of sustainable power 
sources has turned out to be exceptional in hybrid systems. Solar PV and wind 
systems have been promoted around the globe on a comparatively large scale. 

This chapter proposed a hybrid energy system consisting of a wind turbine and 

a PV source to supply continuous power to the load. A controller is used to track 
the maximum power from renewable sources. Simple and economic control 

with a KY converter, which is a DC/DC converter, is used for maximum power 
extraction from the wind turbine and PV array. Power from the dual PV array can 
be transported to the utility grid. At the same time, the MPPT attribute can be 
acknowledged for the dual PV system. The substantial range of PV array voltage, 
current, and temperature variations was caused by various insolations. This solitary 
dual generation system can adequately separate extremes of power from the two PV 
sun-oriented energy sources. Boosting the DC voltage (18-230 V) to an adequate 
level utilizing the dual KY boost converter and obtaining pure AC voltage (230 V) 
from the inverter are key to understanding the above targets. A model of a hybrid 
dual solar system arrangement utilized with a control strategy has therefore been 
created. The KY converter will offer low THD (6.58%), high power factor under 
nonlinear stacking conditions, and pronounced unique reactions under transient 
stacking conditions. 
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Chapter 8 


Survey on Photo-Voltaic Powered 
Interleaved Converter System 


M.L. Bharathi 


Abstract 


Renewable energy is the best solution to meet the growing demand for energy 
in the country. The solar energy is considered as the most promising energy by the 
researchers due to its abundant availability, eco-friendly nature, long lasting nature, 
wide range of application and above all it is a maintenance free system. The energy 
absorbed by the earth can satisfy 15000 times of today’s total energy demand 
and its hundred times more than that our conventional energy like coal and other 
fossil fuels. Though, there are overwhelming advantages in solar energy, It has few 
drawbacks as well such as its low conversion ratio, inconsistent supply of energy 
due to variation in the sun light, less efficiency due to ripples in the converter, time 
dependent and, above all, high capitation cost. These aforementioned flaws have 
been addressed by the researchers in order to extract maximum energy and attain 
hundred percentage benefits of this heavenly resource. So, this chapter presents a 
comprehensive investigation based on photo voltaic (PV) system requirements with 
the following constraints such as system efficiency, system gain, dynamic response, 
switching losses are investigated. The overview exhibits and identifies the require- 
ments of a best PV power generation system. 


Keywords: solar, renewable energy, converters 


1. A survey on soft switching technique 


The efficiency of our conventional boost converters is mainly affected by the 
power loss during switching. The hard switching performed by the switches actu- 
ally causes major power lose and hard-switching of the semiconductor switches, 
generates voltage and current stress during turn-on and turn-off process. In 
order to overcome the losses during the switching stress, plenty of methods have 
been derived. Among these methods, however, soft-switching is considered to 
be the most popular and most efficient method. The soft switching technique is 
introduced in early 90s and soft switching means smooth turn on and turn off of 
switches. 

During soft switching, zero voltage switching (ZVS) and zero current 
switching (ZCS) is possible during turn-on and turn off of the power switches. 
Interleaved converters have showed a significance improvement in the effi- 
ciency when compared to the conventional boost converter which is evident 
from our previous research work. Some of the soft switching techniques 
adopted till now are summarized briefly in Table 1. This section discusses and 
reviews the performance of few soft switching topologies that were introduced 
in the recent past. 
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From the above literature, it is evident that resonant soft switching technique 
with power conditioning system could be considered as the best due to its reduced 
switching losses and its ability to attain high efficiency with no additional switches 
for soft switching. The objective of high efficiency and better performance is 
achieved by connecting PV power conditioning system (PVPCS), MPPT using PV 
simulator and current balancing control technique. Most of the above-mentioned 
work deals with soft switching with an external circuit to achieve ZVS and ZCS 
condition of the switches for boost, buck and also for buck-boost. Since switches 
experience less voltage stress across them, the conduction loss is reduced to a great 
extent which obviously improves the efficiency. 

Though, lossless switching have obtained improved efficiency care should 
be taken to reduce the converter size since extra external device are used for 
soft switching and also make the converter cost effective. Interleaved converters 
have filled this gap of over size of converters by avoiding external device rather, 
it designed and emphasized in such a way that two or more converters are put 
together to obtain improved efficiency and in turn voltage gain. 

Cha et al. [1] proposed a solar powered high efficiency soft switching boost 
converter. Soft switching is achieved by adding an additional switch of 30 kHz 
switching frequency IGBT, a couple of diodes and single capacitor to the conven- 
tional boost converter. By adopting this method, the size of system is reduced. 

Wuet al. [2] eliminated hard switching by adopting active clamping technique 
in 190 kHz switching frequency MOSFET. Where, ZCS cell is included with the 
conventional converter of 500 W, with 150 V DC input voltage. It consists of an aux- 
iliary inductor and a clamping diode. The output of 350 V DC voltages is obtained. 

Das et al. [3] achieved steep conversion ratio in the conventional converter of 
200 W and 66 kHz switching frequency with the help of coupled inductor and 
eliminated switching losses by adding auxiliary circuit consists of a switch and 
passive elements. The converter can be operated in both buck and boost mode with 
input voltage of 30 V DC and 200 V DC voltage. 

Park et al. [4] improved efficiency of our conventional converters of 700 W by 
eliminating the losses by soft switching. The converter input voltage is 230 V DC and 
output voltage is 380 V DC and has switching frequency of 30 kHz. A simple aux- 
iliary resonant circuit (SARC) adopts ZCS during turn on and ZVS during turn off. 
SARC consist of two diodes, an inductor and a capacitor. It is easy control since both 
the switches are controlled by PWM signal. MPPT controller is applied to the circuit. 

Lee et al. [5] used soft switching technique and interleaving technique in a buck 
converter to improve the step-down conversion ratio where the efficiency of the 
converter in high with an input of 200 V DC and output of 24 V DC with a switch- 
ing frequency ranges from 65 kHz to 300 kHz. Here, the active switches connected 
in series and a capacitor is connected in the power path. The two active switches 
are connected in series with a coupled capacitor in the power path to yield higher 
conversion ratio. Since switches are connected in series the voltage and current 
stress across the switches are less and easy to control. High switching frequency, less 
voltage stress and reduced switching losses are achieved. Capacitive discharging is 
lowered. High conversion ratio and smaller output current ripple is obtained. An 
output voltage ripple of 250 mV is witnessed. 

Amorndechaphon et al. [8] connected an active resonant snubber with a DC-DC 
boost converter in grid connected PV system in order to reduce the switching losses 
and switching stress. Circuit components for soft commutation are derived from 
equation of each operating stage. The proposed active resonant snubber experiences 
reduced switching losses as well as stress. 

Lee et al. [9] handled losses in a conventional boost converter of 150-340 V DC 
voltage due di/dt, dv/dt in switch with 50 kHz switching frequency and EMI with 
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Power stage Methodology Additional Achieved objective Reference 

components 

Boost Additional switch Switch, inductor, High efficiency, Gil [1] 
of 30 kHz (IGBT) capacitor couple reduced system size 

of diode 
Non-isolated Active clamp ZCS Clamping diode Circulating loss is Wu 
DC-DC boost cell (MOSFET) minimized, reduced et al. [2] 
voltage stress in 
diodes, ZCV in 
all switches, high 
efficiency 

Bidirectional Single common Coupled Step voltage Das 

DC/DC PWM auxiliary circuit inductors, single conversion ratio et al. [3] 

buck-boost active switch 

passive elements 

Boost Simple auxiliary Auxiliary switch, Reduced switching Sang- 
resonant circuit diode, resonant losses, reduced Hoon 

inductor, resonant voltage and current Park 
capacitor stress of the switching et al. [4] 
devices, easy control 

Boost Current transition Additional cell Reduced di/dt & dv/ Lee 
PWM in main switch, dt, EMI is reduced, et al. [5] 

resonant inductor current stress is 
reduced efficiency 
improved 

Boost Zero voltage DC-link in Voltage stress in Lee et al. 
switching resonant parallel switch is minimized, (2001) 

sub-harmonic 
problem due to 
DPWM is eliminated 

Boost Resonant soft Not required Reduced switching Doo- 
switching losses, high efficiency Yong-Jong 
technique, power no additional switches (1991) 
conditioning are needed for soft 
system, Interleaving switching 

Buck-boost Zero-voltage Resonant Voltage and current Doo-Yong 
switching, zero inductor, resonant stress in switches are Jung et al. 
current switching capacitor, reduced, input and (2001) 
interleaving two parasitic output current ripples 
technique capacitors and an are reduced, size and 

auxiliary switch cost is reduced 

Fly back Single capacitor Two turnoff Smooth turnoff of Chen 

snubber snubbers (with switch, with reduced et al. [6] 
inductor, system complexity 
capacitor, diode) 

Fly back Dual series resonant Clamp capacitor, Switching losses and Tseng 
active clamp auxiliary switch reverse recovery et al. [7] 
technique losses is reduced, high 

efficiency is obtained 
Table 1. 


Literature summary of soft switching technique. 


ZCT-PWM. An additional cell is connected across the main switch reduced di/dt 

obtained by resonant inductor by soft turn on of main switch. 
Doo-Yong Jung et al. [10] connected a photovoltaic power conditioning system 

to a conventional 1.2 KW ILBC for PV power generation application. The author 
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worked as an improvement in ILBC mentioned in literature in Caris et al. [11] by 
connecting PVPCS (PV power conditioning system). Interleaved resonant soft 
switching method is implemented to reduce switching losses. MPPT using PV simu- 
lator and current balancing control is done. The work achieved 97.28% efficiency. 

Chen et al. [6] three stage interleaving topologies in boost converter handles 
both voltage and current stress across the switches also get rid of the input and- 
voutput current ripples in the converter. The 600 W converter works in a couple 
of operational mode depending on the duty cycle with a switching frequency of 
50 kHz and 250 V DC input and 400 W DC output voltage. 

Tseng et al. [7] proposed an isolated interleaved buck converter (ILBC) with 
lossless turn off switches of 50 kHz using a snubber circuit. The snubber is single 
capacitive turn off snubber connected to 150 V input voltage and 100 V DC output 
voltage respectively. Work results in less number of components and complexity. 

Lee et al. [12] proposed a 400 W fly back converter which is designed in such 
a way that the series resonance help in PWM, dual series circuit does the ZCS turn 
off of the diode, active clamping is responsible for magnetizing inductance and 
ZVS turn on. Thus, high efficiency is obtained by reducing the voltage stress across 
switch with 100 kHz switching frequency hence reducing switching losses and 
reverse recovery losses. 

Table 1 summarizes the few soft switching topologies discussed in the above 
Section 1. 

Below are the few advantages of this technique that have been witnessed during 
this literature, 


The semiconductor switches are handled in a smooth fashion during its opera- 
tion to heavily bring down the switching stress. 


In addition, the switching losses (i.e.) due to di/dt and dv/dt are eradicated toa 
great extent. 


Obviously, the improved system efficiency brings down the cost of the con- 
verter system. 


It is also understood that by soft switching, EMI problems and 


Sub-harmonic problems have also been treated by the authors of the literature. 


Therefore, all benefits put together produce a good conversion efficiency to the 
converter. 

It is understood from the above literature that the soft switching is done by 
introduction of an external or auxiliary circuit in parallel to the main switch to 
make the turn on and turn off in a smooth fashion without power loss. Therefore, a 
good conversion efficiency has been obtained for the converter. 


2. Review on interleaving technique topologies 


Ultimately to extract maximum energy from the solar panel a converter is con- 
nected between the panel and the load. If high step converter is connected to the 
panel then the efficiency will be high. This includes ripple reduction with different 
conversion stages for different load conditions with different filters. A single- 
phase boost converter switches due to its extreme duty cycle faced high conduction 
losses and reverse recovery problem. The voltage stress in the semiconductor 
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switches played a major role in the reduction of efficiency. So, it is evident that 
efficiency and the cost of the converter are proportional to each other. In addition 
to that they also suffer high current ripple which does not greatly support MPPT 
algorithms. Interleaved converters have showed a significance improvement rather 
than the conventional boost converter which is justified by our previous work. 
With all the above statements, it has been concluded that the system efficiency is 
improved due to the fact that interleaving technique reduces the ripple current. 
The researchers yet have handled many various topologies and few of them are 
listed below. 

Numerous interleaving topologies have been evolved by the researchers of the 
previous decade, many such topologies have been investigated in this section and 
their methodologies are summarized in Table 1 with some cited reference. 

Garcia et al. [13] proposed Interleaved double dual boost converter IDDB of 200 
KW power and 60 V DC and 360 V DC input and output voltage respectively. It is 
a non-isolated step up DC-DC converter which gains high voltage and suits electric 
vehicle, renewable energy, high power application. Six phase IDDB is verified 
approximately. 

Tseng and Huang [14] proposed interleaved converters of 1000 W output power 
to obtain high step up gain without operating at extreme duty ratio. High step up 
converter suitable for renewable energy system with a voltage multiplier module. 
Proposed converter has reduced current stress, reduced input current ripple, and 
reduced conduction loss, extended life time of the source and improved efficiency. 
The converter has an Input voltage of 40 V, output of 380 V and efficiency of 97.1% 
is obtained. 

Lai et al. [15] included a forward energy delivering circuit and a voltage dou- 
bling circuit with the conventional ILBC to achieve step up ratio in turn high effi- 
ciency which is applied to micro grid. Here, interleaving technique is implemented 
for a high step up converter. Closed loop control of 450 W, 2 paralleled modules 
with 24 V input and 200 V output is implemented. Module efficiency up to 95.8%. 
Steady-state analysis has been done. 

Kejin et al. (2016) proposed a novel bi-directional three-phase DC/DC converter 
for high-power application. The proposed converter yields high conversion ratio 
than other conventional converters. The current and voltage stress of the power 
switches has been taken care to yield high efficiency. 

Akin [16] compared conventional boost converters and interleaved PFC boost 
converters. They are compared by their power factor (PF), total current harmonic 
distortion (THD) and total efficiency. ILBC shows an improved efficiency than 
conventional. Conventional boost converter has 97.8% efficiency and 4.88% THD 
PF ILBC 98% efficiency and 1.93% THD values. 

Table 2 summarizes the few interleaving topologies discussed in the above 
Section 2. 

Thus, from the literature on interleaving techniques we can conclude that 
interleaving technique is 


e The best solution to get rid of input and output current ripple as well as input 
and output voltage ripple. 


e Due to reduced switching loss, conduction loss and voltage stress in the 
converter switches the system efficiency is greatly improved. 


e Since the extreme duty cycle of the semiconductor switches have been 


minimized, the life time of the switches are increased and in turn resulted in 
reduced maintenance cost of the converter. 
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Power stage 


Methodology 


Achieved objective 


Reference 


Buck 


Two active switches in 
series, coupling capacitor 
in power path 


Low switching losses, improved 
step-down conversion ratio, 
reduced output current ripple 


Lee et al. [5] 


Non-isolated boost 


Interleaved double dual 


High voltage gain, suitable for 


Garcia et al. 


boost (IDDB) high power applications [13] 

Cuk and boost Comparison Cuk-design complexity, Newton 
boost-high power application, (2010) 
improved switching frequency 
reduced ripple 

Boost High step up gain without Reduced current stress, reduced Tseng and 

operating at extreme input current ripple, reduced Huang [14] 
duty ratio with voltage conduction loss, extended life 
multiplier module time of the source and improved 
efficiency of 971% achieved 
Boost, interleaved Comparison PFC ILBC-improved power Akin [16] 


power factor 


factor (PF), reduced total current 


correction boost harmonic distortion (THD) of 
(PFC ILBC) nearly 3%, improved efficiency of 
0.2% than conventional converter 
Boost Interleaving, adding High step up ratio in turn high Lai et al. [15] 
forward energy delivering efficiency of 95.8% and also 
circuit and a voltage applied to Micro grid 
doubling circuit to ILBC 
Boost Adding coupled inductor Reduces the current stress, Tseng and 
to ILBC current ripple and conduction Huang [14] 
losses 
Boost Bi-directional three-phase High conversion ratio, current Jin and Liu 
and voltage stress in switches are [17] 
reduced 
Table 2. 


Literature summary of interleaving technique. 


e Interleaved converters are superior due to its simple structure and low 


complexity. 


e From the above literature, it is also evident that they can serve for wide range of 


applications. 


Though many techniques have been evolved and implemented in the boost 
converter still lag exists in the research in extracting 100% energy from the solar 
cell. Interleaving technique being one of the best solutions to get ripple free boosted 
voltage which results in improved efficiency. But, above all care should be taken in 
load matching of the converters. 


3. Review on special converters 


In literature, researchers have compared SEPIC converter, ZETA, Cuk converters 
with ILBC. SEPIC converter means single ended primary inductance converter and 


ZETA converter. 


Newton [18] applied interleaving concept to the Cuk converter and analysis of 
system advantages over ILBC is compared. Cuk converter provides isolation without 
addition of extra active devices. Isolation gives different input and output facility. 
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Rech et al. (2001) compared different hybrid converters and new hybrid con- 
verter is proposed for high power application. The number of level of power conver- 
sion is reduced and THD is also reduced. 

Yan Li et al. (2014) suggested buck-boost mode is preferred in applications 
where the input current needs to continuous and the output with less ripples. This 
gives a long-extended battery life. Non-isolated push pull converters (NIPPC’s) are 
compared. Advantages of NIPPC’s are less current stress, high efficiency and good 
EMC performance. Different types of configuration of push pull converters with 
single and coupled inductors are compared. 

Kosai et al. [19] proposed a silicon carbide (SiC) Power semiconductor tech- 
nology results in SiC rectifiers and controlled semiconductor devices (JFET and 
MOSFET), where the power converters can be operated in 200°C. 

Qahoug et al. [20] proposed that by using only a single output voltage value, the 
controller can track the MPP of every solar panel with improved tracking speed, 
cost effective and also reduced system size. This technique works efficiently on 
non-uniform loads and mismatching conditions due to partial shading with single 
voltage sensor. 

Darwish et al. [21] discusses about the inherent non-linearity of Cuk converter 
which is the reason for distortion of output voltage and current. A three phase Cuk 
converter is proposed. The buck-boost nature of the converter provides flexibility 
for standalone and grid connected applications. 

Hong-Tzer et al. [22] introduced a single phase Cuk bridgeless AC-DC power 
factor correction rectifier which produces positive output voltage. This converter 
converts negative output voltage to positive value unlike other traditional topolo- 
gies. The efficiency is improved by reducing the conduction loss of the power 
switches and also the power factor is maintained at 0.99 all input and output 
conditions. 

Yifannyu et al. (2016) optimized the power imbalance problem in PV system 
by zero injection sequence injection method. They also introduced and applied the 
above-mentioned method in a multilevel-cascaded h-bridge converter to test the 
feasibility and effectiveness of the introduced method. 

Forest et al. [23] is implemented isolating intercell transformer (ICT) in 
an interleaved converter. Two ways are implemented with coupled inductor or 
transformer. One is bidirectional converter buck or boost, another is bidirectional 
converter buck-boost. Advantages of this concept are removal of discrete inductors 
due which magnetic loss is less, high step up ratio and minimization of switching 
losses. 

Khan et al. [24] proposed a synchronous buck converter with high conversion 
ratio and reduced system size. Gallium nitride (GaN) field-effect transistors are 
used in the converter. MPPT is implemented for all irradiation conditions where, 
maximum power extraction is achieved. 

Kaouane et al. [25] obtained maximum energy transfer from generator to load 
side then other conventional converters by designing a SEPIC converter. High 
efficiency is obtained also during discontinuous conduction mode using PWM 
controlled switches. 

Table 3 summarizes the few special converters involved in PV Power Generation 
discussed in the above Section 3. 

ILBC yields better efficiency with simple design structure which serves for 
wide range of applications. Moreover, it is evident that it suffered with less ripples 
when compared the rest of the two converters with increased switching frequency. 
And also, they are compact in size though additional components are required for 
the control of switches. Finally, this section dealt with some other type of special 
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Converter Methodology Summary Results Reference 
Cuk converter Comparison of Isolation gives Cuk converter Singh et al. 
Cuk and ILBC different input and provides isolation (2005) 
output facility without addition of 
extra active devices 
Hybrid Comparison Applied in high Number of level is Rech et al. 
converters different hybrid power application reduced and THD is [26] 
converters and reduced 
new hybrid 
converter 
Non-isolated Push pull Buck-boost mode Long extended Yan Li et al. 
push pull converters is preferred in battery life, less (2014) 
converters with single applications where current stress, high 
(NIPPC’s) and coupled the input current efficiency and good 
inductors are needs to continuous EMC performance 
compared and the output with 
less ripples 
MPPT controller Single output Woks efficiently in Improved tracking Jaber Abu 
voltage value, non-uniform loads speed, cost effective Qahouq 
the controller and mismatching and Reduced system et al. [20] 
can track the conditions due to size 
MPP partial shading with 
single voltage sensor 
Cuk converter Buck-boost To overcome Provides flexibility Yang 
nature inherent non- for standalone and et al. [22] 
linearity in turn grid connected 
distortion of output applications 
voltage and current 
Single phase bridgeless Converts negative Improved efficiency Yu 
Cuk AC-DC rectifier output voltage to reduced conduction et al. [27] 
positive value unlike losses, power factor 
is maintained at 0.99 
all input and output 
conditions 
Multilevel- Zero injection To optimize the Reduced switching Khan 
cascaded sequence power imbalance loss, high conversion et al. [24] 
h-bridge injection problem in PV efficiency, simple 
converter method system structure 
Synchronous Gallium nitride For all irradiation High conversion ratio Kaouane 
buck converter (GaN) field- conditions where, and reduced system et al. [25] 
effect transistors maximum power size 
are used in the extraction is 
converter achieved 
SEPIC converter DCM mode — Maximum energy Simonetti 
using PWM transferred from etal 
controlled generator to load 
switches side, High efficiency 
Interleaved Isolating Two ways are Removal of discrete Forest 
converter intercell implemented with inductors due which François 
transformer coupled inductor or magnetic loss is less, etal. [23] 
(ICT) transformer (buck or high step up ratio 
boost, bidirectional and minimization of 
buck-boost) switching losses 
Table 3. 


Literature summary special converters involved in PV power generation. 
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converters like Cuk, SEPIC, etc. Here some work on comparison of ILBC’s with 
SEPIC and ZETA converter, pros and cons between them, some new techniques 
implemented in special converters are also discussed. Table 3 summarized the 
aforementioned features. 


4. Review on MPPT technique 


As we know PV panel is a non-linear device and in order to be able to extract 
maximum output power despite changing insolation due to partial shading and 
changing temperature due to weather conditions, maximum power point tracking 
(MPPT) is required. However, MPPT works based on maximum power transfer the- 
orem which goes like this, whenever source impedance is equal to the load imped- 
ance maximum power is delivered to the load. This can be done by adjusting the 
duty cycle of a DC-DC converter between the source and the load. The researchers 
have approached various control techniques to achieve this that is listed in Table 4. 

Doo-Yong-Jong et al. (2012) imposed MPPT to half bridge PV DC-DC boost 
converter and full bridge inverter applied to grid which adopts variable step-size 
algorithm. Grid current control is done with 4th order linear IIR filter controller. 

Schuck and Pilawa-Podgurski [28] attended the problem faced by PV input 
interleaved converter with MPPT to work in asymmetric fashion due the difference 
in load and source combinations. The author proposes a control technique based 
on and harmonic elimination and ripple reduction under asymmetric source-load 
conditions. 

Femia et al. [29] proposed distributed MPPT rather than the conventional MPPT 
to overcome the problem of model mismatch and partial shading in PV modules. 
Though the author showed improvement in the efficiency is a research gap of 
complexity in the design. 

Qin et al. [30] formed an algorithm to extract maximum power from PV mod- 
ules with differential power processing converters. Fast tracking is achieved with 
less number of perturbations for long strings. 

Park et al. [4] designed a MPPT controller to track maximum power from PV 
module during fast transient changes due to environmental conditions with good 
tracking speed and accuracy. The author analyzed different dynamic PV panel 
data with high number of samples than before which includes different whether 
conditions. 

Pilawa-Podgurski and Perreault [31] proposed an MPPT power converter 
with digital control technique for Synchronous buck converter. Obtains local and 
global MPPT in parallel connected PV modules. Direct integration is done hence 
economical. 

Qin and Pilawa-Podgurski [32] proposed a differential power processing archi- 
tecture. MPPT with minimum communication and no local sensing is ensured. 
Partial shading condition is considered. The system ensures extraction of more 
energy which is employed low voltage buck-boost converter. Control method 
achieves global and local MPPT operation. 

Alexander et al. [33] analyzed different MPPT algorithms. Presence of noise 
in MPPT algorithm leads to deviation MPP steady state. Parameters like slowing 
tracking and noise are considered which affect the overall efficiency of the system. 
Optimizations of these system parameters are determined. 

Sundareswaran et al. [34] formed artificial bee colony (ABC) algorithm for MPP 
global maximum power point (GMPP) which enables maximum utilization of solar 
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System Problem Methodology Objective Reference 
achieved 

Boost, Unbalance in Repetitive current Minimal devices Doo-Yong-Jong 

half-bridge grid current control method- are used, regulates et al. (2012) 


connected to 


grid 


current controller 
(repetitive) -based 
on fourth order 
linear phase IIR 
filter. Variable step- 
size technique 


the grid current 


Asymmetric Asymmetric Harmonic Reduced Schuck and 
Interleaved conditions due elimination control Ripples, works Pilawa-Podgurski 
Multiphase to difference technique in asymmetric [28] 
in source condition, first 
and load harmonic ripple 
conditions component is 
reduced by 14.8 dB 
Interleaved Effect of Distributed MPPT Improved Femia et al. [29] 
model efficiency, reliable 
mismatch and 
partial shading 
in PV module 
Series To maximize Differential Fast tracking Qin et al. [30] 
connected the power power processing without extra 
extraction algorithm hardware, 
from series reduced no. of 
connected PV perturbations in 
module each step, suitable 
for long sub 
module strings 
DC-DC, Rapid MPPT controller Improved Park et al. [4] 
DC-AC transient design considering tracking speed 
changes in PV field measurements and accuracy— 
output power of transient effects captures dynamic 
PV panel data of 
high temporal 
resolution 
Synchronous To increase Sub module Ensures both local Pilawa-Podgurski 
buck converter of the energy integrated digital and global MPPT, and Perreault [31] 
captured control technique direct integration 
during partial is done hence 
shading economical 
conditions. 
To reduce the 
overall cost 
Buck-boost extraction of Differential Global and local Qin and Pilawa- 
converter more energy power processing MPPT is obtained Podgurski [32] 
from PV architecture (MPPT 
applications with minimum 
communication and 
no local sensing) 
PV power MPPT under ABC algorithm Ensures fast Alexander 
generation partial shading tracking et al. [33] 
conditions characteristics, 


improved energy 
saving capability, 
improved revenue 
generation 
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System Problem Methodology Objective Reference 
achieved 
Bidirectional MPPT for Parallel MPPT Reduced loss in Sundareswaran 
DC-DC standalone technique power converter, etal. [34] 
converter power improved 
generation efficiency, reduced 
number of power 
processing stages 
Standalone To achieve Neural-Network- Fast and improved Subudhi and 
hybrid system fast, better Based MPPT transient response Pradhan [35] 
(solar and transient control, radial basis 
wind) response and function network 
pitch angle of (RBEN), Elman 
wind turbine Neural Network 
(ENN) 
PV system To minimize FPGA platform to Settling time is Roger Gules et al. 
the settling obtain real time minimized (2008) 
time. weather conditions, 
optimization of 
perturbation period 
in P&O algorithm, 
the settling time is 
calculated by dual 
kalman filter 
Buck converter To improve Fuzzy-PID Achieving high Whei-Min Lin 
solar energy controller MPP under various Member et al. 
conversion environmental (2011) 
efficiency conditions 
PV system To obtain Two phase tracking Maximum power Jiang et al. [36] 
maximum INC MPPT transfer 
power transfer 
to load 
PV system To attain GMPPT technique Less design Sarah and 
GMPP with a simulated complexity, Haque [37] 
annealing (SA) less number of 
parameters per 
iteration in PXO 
algorithm 
Grid To identify Modern predictive Dynamic response Shadmand 
connected PV maximum control (MPC) of the system is et al. [38] 
system power point improved 
(MPP) rapidly 
PV system Uncertainty Double integral Easy to control Raseswari 
due to rapidly sliding mode where, switching et al. [39] 
changing controller (DISMC) frequency is kept 
weather MPPT system constant 
condition and 
variable load 
condition 
PV system To avoid Gray wolf Better transient Seyed 
steady-state optimization response Mahmoudian 
oscillation, method for MPPT etal. 
to improve tracking 
tracking 
efficiency 
Table 4. 


MPPT techniques imposed in solar powered converters. 


137 


Advanced Statistical Modeling, Forecasting, and Fault Detection in Renewable Energy Systems 


energy under varying atmospheric conditions. Two different PV configurations 
under different shading patterns are tested. This ensures fast tracking characteris- 
tics, improves energy saving capability, increased revenue generation when com- 
pared with other energy MPPT techniques. 

Subudhi and Pradhan [35] reviewed MPPT techniques developed until January 
2012. Analysis based number of control variables, types of control strategy, types of 
circuitry for PV systems and application based on practical and commercial. Nearly 
26 variety of MPPT technique has been compared. 

Gules et al. [40] proposed a bidirectional DC-DC Converter with parallel MPPT 
for standalone power generation. The power converter works as multifunctional 
MPPT circuit, battery charger, battery regulator and as step up converter. Parallel 
MPPT helps in reducing losses in the power converter in turns improves efficiency 
also reduces the number of power processing stages. Same series connected MPPT 
algorithms can be used. 

Lin et al. [41] constructed a standalone hybrid system consists of solar power, 
wind power, diesel engine and intelligent power controller. Radial basis function 
network (RBFN) is used to achieve fast and transient response and pitch angle of 
wind turbine and Elman Neural Network (ENN) is used to achieve MPPT. 

Jiang et al. [36] proposed a grid connected boost half bridge PV micro inverter 
system which uses only minimal device. Dynamic response is improved when load 
and solar irradiance were changed rapidly. Variable step size MPPT is adopted. 

Ricco et al. [42] used FPGA platform to MPPT algorithms acquire its parameters 
to the real-time weather conditions. The settling time is minimized. The proposed 
controller presents optimizing the perturbation period of the P&O algorithm the 
settling time can be calculated by a dual kalman filter. 

Dounis et al. [43] proposed that Maximum power point tracking (MPPT) with 
Fuzzy-PID controller to regulate to output power of PV system for a DC/DC buck 
converter. The proposed system improves the solar energy conversion efficiency. 
Fuzzy-PID is highly effective in achieving high MPP under various environmental 
conditions. 

Hsieh et al. [44] proposed two phase tracking INC MPPT which either constant- 
frequency with variable-duty cycle (CFVD) or variable-frequency with constant 
duty control(VFCD) to perform with reference to a commonly known threshold- 
tracking zone (TTZ) thus obtained maximum the power transfer to load at MPP. 

Wang et al. [45] identified that the mismatches the PV system affects the major 
portion of energy production. Then, proposed a MPPT system which is simple with 
unified output. 

Lyden et al. (2015) proposed that the global maximum power point track- 
ing (GMPPT) technique with a simulated annealing (SA) for PV energy system 
which experiences the partial shading condition (PSC). This attains GMPP or near 
GMPP. Less design complexity when compared to P&O technique and less number 
of parameters per iteration to iteration. Shadmand et al. [38] proposed a MPPT 
technique called as modern predictive control (MPC) for a grid connected PV 
system. The technique identifies the Maximum Power Point (MPP) rapidly than 
other earlier techniques. The dynamic response of the system is improved than 
other conventional techniques. 

Pradhan and Subudhi [39] handled the uncertainty problem due to rapidly 
changing weather condition and variable load condition in a PV system by imple- 
menting double integral sliding model controller (DISMC) MPPT system. The 
system is easy to control where, switching frequency is kept constant. Mohanty 
et al. [46] used Gray wolf optimization method for MPPT tracking. The limitations 
such as steady-state oscillation, lower tracking efficiency and transient condition 
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has been fixed in this method. Where, our other traditional methods such as P&O, 
improved P&O failed to do so. 

Table 4 summarizes the few MPPT techniques imposed in solar powered 
converters discussed in the above Section 4. 

It is clearly evident from the above literature review that various problems have 
been attended and solved by implementing MPPT controller to a system. To sum- 
marize some of the key problems addressed, 


e Uncertainty due to partial shading and changing weather conditions. 


e Asymmetric source-load mismatch, model mismatch especially in interleaved 
converters. 


e Rapid transient changes in PV power generation system. 
e Steady state oscillations are avoided by minimizing the settling time. 
e Unbalanced grid current is handled. 


The MPPT techniques used by the researchers in the previous decades have been 
discussed in the above section. The Problem, objectives, various techniques and 
advantages of MPPT imposed PV system are summarized and discussed in Table 4. 
The MPPT applied PV power generation system have witnessed the following 
advantages, which is understood from the literature. 


e High conversion efficiency because of maximum power extraction especially in 
series connected models. 


e Reduced ripples and harmonic distortions in the system. 
e Fast maximum power tracking in turn better tracking efficiency. 


Supports wide range of applications, which includes stand-alone power genera- 
tion and also in pitch angle control for wind turbines. 


5. Review on controller for ILBC 


Since the boost converters are nonlinear in nature, the controller maintains the 
stability during the source-load variation. In order to balance the fast-transient 
response, it is vital to select an appropriate controller to overcome over shoot damp- 
ing and achieve the steady state response. Different types of controllers used by 
the researchers so far till 2016 have been discussed here. The extended summary is 
displayed in Table 5. 

Zdravko et al. [51] proposed a sensor less multiphase average current controller 
which a self-tuning current estimator. The average inductor current is calculated 
with the duty cycle, input and output currents, based on these parameters the active 
phase gets automatically changed. 

Cheng et al. [47] introduced Hybrid ripple concept for ripple-based control 
which may suffer from large jitter issues and instability. It is achieved by feed back- 
ing both inductor current ripple and compensated external ramp to the modulator. 
Reduce jitter and improved stability is achieved. Fast load transient performance 
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with the proposed control scheme of controller is achieved. Research gap: system 
becomes more complicated with two different ramp compensations. Wang and Li 
[48] developed a RPC (ratio pre-setting controller) method to enhance current 


Controller Technique Methodology Summary Results Reference 
Multiphase Self-tuning Average Using these Sensor less ZdravkoLukié 
average current inductor parameters controller (2009) 
current estimator current is active 
controller calculated phase gets 
with the duty automatically 
cycle, input changed 
and output 
currents 
Ripple- Hybrid Feed backing Research Reduce Cheng [47] 
based ripple both Inductor gap-system jitter and 
control concept current becomes more improved 
ripple and complicated stability. 
compensated with two 
external different ramp 
ramp to the compensations 
modulator 
Ratio To enhance The Three types of Y-connection Wang and Li [48] 
pre-setting current relationship three-phase type 
controller sharing between bidirectional converter 
(RPC) capability phase shift DC-DC evidence 
angle and the converters are less current 
transformer introduced mismatch 
current is and other due to its 
arrived two types are current 
compared sharing 
capability 
Power Artificial Output Battery as load Fast settling Jiteurtragool 
control Neural voltage is time and et al. [49] 
system Network regulated low over 
(ANN) which is shooting 
received 
from the PV 
array from 
changing 
ambient 
temperature 
and irradiance 
IC IC is Controls the 4 Phase Use less Pulvirenti 
controller embedded converter as ILBC for PV number of et al. [50] 
with power well as takes distributed external 
MOS care of MPPT power components 
algorithm conversion and less 
overall cost 
State State Based on Method Good Rafael Pena- 
trajectory trajectory boundary involves four dynamic Alzola (2015) 
controller control control switching response 
states, and during 
three state starting and 
equations load change 
Table 5. 


Literature summary of analysis of controllers and implemented applications for PV generation system. 
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sharing capability. Here, three types of three-phase bidirectional DC-DC convert- 
ers are introduced and other two types are compared. Where, Y-connection type 
converter evidence less current mismatch due to its current sharing capability. The 
relationship between phase shift angle and the transformer current is arrived. 

Jiteurtragool et al. [49] proposed a power control system for DC-DC converter 
using ANN. Output voltage is regulated which is received from the PV array from 
changing ambient temperature and irradiance. With battery as load hardware 
implemented with fast settling time and low over shooting. Pulvirenti et al. [50] 
proposes IC controller 4 Phase ILBC for PV distributed power conversion. The 
system uses less number of external components and less overall cost is assured. 
An IC is embedded with power MOS which controls the converter as well as takes 
care of MPPT algorithm. Sean and Williamson [52] compared various state space 
model for power amplifiers and concluded there is a considerable difference in the 
measured frequency response. Three different models are compared in this chapter. 
Rafael Pena-Alzola et al. (2015) handled instability problems in ILBC’s during large 
bandwidth requirements by state trajectory control which is based on boundary 
control where, conventional controller failed to do so. The system provides a good 
dynamic response during starting and load change. The method involves four 
switching states, and three state equations. 

Table 5 summarizes the few analysis of controllers and implemented applica- 
tions for PV Generation system discussed in Section 5. 

Implementation of controller is also being a best method to eliminate ripple. Not 
only for ripple reduction, the controller also helped to improve the systems overall 
performance by maintaining good dynamic response. 


6. Ripple minimization technique 


When input current is equal to the inductor current then the converter works in 
continuous conduction mode (CCM). But the drawback in this style of working is 
the diode reverse recovery current. When the inductor current does not match with 
input current then is said to be working in pulsating mode. In this mode, the input 
and inductor currents are equal only when it is switched ON. Boost converter works 
in CCM mode but buck converter works in pulsating mode. So, in pulsating mode of 
operation the inductor current turns zero when the switch is OFF. Coupled induc- 
tor can bea solution for this issue but then integrated magnetic components have 
caused EMI problems so coupling coefficient should be designed well in advance. 
So, pulsating type interleaving technique can be adopted to achieve current ripple 
reduction. 

Kim [53] proposed reconfigurable solar converter (RSC) which helps in single 
stage power conversion. The proposed system is nothing but a single stage, 3 phase 
grid tie solar PV converter for DC-DC or DC-AC generation used for PV battery 
application. Numbers of conversion stages are minimized. Therefore, improved 
efficiency, reduced cost, weight and volume are ensured. 

Garcia et al. [54] handled the tolerance in the inductance on the output current 
and the voltage ripple. The author concluded that Multiphase is advantageous for 
its filter less input and output. Where, in the other side tolerance in semiconductor 
switches is undesirable. 

Caris et al. [11] proposed a method to calculate the phase-shifts to avoid harmonics 
generated due to the cell inductors in the interleaved converters. Interleaved topolo- 
gies used to improve power converters to obtain higher power ratings, higher band- 
width and lower current ripple. On other hand, it suffers with cell inductor tolerance. 
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Yu Gu et al. (2013) proposed ILBC with RPN (ripple cancelation network). 
RPN is connected with a normal ILB, it consists of additional inductors, a couple of 
capacitors and coupled inductor. Input ripple current cancelation with minimum 
loss and 96% efficiency is achieved. 

Caracas et al. [55] implemented solar water pump where, a 3-phase induction motor 
runs directly from the PV cell without the help of battery. This is a system based on reso- 
nant two inductor boost converters (current resonant converter) and full bridge 3 phase 
voltage source inverter. Also, non-isolated recovery snubber is added with the system 
in order to control hysteresis loss and duty cycle improvement in order to improve the 
efficiency. Advantage of resonant two inductor boost converters is high voltage gain and 
low input current ripple. Prolonged lifetime due to the absence of battery, use of induc- 


tion motor is more efficient, maintenance free when compared with DC motor. 


Peng [56] ILBC with self-balancing capability is designed. Smoother and less 
distorted power conversion is achieved. More than two level will be convincing 
to balance DC voltage level without additional circuits. Advantages of system are 
magnetic less, compact, high efficient, zero EMI, low cost power conversion. 

Ramos-Paja et al. [57] analyzed the pre-filter optimal operatinvg conditions 
where the input ripple is mitigated for better conversion ratio of ILBC. The inductor 
losses are predicted in terms of efficiency and input current ripple then a pre-filter 
is designed accordingly. 


Problem Methodology Results Remarks Reference 
Load side ripple Minimized by Research Can be applied to Kim et al. 
in parallel increasing the gap-increase in smart grids [53] 
multiple in input switching frequency switching loss 
source 
More number of Reconfigurable solar Single stage Improved efficiency, Garcia 
conversion stage converter (RSC) conversion reduced cost, weight et al. [54] 
and volume are 
ensured 
Tolerance in Digital control, Filter less input Research gap- Caris et al. 
the inductance multiphase and output undesirable side [11] 
on the output tolerance in 
current and the semiconductor 
voltage ripple switches 
Harmonics Calculated phase shift Higher Guand 
in multilevel bandwidth, Zhang 
converter ow current [60] 
ripple 
Input ripple Ripple cancelation Minimum Consists of additional Peng [56] 
current network (RPN) oss and 96% inductors, a couple 
efficiency of capacitors and 
coupled inductor 
Power mismatch ILBC with self- More smoother Magnetic less, Ramos- 
balancing capability and less compact, high Paja et al. 
distorted power efficient, zero EMI, [57] 
conversion low cost power 
conversion 
Input current Pre-filter optimal Better Jiang et al. 
ripple mitigation operating conditions, conversion [59] 
predicted inductor ratio in ILBC 


losses 


Table 6. 


Literature summary of ripple minimization techniques in PV power generation. 


142 


Survey on Photo-Voltaic Powered Interleaved Converter System 
DOI: http://dx.doi.org/10.5772/intechopen.90111 


Tseng et al. [58] included voltage multiplier for managing asymmetrical inter- 
leaving condition without changing the duty cycle extremely is proposed. Voltage 
multiplier is nothing but conventional converter with a coupled inductor. Reduces 
the current stress, current ripple and conduction losses. 

Jiang and Wang [59] discussed about the current ripples in converters and pro- 
posed of three phase PWM converter. The author concluded bigger current ripples 
are generated by discontinuous PWM when compared to space vector PWM. This 
is mostly suitable for variable switching-frequency control for three phase PWM 
converters applications. 

Table 6 summarizes the few ripple minimization techniques in PV power 
generation discussed in the above Section 6. 

There is high frequency ripple caused by the switching frequency and the low 
frequency ripple caused by AC load. For the symmetrical 3-phase load, low fre- 
quency ripple is not experienced. Rather, the single-phase or unbalanced 3-phase, 
will experience a double line frequency ripple. 

Mainly in interleaving converters source and load mismatch in various conver- 
sion stages have occurred. This is handled by the researchers in two ways either, 


e By tunning the switching style or, 


e By some unique ripple mitigation methods and fixed it with some additional 
integrated circuit. 


In the above two methods the later suffers with switching stress and the former 
with increased volume of the system size. This method can be considered as bridge 
to fill the research gap when interleaving technique cannot be applied and consid- 
ered as limitation. This section discusses about the ripple minimization techniques 
used so far and Table 6 exhibits the summary of the same. 


7. Power conditioning PV simulator and performance evaluation 


Below are some set of references which carried few sophisticated tools for bet- 
terment of the PV system that involves 


¢ Evaluation 
e monitoring 
e As well as rectification 


Koran et al. [61] proposed high efficiency and fast transient response time that 
are essential for PV simulator. PV source simulator-evaluate the power conditioning 
system and MPPT algorithms, high efficiency photovoltaic source simulator both 
analog and digital based is proposed. 

Kim et al. [62] compared several parameters like efficiency and power factor 
performance. Several topologies are compared, conventional average current 
mode control boost PFC, an interleaved boost PFC, a back-to-back bridge- 
less boost PFC, semi-bridgeless boost PFC, are assessed through loss analysis. 
Improved power factor correction (PFC) topologies suitable for a high density are 
identified. 

Singh et al. [63] improved power quality by a control strategy for achieving 
maximum benefits for grid interfacing inverters. Perform a multi-function device 
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by incorporating active power filter functionality. Power converter to inject power 
generated from RES to the grid. 

Rosario et al. (2005) analyzed series—parallel L, Cp, Cs resonant inverter for 
effect of tolerance, to ensure circuit performance-with no feedback-good repeat- 
ability. The author concluded that the capacitor value affect the sensitivity values 
and done statistical study by Monte-Carlo method. 

Roh et al. [64] proposed two phase interleaved critical conduction mode (CRM) 
power factor correction for boost converter. Variation tolerant phase shifter (VTPS) 
which ensures accurate 180-degree phase during interleaving operation. This VTPC 
circuit can be applied to any type of ILBC’s. 

Stanth et al. (2013) proposed a converter which is connected in parallel with the 
PV string to improve the energy captured during the shading or power mismatch. 
The parallel configuration is proposed to handle mismatch of power, its turns off if 
no mismatch. The converter handles only, mismatch in the adjacent string, turns off 
during no mismatch, low average power handling and low operational duty cycle 
these can be considered as research gap. 

Oscar et al. [65] developed a multiphase converter with buck converter as base. 
Natural ZVS is achieved in buck converter the instantaneous inductor current itself 
will manage the current imbalance in the phases. Better current balance and ZVS 
in both transitions are achieved but with high current ripple. Higher conduction 
losses, high number of phase and high current ripple are found as research gap. 

Liserre et al. [66] introduced generalized hybrid modulation (GHM) technique 
to overcome the drawback in the conventional method that these techniques 
succeed in canceling the sideband harmonics in multilevel voltage only when the 
DC-link voltages are equal. PSC-PWM technique proposed, a method to obtain the 
proper shifting angles to avoid the reduction in overall output voltage. 

Shibin Qin Andrew et al. (2014) proposed a differential power processing 
(DPP) which combines micro inverter and DPP converter. The architecture con- 
figures DC-DC parallel with PV string. Power yield is improved by overcoming the 
mismatch in sub module. Since the converters are connected in parallel low power 
rating is sufficient. Low power losses are evident and suitable for sub module 
MPPT. 

Tolbert et al. [67] introduced a procedure to find all sets of switching times 
(angles) for each converter in a multilevel converter to obtain desired output 
voltage. Here, 5th and 7th are eliminated by converting the transcendental equa- 
tions that identifies the elimination of the harmonics into a set of polynomial 
equations. 

Zhang and Yu [68] investigated on the feature of interleaved PWM DC-DC 
converter and analyzed their methodology. Practical design is presented such as 
ripple quantification, zero ripple necessary condition and sufficient condition, fault 
detection, and waveform shape prediction. 

Sangswang and Nwankpa [69] quantified random noise in PWM switching 
DC-DC boost converter. The presence of noise creates uncertainty in switching 
time. Performance analysis for the parameters such as stochastic performance 
index, MPPT and critical energy are done. 

Zhang et al. [70] connected pair of additional diode and capacitor/inductor 
are connected in X-shape to a conventional converter (diode assisted converter). 
Comparison based on silicon devices, passive components, EMI and efficiency are 
done. Input current ripple is reduced, improved power conversion and buck-boost 
capability are ensured. 

Ray et al. [71] done a multimode analysis of the ILBC for switching duty cycle 
from 0 to 1, for both CCM and two DCM modes. This analysis helps designer to 
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optimize the inductor value and coupling factor for different performance values. 
Also, could optimize the number of power stages for varying input and loads. DCM 
shows less switching losses. 

Li and He [72] reviewed on challenges based on high step up, low cost, high 
efficiency converters are summarized. To extend the voltage gain, to avoid extreme 
duty cycle in order to reduce current ripple, to avoid power devices losses by soft 
switching and to avoid reverse recovery problem of diode are the problems summa- 
rized. Limitations of conventional boost converters are analyzed. Parallel connected 
and non-isolated high step up converters for its low cost and improved efficiency. 

Lefevre and Mollov [73] discussed about the asymmetric voltage operation while 
handling multiphase boost converter. Conversion efficiency of 99% is achieved by 
introducing a flying-capacitor in the existing topology. 

Nigsch et al. [74] proposed a bridgeless single stage PFC converter for efficient 
solution to fulfill the harmonic standards. Efficiency is improved by removing the 
input rectifier and providing galvanic isolation by allowing single switch active at a 
time. Around 95% of efficiency is obtained and single stage approach makes it a low 
cost high-power density PFC. 

Xiong Li et al. (2016) proposed a space vector pulse-width modulation 
(SVPWM) to regenerate phase current using neutral-point current sensing, 
scheme. Complementary active vectors replace the zero and offset vectors. Hybrid 
switching pattern are used, to synthesize the reference voltage vector. Phase cur- 
rents can be reconstructed using SVPWM scheme. 

Yuan-Chuh et al. (2016) achieved desired generator performance by one- 
cycle voltage regulation (OCVR) for permanent-magnet synchronous generator 
(SPMSG). Design of voltage controller is simplified and computational time is less 
when compared to other conventional controllers. 

Table 7 summarizes the few power conditioning in PV power generation 
discussed in the above Section 7. 

The power quality monitoring has been done to enhance the system perfor- 
mance, the power quality improvement includes the following, 


e Power factor correction 

e Power mismatch in multiphase converters 

° Fine-tunned control strategy for maximum benefits 

e Optimized MPPT algorithms 

These set of dynamic response tuning of the PV system are overviewed in this 
section and power quality issues are also discussed here. The power conditioning, 
power quality improvement and performance evaluation works of the previous 
decade has been discussed in the above section and the objective, evaluation param- 
eters, summary of discussions are summarized in Table 7. 
8. Inference 

This research work has compared and identified an efficient and suitable 
converter extract solar energy from the PV panel in an effective way to obtain 


high conversion ratio among the conventional converters. The overview exhibited 
and identified the requirements of a best PV power generation system, where the 
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Objective Evaluation Summary Remarks Reference 
parameter/ 
methodology 
PV source Efficiency and Both analog Improved Koran et al. [61] 
simulator transient response and digital efficiency, 
based are fast transient 
proposed response 
Comparison, Loss analysis of Efficiency and Kim et al. [62] 
power factor converters—current power factor 
correction mode control boost, correction 
interleaved boost, 
buck-boost and semi 
bridge boost 
Power quality Performs asa Control Inject power Mukhtiar Singh 
improvement multifunctional strategy from RES to et al. [63] 
device by adding an to achieve grid 
active power filter maximum 
benefits 
from grid 
interfacing 
inverter 
Output current Effect of tolerance, Capacitor Good Rosario et al. 
sensitivity performance without values effect repeatability (2005) 
analysis feedback the sensitivity 
value 
Power factor Variation tolerant Two phase Accurate Roh et al. [64] 
correction phase shifting interleaved 180-degree 
technique (VTPS), converters phase shift 
critical conduction and can be 
mode technique applied to any 
(CRM) type of ILBC’s 
Improvement Resonant switched Converter Research gap- Stanth et al. 
of dynamic capacitor converter imposed in handles only (2013) 
response parallel with adjacent string 
PV string mismatch, 
turns off 
during 
mismatch, low 
average power 
handling 
and low 
operational 
duty cycle 
Current Current self-balance Multiphase Research Oscar et al. [65] 
imbalance in the mechanism in CCM buck gap-high 
phases mode converter conduction 
losses, more 
no of phases, 
high current 
ripple 
Sideband Generalized hybrid To obtain Reduction in Liserre et al. [66] 
harmonics modulation technique, exact phase overall output 
in multilevel carrier PWM shifting angle power is 
voltage, technique avoided 


multilevel phase 
shifting in non- 
equal DC link 
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Objective Evaluation Summary Remarks Reference 
parameter/ 
methodology 
Mismatch in Differential power Combines Less power Shibin Qin, 
sub-module processing (DPP) micro inverter is only Andrew et al. 
DC-DC parallel and DPP required since (2014) 
PV string converter connected in 
parallel 
Harmonic To find switching In multilevel 5th and 7th Tolbert et al. [67] 
elimination times (angles) converters harmonics are 
where, eliminated 
non-equal DC 
source 
Ripple Analysis of interleaved DC-DC Ripple Zhang and Yu [68] 
quantification, DC-DC PWM converter minimization, 
zero ripple converter fault 
necessary detection, 
condition wave shape 
prediction 
Uncertainty Analysis of critical PWM Anawach 
switching time energy, stochastic switching Sangswang and 
due to noise performance index DC-DC boost Nwankpa [69] 
converter 
To optimize Multimode analysis of Multimode Helps to Ray et al. [71] 
the number of the ILBC for switching analysis of the optimize the 
power stages for duty cycle ILBC inductor value 
varying input and coupling 
and loads factor for 
different 
performance 
values 
Asymmetric Flying capacitor is Multiphase 99% Liand He [72] 
voltage operation introduced boost conversion 
converter efficiency 
To fulfill Removing the input Bridgeless 95% of Guillaume Lefevre 
harmonic rectifier and providing single stage efficiency is et al. (2010) 
standards galvanic isolation by PFC converter obtained, low 
allowing single switch cost high- 
active at atime power density 
PFC 
To synthesize the Space vector pulse- — — Nigsch et al. [74] 
reference voltage width modulation 
vector (SVPWM) 
Performance One-cycle voltage Permanent- Simplified Liet al. [75] 
enhancement regulation (OCVR) magnet design of 
synchronous voltage 
generator controller 
(SPMSG) and less 
computational 
time 
Table 7. 


Literature summary of power conditioning in PV power generation. 


investigation involves the work till 2016. The above survey vigorously focussed 
on the analysis and requirements of an efficient PV power generation system. In 
order to yield the maximum energy from the PV panel which is highly nonlinear in 
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nature, in spite of its heavy capitation cost and low conversion efficiency a well- 
developed PV power generation system is essential. 

The boost converter which is intended between the PV panel and the load plays 
a major role in this process of extracting maximum energy. It should boost and 
regulate the voltage from the source, should have less switching losses semiconduc- 
tor switches, reduced ripples, should act as an isolation between source and the load 
and finally needs to be free from EMI. These forth-coming conclusions are derived 
from the above bunch of literatures that can be considered while designing a PV 
power generation system. 


e Si (crystalline) is considered as the best technology which produces nearly 
24.7% of efficiency when compared with rest of the evolved solar cell technol- 
ogies. Where, in this technology the electron-hole pair recombination seems to 
less than any other technology. 


e Though, soft switching technique fetched a considerable improved status 
in conversion efficiency there are few limitations, therefore care should be 
taken to reduce the size and cost of the system due to extra external circuit for 
control of switches. 


e Though, interleaving technique is best solution to get rid of input and output 
current ripple without the extra external circuit which is considered as the 
drawback of soft switching technique mentioned above, care should be taken 
for the source load mismatch. 


e In order to overcome the above-mentioned limitation of source load mismatch- 
ing in interleaved converters maximum power point tracking MPPT can be 
enabled to the system. With MPPT reduced ripples and harmonic distortions 
in the system high conversion efficiency is obtained. Care should be taken to 
choose the best MPPT technique. 


e There is high frequency ripple caused by the switching frequency and low 
frequency ripple caused by AC load. For the symmetrical 3-phase load, 
no low frequency ripple occurs. For single-phase or unbalanced 3-phase, 
there will be a double line frequency ripple. In this survey, high frequency 
ripple caused by switching frequency is concentrated. This is handled by 
the researchers in two ways either by tuning the switching style or by some 
unique ripple mitigation methods and fixed it with some additional inte- 
grated circuit. 


e Implementation of controller is also a best method to eliminate ripple. Not only 
for ripple reduction, the controller also helps to improve the systems overall 
performance by maintaining good dynamic response. 


e ILBC is better than the other special purpose converters like ZETA and SEPIC 
in terms of its system efficiency, compact size, ripple reduction, design 
simplicity. 


The survey mainly focused on the requirements to build an efficient PV system 
with extraction of maximum energy and better dynamic response. The general dis- 
cussion on objectives, techniques and advantages of each sections were also exhib- 
ited. The detail requirements of different constraint are presented in a systematic 
way in order to facilitate the reader to find a better solution for their application and 
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constraints. It is essential to enhance these methodologies used so far and provide 
sufficient guidelines for selection criteria to support the researchers of PV power 
generation to opt and design a best PV architecture. 
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the Grid via Three-Time-Scale 
Singular Perturbation Technique 
with Performance Analysis 


Youssef Mchaouar, Abdelmajid Abouloifa, 
Ibtissam Lachkar and Mohammed Fettach 


Abstract 


This chapter addresses the problem of controlling single-phase grid-connected 
photovoltaic system through a full bridge inverter with L-filter. The control objec- 
tives are threefold: (i) forcing the voltage in the output of photovoltaic panel to 
track a reference. This reference has been obtained from the maximum power point 
tracking strategy; (ii) guaranteeing a tight regulation of the DC-link voltage; and 
(iii) ensuring a satisfactory power factor correction (PFC) at the grid such as the 
currents injected must be sinusoidal with the same frequency and the same phase as 
the grid voltage. The considered control problem entails several difficulties includ- 
ing: (i) the high dimension and strong nonlinearity of the system; (ii) the changes in 
atmospheric conditions. The problem is dealt with by designing a synthesized 
nonlinear multi-loop controller using singular perturbation technique, in which a 
three-time-scale dynamics is artificially induced in the closed-loop system. A formal 
analysis based on the three-time-scale singular perturbation technique and the 
averaging theory is developed to proved that all control objectives are asymptoti- 
cally achieved up to small harmonic errors (ripples). The performance of the 
proposed approach and its strong robustness with respect to climate changes 
are evaluated based on the various simulations results carried out under 
Matlab/Simulink software. 


Keywords: single-phase grid-connected photovoltaic system, nonlinear control, 
three-time-scale singular perturbation technique, MPPT, power factor correction, 
averaging theory, stability analysis 


1. Introduction 
Due to dramatic increase in energy consumption and thrust to reduce carbon 
and greenhouse gas emissions from the traditional electric power generation sys- 


tems, photovoltaic (PV) power generators have gained a great popularity in recent 
years. Indeed, photovoltaic systems produce electric power without harming the 
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environment, transforming a free inexhaustible source of energy, solar radiation, 
into electricity. Furthermore, the major advantage of the photovoltaic systems is to 
meet the basic power requirement of non-electrified remote areas, where grid 
power has not yet reached. Also, there are other advantages such as the declining 
cost and prices of solar modules. On the other hand, the importance of PV systems 
in the solar industry makes these systems more efficient and reliable, especially for 
utility power in distributed generation (DG) at medium and low voltages power 
systems [1]. All these considerations assure a promising role for PV generation 
systems in the near future. 

On the other hand, many technical problems, such as untimely failures, could be 
found on electronic systems related in particular to the transfer and conversion of 
this energy to the network. Today, most conversion systems often suffer from low 
yields in real production sites. To meet the requirements of the new international 
standards on expected performance on associated conversion systems, it is impor- 
tant to make a research effort to solve the many control problems associated with 
the static power converter and bring this area to a degree of sufficient maturity to 
make them industrial products in their own right. One of the difficulties caused by 
the use of a photovoltaic conversion chain is focused on the problem of non-perfect 
control of the chain between the photovoltaic generator itself and the continuous or 
alternative type of load. 

The efficiency of a PV plant is affected mainly by three factors: the efficiency of 
the PV panel, the efficiency of the static power converter and its control, and the 
efficiency of the maximum power point tracking (MPPT) algorithm. 

PV grid-connected systems represent the most important field applications of 
solar energy [2-4]. In general, the power converter interface, from PV module (the 
DC source) to the load or to the grid, consists of two-stage converters: The first- 
stage DC/DC converter is usually used to boost the PV voltage and to implement the 
maximum power point technique. While the second stage is used to convert this 
power into high-quality AC voltage, with power factor correction (PFC) respecting 
to the power supply grid (i.e. sinusoidal and in phase with the AC supply voltage). 

Maximum power point tracking (MPPT) is required to match the PV array 
power to the environmental changes achieving to extract the maximum power 
output from a solar cell [5]. To this end, different MPPT techniques have been 
proposed such as incremental conductance [6], perturbation and observation 
(P&O) [7, 8], the hill-climbing, and some other special methods, such as neural 
networks, fuzzy logic technique [9]. Among all available techniques, a simple and 
effective MPPT of incremental conductance algorithm is applied to attain the max- 
imum power of PV array in different solar irradiance and temperature condition 
parameters. 

In order to provide a stable controller of DC/DC and DC/AC converters, many 
linear control methods have been proposed using many methods, such as a fuzzy- 
proportional integral controlled [10], a simple PR controller [11] where the perfor- 
mances have been illustrated by experimental result. However, in both proposed 
controllers, the problem of maximizing PV power transfer is not accounted for in 
the controller design. In contrast to linear control, nonlinear approaches can opti- 
mize the dynamic performance of system, such as sliding mode [12], fuzzy-sliding 
mode [13], feedback linearization [14], singular perturbation technique [15], and 
many others works [16-18]. In light of the previous descriptions, no theoretical 
analysis is made to formally prove that the closed-loop control performances are 
actually achieved. 

In this chapter, a multi-loop nonlinear controller is designed and developed via 
singular perturbation technique (Chapter 11 in Refs. [19, 20]), as was shown in 
Refs. [21, 22], where three-time-scale dynamics is artificially induced in the closed- 
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loop system. The control objectives are threefold: (i) achieving the MPPT for the PV 
array; (ii) ensuring a tight regulation of the DC-link voltage; and (iii) ensuring a 
grid connection with unity power factor (PF). These objectives must be met despite 
changes of the climatic variables (temperature and radiation). A theoretical stability 
analysis, for the closed-loop system, is provided using the three-time-scale singular 
perturbation technique [23, 24] and averaging technique (Chapter 10 in Refs. 
[19, 25]. The three-time-scale analysis allows to construct a suitable composite 
Lyapunov function candidate for the closed-loop photovoltaic system, and the 
stability properties of the resulting subsystems are analyzed providing mathemati- 
cal expressions for the upper bounds of the singularly perturbed parameters. 
Compared to previous works, the contribution of the new nonlinear controller 
enjoys several interesting features including the following: 


Several control objectives are simultaneously taken into account such as: 
MPPT, DC regulation, and PFC, whereas only some of these objectives have 
been tackled in previous works [10, 11]. 


A theoretical analysis will prove, using three-time-scale singular perturbation 
and averaging technique, that the desired multiple objectives are achieved. 
Such a formal analysis was missing in the previous works [12-15]. 


The nonlinearity of the controlled system was preserved in the controller 
design in order to keep all the properties of the studied system, whereas it is 
partly or totally ignored in previous controllers [14]. 


By including of three-time-scale dynamics in the full-order closed-loop system 
can ensure to achieve desired properties, such as robust zero steady-state error 
of the reference input realization, desired output performance specifications 
(overshoot, settling time), and insensitivity of the output transient behavior 
with respect to parameter variations and external disturbances. 


The content of this chapter is outlined as follows: in Section 2, the grid- 
connected PV system is described and modeled. Section 3 is devoted to the cascade 
nonlinear controller design and its performances are formally analyzed in Section 4. 
The global performance of the closed-loop photovoltaic system will be illustrated by 
numerical simulation using MATLAB/SIMULINK tool in section 5. A conclusion 
and a reference list end the chapter. 


2. System description 


This section describes the modeling of photovoltaic system connected to the 
grid. The power circuit topology used in the proposed single phase grid connected 
to the photovoltaic array is shown in Figure 1. It consists of the following compo- 
nents: (i) a photovoltaic array which consists of an arrangement of N,-series and 
N,,-parallel strings; (ii) an input capacitor C,, and a DC-DC boost converter used to 
increase the voltage level and achieve MPPT for photovoltaic array; (iii) a DC link 
capacitor Cy,; and (iv) a single-phase full-bridge inverter including four power 
semiconductors with L filter that is used to provide the energy to the grid and 
ensure power factor correction. 

Typical (Ip-Vp) characteristics of solar cells arranged in N,,-parallel and N,- 
series can be found in many places (see, e.g. [26]). The PV array module considered 
in this paper is of type KC200GT. In this chapter, a simple and effective MPPT of 
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Boost Converter Full-Bandge Inverter L-Fiber 


Figure 1. 
Single phase grid-connected PV system. 


incremental conductance algorithm is applied to attain the maximum power of PV 
array in different solar irradiance and temperature condition parameters. 

By analyzing the circuit and applying the well-known Kirchhoff laws, the system 
of Figure 1 can be described by the following set of differential equations: 


di; Vov Vde fi. 
dt Li G a L; L ( ) 
AVyy T z ; 
z = Ga (ipo a i) (2) 
di, Ve Te, | Vde 
ae iy, Eye T (3) 
Avac 1, li 
= | 4 
Ai G,, te (1 — u1) CG, (4) 


where vp, and ip, are, respectively, the photovoltaic generator voltage and cur- 
rent. vg, and iy, are, respectively, DC link voltage and current. i; designates the input 
current chopper, Cy, is DC link capacitor, r; and r, are, respectively, the equivalent 
series resistances (ESR) of input inductance L, and the filter inductance L,. v and i, 
are, respectively, the voltage and current of the grid. Here, the grid voltage is 
defined by v; = Eş sin (wt), where E, and œ; denote the constant amplitude and the 
constant angular frequency. The switching functions „4 and y, are defined by: 


1 if so is ON 
= 5 if so is OFF © 
1 if (s1,84) is ON and (s2,s3) is OFF 
Ma ={ if (ss4) is OFF and (s2,s3) is ON (6) 


The instantaneous model (1)—(4) cannot be used directly for the development of 
continuous control laws since it involves, as input variables, the binary signal u, and 
#5. To overcome this inconvenience, the average model is used [27]. Therefore, the 
state variables i;, vp», ig, and vg, are replaced by their average values x1, x2, x3, and 
x4 over a cutting period. The control inputs u1 and u2 denote the average values of 
41 and po, respectively. 


dx, 1 1 ri 1 (7) 
dt L; “ L; a Li ae Li a! 
dx ip 1 
= ees 8 
a Co Gy (8) 
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dx; Yg Vg 1 
dt Eep ie (9) 
dx4 1 1 1 
= 10 
di G3 - T + oF (10) 


3. Controller design 
3.1 Input inductor current regulation and PFC objectives 
3.1.1 Control law design 


The first control objective is to enforce the photovoltaic voltage x to track, as 
closely as possible, the optimal point Vy (called Regulator 2). However, it is well- 
known that the boost converter has a non-minimum phase feature. Such an issue is 
generally dealt with by resorting a cascaded loop design strategy that starts with the 
input current loop (Regulator 1), as it is shown in Figure 2. More specifically, the 
controller makes the input inductor current x; to track a reference signal x, the 
latter is determined from (Regulator 2). 

In parallel with the input current controller, the network current controller 
(Regulator 3) will be designed for power factor correction requirement that 
amounts to forcing the network current x3 to match the reference signal of the form 
x} = p sin (a,t). It means that the grid current x3 should be sinusoidal and in phase 
with the AC grid voltage v,, with £ is a signal. In fact, the latter is allowed to (and 
actually will) be time-varying but it must converge to a positive constant value. For 
both objectives (i.e. the input current x; regulation and PFC objectives), let us 
consider Eqs. (7), (9) in the following form: 


X13 = MıU + My (11) 
x4 ri X2 X4 
x4 Li ~ t Li Li 
where X13 = f »Mi=| x, pM = Tg v » and 
3 = = ——X3 — S: 
G S 


u=(%). 
u2 

Remark 1. Under condition that the DC-link voltage remains all the time 
positive, it can be showed that 0 < M1 min < M1 < Mimax < œ and |M2| < Momax < 00 


are satisfied. 
Let us introduce the following current tracking errors: 


e&eu=xř —x; (i=1,3) (12) 


Then, the reference model can be constructed in the following form 
= D;(xř,xi) (i=1,3) (13) 


where T1 and T; denote the time constants and they are selected based on the 
desired settling time, respectively, for the currents x; and x3. Based on (13), the 
realization errors of the desired behaviors of x; and x3, namely Q4 and Q3, are 
given by 
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Boost Single-phase 
Converter Inverter 


First Equation (11) 


al 


Figure 2. 
Schematic diagram of the proposed controller for single phase grid-connected PV system. 
Q; = D; (xe ,Xi) = xy (i = 1, 3) (14) 
Therefore, the control problem lim e;(t) = 0 (i = 1,3) corresponds to the 
— 00 
insensitivity condition defined by 
Q,;=0 (i=1,3) (15) 


Doing so, the behaviors of x; and x3 with prescribed dynamics of expression (13) 
will be fulfilled. Replacing in Eq. (14) x; and x3 by their expressions (7), (9), and 
since the requirement (15), one gets 


D'—-M,—M,U =0 (16) 


Hence, there is an isolated root called the inverse dynamic solutions for U given 


by 


L; ($ ( Tj i 1 1 
| id T PTR 
yid = & ee ae (17) 


Le e3 X4 T Ug + 19%X3 
X4 T3 i Le 


The control variable, namely U, has emerged in Eq. (16). At this point, an 
appropriate control law with the first derivative in feedback has to be found, so that 
the (e1, e3)-systems are made asymptotically stable. As the objective is to drive the 
error to zero, it is natural to choose the Lyapunov functions candidate 


V3\ _1/ (Qi())” 
& “2 ae Me 
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From expressions (14) and (16), it can be easily checked that the following time- 
derivatives 


0V3duy Q40Q, du, Qıx4 dur 
d s Ou dt = Ou, dt = L; dt (19) 
dt\ V4) | Vadum | | Q30Q3du. | Q3x4 duz 
ðu dt Ou, dt Lg dt 
are made negative-definite using the following control laws: 
p es 
u EVE 
a ') |k í ee 
u 
2 13 o 
E1E2 


At this point, kı and k3 are any design parameters, £1 and e, are small positive 
parameters. In view of M, > 0, it is easily seem from expressions (11) and (14) that 


kr 4 
d (V3 €, Li 
nA = — <0 Qi 0 21 
S k3 x495 ee) 1) 
€€, Ly 


for a positive values of k; > 0 (i = 1,3). Therefore, from expressions (14) and 
(20), the discussed nonlinear control laws are formulated as follows 


k ($ =) 
d [m Ti dt 
£1627 (") = : es dx; (22) 
3 T3 dt 


3.1.2 Singular perturbation system of the inner current loops 


Consider the closed-loop system of inner loops composed of the Eqs. (7)-(10) 
and the control laws (22), which can be rewritten in the following form 


dZ 


€1€9 — = A(X, Z,t 23 
aZ =h Z) 03) 
dX 
with Z = (z1 22)" = (uyu2)", X = (x1x2X3x4) , 
X4 "iy, , X2 X4 
GD be de. L 
e x4 ri X2 X4 xı ip 
k | | = us 
1G L; a Le L; r) Cp Cp 
h(X,Z,t) = , f (X,Z,t) = A 
e83 xX Vg g g 4 X4 
k H H =- y Haz 
1G u E is) ae a 
x41 xı X3 
Cac Cac Gr” 
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Now, we go to the fast time 7, = t/e,€2. Then for £1 = 0, the ultra-fast dynamic 
subsystem (UFDS) is given by 


a _ 4(X,Z,t) (25) 
dt 
dX 
oe 2 
a 0 (26) 


After the rapid decay of transients in expression (25), the steady state (more 
precisely, quasi-steady state) tends toward an equilibrium Z° = h(X). The manifold 
defined by Z’ is called the slow manifold, which is given by 


Li (3 x4 — X2 + a) 
zi T L; 
z= ( ) E a (27) 
25 L; e3 sA Vg + 1gX3 
X4 T3 L; 


Remark 2. 


i. During the fast transient in expression (25), the variables X are treated as the 
frozen parameters. 


ii. The equilibrium point given by expression (27) involves a division by the DC 
link voltage x4, from a practical point of view this division is not a problem 
because the DC link voltage remains all the time positive for the power 
converter to work correctly. 


By substituting of this equilibrium Z° into Eq. (24), the slow dynamic subsystem 
(SDS) of inner loops takes place on the slow manifold, according to the equation 


dx Cp Cp ( 2 8) 


Lox3 (2 Ve + =) Lixı (= $ rixı — 2) 
Cack 4 T3 Le CacX4 Tı Li 


Proposition 1. Consider the closed-loop system composed of Eqs. (23) and (24). 
For 1 — 0, the system takes the singular perturbation form where the UFDS is 
defined by equation (25), while the SDS of inner loop is defined by equation (28). 
Under the considerations given by Remark 2, one has the following properties 


i. If the design parameters k; (i = 1, 3) are positives, the UFDS (25) will be 
exponentially stable, and Z converge exponentially fast to Z°. 


ii. The behaviors of x; (i = 1,3) are prescribed by the stable reference equations 
of the form dx; /dt = (ae? — xi) /T;. Then, the requirements lim ei = 


lim (x;* — x;) = 0 are maintained. 
t—=œ 
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3.2 MPPT and DC bus voltage regulation objective 


The second step consists in completing the inner control loops by outer control 
loops for PV voltage (Regulator 2) and DC-link voltage (Regulator 4). The aim is 
now to enforce the photovoltaic voltage x and the DC-link voltage x4 to track, 
respectively, the optimal point x} = Vj, and a given reference voltage xj = Vo 
such that tuning laws for the ratio J and x¥ must be designed. According to the 
three-time-scale design methodology that is employed in this work, the general 
formulation of the three-time-scale singular perturbed systems requires the system 
to possess three different time scales. To this end, the voltage loops will be slow 
compared to the transients of the current loops. Therefore, the design parameters 
for voltage loops, in particular (e2, T2, and T4) must satisfy: 0 < e¢2 « £2 «1, and 
(Ti, T3) < (T2, T4). In addition, the steady states for the current x; and x3 yield, 
respectively, x; and xy. Therefore, the SDS of inner loop given by equation (28) 
will be reduced to 


ipv B xi 
dX, d (*) =: as (29) 
dt dt = oh N (gp \ Pee »\? 
X4 XX1 — VgX3 — fi (xi ) Tg (x3 ) 


Cacx 4 


Now, the first step is to establish the relation between the ratio 2 (which acts as 
the control input of the outer loop) and the DC-link voltage x4 (representing the 
output of the outer loop). 


3.2.1 Relation between (p) and (x4), and the control law 


The relation between the ratio # and the DC-link voltage x4 is the subject of the 
following proposition. 

Proposition 2. We consider the second equation of expression (29) and the power 
factor correction requirement defined by x} = f sin (a,t). 


i. The relation between x4 and £ is described in the following results 


2 2 292 
* * 
dx4 _ Xx] =i ar) vp rozp 
dt CacX4 Cacx4Eg Cacx aE 


(30) 


ii. Therefore, the squared-voltage x,,4 = Ke varies, in response to the tuning ratio 
B, according to the following first-order time-varying nonlinear equation: 


2 2 2 g2 
dXr,4 _ g2 =f; (wf) 2v,P 2rgveP 


31 
dt Cac CacEg Cack, 2 l } 


The second step is to establish control laws for the outer loops, in which 
Y = (y,y,)' = (xž P)” represent the new control inputs, while X, = (x„2%,4)" = 


T ; ; : f 
(x2x3) represent the new output variables. To this end, introduce the following 
tracking errors 


ej =x] —x,; (j=2,4) (32) 
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Then, let the desired behavior of X, be assigned by 


Bei OB (ge pen 
z =F Dhe) (j=2,4) (33) 


The error of the desired dynamic realization it follows 
e jp (F=24) (34) 


Then, the insensitivity condition is given by 
Q;=0 (35) 


Similar to the previous subsection and bearing in mind the fact that p and their 
first derivative must be available, and in order to meet the requirement (35), we 
should apply the control law given by the following structure 


a « 7 Gea — dx7,/dt) ) (36) 
at \y,)  \kalea/Ta — dx;,4/dt) 


3.2.2 Singular perturbation system of the outer voltage loops 


Combining expressions (31) and (36), and the first equation of expression (29), 
one obtains 


é2 Vy lov 
b( + ) 
Y T Cov 
met = 9(X,,Y,t) = on Te (37) 

t k e4 on 

a (= s) 
dX. ceed 
Sa S"& Vat) = | Co Co” (38) 

TI 


with f} =2 manny _ Bete yates The fast dynamic subsystem (FDS) is 
= Cai CucE? : y: y 


obtained by transforming the slow time-scale t to the fast time-scale t2 = t/e€2, then, 
by setting £2 = 0 


dY 


aes =g(X,„Y,t) (39) 
dX, 
ae = (40) 


Notice that expression (39) has an isolated equilibrium at Y° = (y$ y$), which 
will be determined (in the mean) in Appendix. As the FDS (39) is nonlinear, the 
stability properties of its equilibrium can be checked through the analysis of the 
Jacobian matrix of the linearized version defined as follows 


k2 
0 
Ar = Cp (41) 
kaAl, k4A" 
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Aren y5, 
Ca? 
we conclude that all eigenvalues of Ap satisfy Re (%4,2) < 0 for the negative values of 
k2 and k4. Therefore, Ar is Hurwitz matrix. By substituting of the equilibrium Y° 

into expression (38), the reduced SDS of outer loops takes place on the slow 
manifold, according to the equation 


dX, = €2/T2 
dt Ge a 


2. 
. 1 4 x aN W 
with Ap Gis Cè Ap = ‘aE + 


. Taking into account that y5 is positive, 


Proposition 3. Consider the system closed loop composed of expressions (37) and 
(38). For £2 — 0, this system takes the singular perturbation form, where the FDS is 
given by expression (39) and the reduced SDS of outer loop is given by expression 
(42). One has the following properties 


i. If the design parameters kz and k4 are negative, the FDS (39), will be 
exponentially stable and Y converge exponentially fast to Y°. 


ii. The behaviors of x,,; (j = 2, 4) are prescribed by the stable reference 
equations of the form dx,,;/dt = (x 7 —%r, j) /T;. Then, the requirements 


lim e; = lim (x3 — %,,) = 0 (j = 2,4) are maintained. 
t-00 t= oœ J É 


4. Control system analysis 


The objective of the global stability of closed-loop system can be analyzed in the 
following theorem. It is shown that the control objectives are achieved (in the 
mean) with an accuracy that depends on the network frequency @, and the small 
parameters ¢® (i= 1,2). 

Theorem. Consider the overall control system composed of the Pv panel, boost 
DC-DC converter and DC-AC inverter, described by the model (7)-(10), in closed 
loop with the multi-cascade multi-loop composed of: 


¢ The inner regulators (23), where (€122, k1, k3, T1, T3) are the design 
parameters; 


¢ The outer regulators (37), where (e2, k2, k4, T2, T4) are the design parameters. 
Then, one has the following property 


i. The augmented state vector Z = (u1 uz)" = (z1 22)", Y = (xf B)’ = (719) 


T T À ; 
and X’ = (x{x4x4x/,)" = (x1x2x3x4) undergoes the following state equations 


£1€.Z(t) = h(X', Y,Z,t) (43) 
e&f (t) =g(X',Y,Z,t) (44) 
X(t) = f(X',Y,Z,t) (45) 


with 
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L; Li Li 
u(i Zz1y/xh 7 rix] — xh + f/x) ome ipv 
Tı Li Li Cov Cov 
- VEA Ve + ToX! f= x4 Vg To a i 
k3 €3 45 £ E L @ TE S 
T3 L Ly o g g 


5X21 ri weya — 2WyVeV> 
Cac CuEg  CacEy 


g = 
ala wy £ rovg Y3 5X2) — "i 
CacEg CacE; Cac 
ii. Let the control design parameters be selected, such that the following 
inequalities hold T1 < T2, T3 < T4, kı > 0, k2 < 0, k3 > 0, k4 <0, and 
0 < ere: Ke: X1. 


Then, there exist oo constants p®, ¢® (n = 1,2), such that for and 
0<p<p®0<e,<e® «1 (n =1,2) the sees (43)-(45) has an asymptotically 
stable (21/a,) periodic solution X(t, €„, p), Y(t, €n, P), Z(t, €n» p) that continuously 
depends on ¢, (n = 1,2) and p = 1/a,. Furthermore, when (¢, = 0 and p = 0), one 
has 


lim ZE En Pp) = XG lim ot Enp)=Yo, lim PAG Ens P) = Zò, 


€n—0,p—0 En 0, p> 0 En—0, p> 


(46) 


T 
T k Toy 
} x* * * * * Y * # nd Z* * * with 
where A 9 (xi X70 %3,0 xo) > (i 02 A »4 0 (Zio zžo) it 


* yy! * _ yl Gs oes X20 X10 ToX3.0 > 
X20 == Vm: X40 E Vlo 210 1- Veia + Viio > Zp oc = Vio 3 Ji 0 = 1pv,0> 
2 g 3 2 
PENI (E412) 2r (223 gipo.0—"i(ipn,0) ) 
J2,0 = Wr ? 
See Appendix for the proof. 


5. Simulation and discussion of results 


The experimental setup is described by Figure 2 and the nonlinear controller, 
developed in Section 3, including the control laws (23) and (37), will now be 
evaluated by simulation in MATLAB/SIMULINK platform using the electrome- 
chanical characteristics of Table 1. 

The numerical values used for the design parameters are chosen as follows: 
e1 = 5.81 x 107°, e = 7.21 x 1077, T1 = 2.1 x 104, T, = 9.3 x 10%, 

T3 = 9.51 x 10°, T4 = 1.43 x 107, kı = 7.1 x 107, ky = —2.56 x 107%, 

k3 = 2.3 x 10°, K4 = 8.1 x 10°. These values have proved to be suitable based on 
several trials respecting the singular perturbation technique. In this simulation, we 
consider the KC200GT type of PV array module with N, = 54 and N, = 6. 

The performances of the proposed controller are illustrated by Figures 3-5. 
Figure 3 shows that the DC-link voltage x, is well regulated and quickly settles 
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Parameters Symbol Values 
Network E/f 220./2V /50Hz 
Boost Cpv 0.3mF 

Li 6mH 

ri 20mQ 
L-filter Lo 5mH 

Yo 50mQ 
PWM switching frequency fowm 14kHz 
DC capacitance Cy 7mF 

Table 1. 


PV system and single-phase grid characteristics. 
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Figure 3. 
DC-link voltage. 


400 


300- AN a rs, — W 


200 / \ a / \ 
| / / f 7 


do 
S 
T 
a 
N 
Wei 

wo 
S 
r 

1 


| 
0.34 0.35 0.36 0.37 0.38 0.39 0.4 
Time(s) 


Figure 4. 
Power factor correction checking. 


down after each change in the signal reference (stepping from 600V to 750V at 

t = 0.4s). The wave frame of the output current x3 is showed in Figure 4. The 
current is sinusoidal and in phase with the network voltage complying with the PFC 
requirement. This is further demonstrated by Figure 5, which shows that the ratio £ 
takes a constant value after transient periods following the changes in reference 
signals. 
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0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 


Figure 5. 
Tuning parameter p. 


5.1 Radiation variation effect 


Figure 6a-c illustrates the resulting closed-loop control performances in pres- 
ence of radiation changes. Specifically, the radiation takes a low, medium, and high 
value (equal to 800, 1000 and step to 600 W/m’ at times 0, 0.4, and 0.7 s, respec- 
tively), meanwhile the temperature is kept constant, equal to 298.17K (i.e. 25C’). 
Figure 6a shows that the captured PV voltage varies between 182, 212, and 138 V. 
These values correspond to the maximum points. Figure 6b shows that the DC-link 
voltage regulation is recovered after a short transient period following each change 
of the irradiation. Figure 6c shows that the current amplitude changes whenever 
the radiation varies. It is seen that the output current x3 and the grid voltage v, are 
actually sinusoidal and in phase. Hence, the converter connection to the supply 
network is done with a unitary power factor. 


5.2 Temperature variation effect 


The perfect MPPT is illustrated by Figure 7a. Here, the temperature steps from 
298.15 to 318.15 K, then to 308.15 K while the radiation 4 is kept constant equal to 
1000W /m?. Figure 7b shows that the DC-link voltage x4 is tightly regulated: it 
quickly settles down after each change in the temperature. Figure 7c illustrates the 
current amplitude changes whenever the temperature varies. The current remains 
(almost) sinusoidal and in phase with the network voltage complying with the PFC 
requirement. 


6. Conclusion 


In this chapter, an advanced controller is developed for PV grid-connected 
system. The latter is described by fourth-order nonlinear averaged model. The 
multi-loops nonlinear controller has been designed and developed using three-time 
singular perturbation technique and averaging theory. 

Using the theoretical analysis (via three-time-scale singular perturbation 
technique and averaging theory) and simulation, it is proved that the controller 
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Radiation variation effect: (a) Photovoltaic voltage. (b) DC-link voltage. (c) PFC checking. 


does meet the performances for which it was designed, namely: (i) Maximum 
power point tracking of PV array; (ii) tight regulation of the DC bus voltage; (iii) 


perfect power factor in the grid; and (iv) global asymptotic stability of the all 
system. 
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Figure 7. 
Temperature variation effect: (a) Photovoltaic voltage. (b) DC-link voltage. (c) PFC checking. 


Several simulation results have been made that illustrate the high performances 
of the proposed controller in ideal operating conditions (in the presence of meteo- 
rological constant) and its robustness against radiation and temperature change. 


Appendix (Proof) 


Part 1: Eqs. (43)-(45) are immediately obtained from expressions (37), (38), 
(23), and the first and the third equation of expression (24). 
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Part 2: The stability of the time-varying system (43)—(45) will now be performed 
in two steps using the averaging theory (e.g. Chapter 10 in Refs. [19, 25]) and the 
singular perturbation theory (e.g., Chapter 11 in Refs. [19, 23, 24]). The next step 
consists in using the averaging theory. To this end, let us introduce the time-scale 


change f = wet. Using this time, one gets e1¢2Z(£) = ph (X.Z ; i), 


asec Re epi as of L A 
&Y (t) = g(x, Y, i) , and X (t) = f(X, Y,Z, i) with p = 1/a,. In view of expres- 
sions (43)—(45), it is seem that the functions h te Z, i) g(x, Y, i) : and F(X’, Y,Z, i) 


as functions of t are periodic with period—2z, let us introduce the averaged functions: 


£1€2Z0(E) = plim z |F(X.2, f)at = pho(Xo, Zo) (47) 
ee en a pees aa 
ex¥o(f) = plim = g(x „Y, i) dé = py(Xo; Yo) (48) 
od! ee eee 
Xo(f) = plim = F(X: Y,Z, i) di = of (Xa Yo, Zo) (49) 


Since the systems here studied present equilibrium different from zero and 
in order to satisfy this requirement, a change of variables is introduced such 
that defines the new system in terms of its error dynamics. Therefore, the error 


dynamics are defined by introducing: Žo = (€1,0%,0)" = ((z1,0 =z o) (z2,0 -z% ae 


3 T T z 7 E - 7 
Yo= (Fro Foo) = (10 - Yio) (220 -= ¥o)) , and Xo = (X10 %2,0 %3,0 54,0) = 
((x1,0 — xž0) (%2,0 — x30) (%3,0 — X31) (%4,0 — X4.9)), where the constant xž o, X30» 
X39>X4,0> Jio Zo ZLo and 279 represent the desired average values of the state 
variables. Then, since expressions (43)—(45) and according to expressions (47)- 
(49), the system can be rewritten into its error-dynamics formulation thus defining 
the closed-loop error dynamics as: 


e162Zo = pho (Xo, Zo) =p in Doane (50) 
k (-2-Fo(to.20)) 
2 m ka — 22 —fro(%o.¥o) ) 
£2Y0 = pfo (Xo, Yo) =p 5 a (51) 
k4 (- z -Fro (Žo Žo) ) 
F 1,0 (X0,Zo) 
> AA nia fao (Šo Yo) 
Xo = pfo (Xo, Y0; Zo0) =p L> & (52) 
: faol 0» Zo) 
faol o Yo) 
where 
a os, (%2,0 + x39) — 4/ (%4,0 +7, as i (Zio +210) 4/ (%4,0 +7, 
fio(Xo Zo) 3 = ( i) ri(žno + xio) | = ( io) 


Li Li Li 
(53) 
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Tn > S j v y > +y 
fig (Xo, Yo) = T 50 ( 1 7 1 a (54) 
pv ‘po 


(#40 +x a) 


L 


Deea. ee 
F36(%020) = ii x30) 
g 


(%,0 +230) (55) 


Fo(Xo. Yo) = eta +y) 


(X2,0 + x3) (x0 + is) z ri( (o + sal) 
Cac 


+2 (56) 


In order to get stability results regarding the system of interest (43)-(45), it is 
sufficient (thanks to averaging theory) to analyze the averaged system (50)—(52). 
Now, the asymptotic stability of the resulting three-time-scale photovoltaic single 
phase grid system (50)—(52) is discussed, which is based on the sequential (double) 
time-scale analysis similar to the one presented in Refs. [23, 24], it is an extension of 
the two-time-scale analysis presented [19]. 

The use of theory of the three-time-scale singular perturbations for the stability 
analysis is based on the idea that, for 0 < £1£2 « €2 <1, the trajectories in Xo Yo 
and Žo of the system (50)-(52) can be approximated by three models: the slow 
dynamic subsystem (SDS) of full system, the fast dynamic subsystem (FDS), and 
the ultra-fast dynamic subsystem (UFDS). We can thus find Lyapunov functions 
for each one of the singularly perturbed subsystems. 

For the UFDS, it is necessary to ensure that the dynamic of expression (50) does 


not to shift from the quasi-steady-state equilibrium Žo = hy (X o). Then, the associ- 
ated Lyapunov function candidate is obtained by introducing a change of variables 


Zo =Zo— ho (Xo), so that its equilibrium is centered at zero. By letting t1 = t/e1€2, 
the UFSD in function of Zp is defined as follows 


£ 21,0 
dZ Apes! oA L; ? 
da, 7 o (Ro 2o) =p (57) 
k3 (#40 +xio) 
22,0 
Lg 


in which Xo is treated as a fixed parameter. Thus, the associated Lyapunov 
function can be defined by 


1 7 x 
Vy = 2p (Puii o + Py28}0) (58) 
In view of k; > 0 (i = 1,3) and (%4,0 + Xa a) > 0, it is clear that the solutions Py; 
and Py? can be chosen as 
L; 
Pu = qui (59) 
ey (#40 + Ae) 
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L 
Pi = 2 Gu2 (60) 


2k3 (#40 +230) 


with qy and qy; are positive. Similar to the ultra-fast subsystem, it is necessary 
to ensure that the dynamic of expression (51) does not to shift from the equilibrium 

= ¢ (Xo). Then, by introducing a change of variables Yo = Yo — g/ (Xo), so that 
its equilibrium is centered at zero and by letting t2 = t/e2, the ultra-fast subsystem 
can be rewritten in the following form 


- ko . 
dY fact anA , 
Te, ~ P(o Yo) =p Cp? (61) 
w" (Xo; Yo) 


; Nn rifo reo w Îao 2riCpv (X2,0 , tpv,0 (otio) 
with y” = 2k4 k +36, t G + (2 T, + GC, Ir0|- 


P 
Since expression (61), it is seem that the fast subsystem is aan According 
to the proposition 4 (part i), it is shown that this subsystem can be made asymp- 
totically stable by letting k; > 0 (for j = 2, 4). Moreover, we can find a Lyapunov 
function, which takes the following quadratic form 


1 n ^ A a 
Vp = 2p (Privo +P a Biao) (62) 


with the solutions to the associated Lyapunov given as 


2 
x ipv, X20 +23 
Cou =- 2) Z a) 1r2 
Pn =-— 
2k Roy! [Y q Crk 
2y" m t kC, 


Cac 
PR = ~ Akay? 1r (64) 


Cac (rcn (A z+ e) egsa) 


kay’ (2 F cate) 


Pr = 1ra (65) 


in which Xo is treated as a fixed parameter. Finally, the new slow dynamic 
subsystem is obtained by substituting of the ultra-fast subsystem equilibrium Žo = 


ho (X o) and the fast subsystem equilibrium Žo =f) (X 0) into expression (52) 


3 a S & T Pi oe, x n 
Xo = (X10 %20 %30 %40) =—p(%10/T1 %2,0/T2 %3,0/T3 %4,0/Ts)" 
(66) 


It is easy to define the associated Lyapunov function for SDS as follows 
Vs(%o) =—-XeP.Xo =— (Pax? 2 67 
s(Xo) = pe a (Psixi, o + PsoX2,0 + Ps3*3,0 + PsaX4,o) (67) 
where psy = T1451/2s Ps2 = T 2459/25 Ps3 = T34s3/2 Ps4 = T 4454/2 and 451 50> 


4539 4ş4 are the positive constants. 
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Based on these Lyapunov function candidates, the double application of the 
standard two-time-scale stability analysis is divided in two stages: in the first stage, 
the stability analysis focusses on proving the stability properties of the degenerated 
Š sp-subsystem (slow-fast subsystem). The results obtained will be used in the 
second stage in order to prove the stability properties for the full }°,,-system 
(slow-fast-ultra-fast subsystem). 

X SF Stability analysis 

In the first stage, the standard method for two-time-scale systems is applied in 
which the previously derived Lyapunov functions for the slow and fast subsystems, 
that is, Vs and Vp, respectively, must satisfy certain inequalities. 


* Isolated equilibrium of the origin for the X` o,-subsystem 

The origin (Xo = 0, Yo = 0) is an isolated equilibrium of the }>,,-subsystem 
0 = g,(0, 0) 

0 = f, (0, 0, ho (Xo) ) 


(68) 
(69) 


moreover, Yo = (Xo) is the root of 0 = gy (Xo; Yo) which vanishes at Xo = 0, 
and ||¢, (Xo)| < 91(||Xoll) where 94() is a x function. 


* Reduced system condition for the X` çp-subsystem 


Using the Lyapunov function (67) and substituting pf , (x o Yo, ho(X 0)) yields 


gE 
OVs(Xo Oe Se eee ee - 
(=g) Pfo (Xo, Yo, ho (Xo)) < — am} (Xo) (70) 
0 
where a; < 1 and ¢,(Xo0) = y3 Kio HRE o + 8X3) + Xo. 
* Boundary-layer system condition for the X` .,-subsystem 
Since expression (62)-(65), it is seem that 
Wio N 
uF |Xo, Zo SIG oz ‘ 7 
(megea) pH(Xo, Žo) < — QrI30 cs QrIio (71) 
0 
where Qp and Qp are defined as follows 
4m (iy +) E) 
Qa = rı _ "ilr ( poe ee i EAE (72) 
a= a DERA Y1,0| T |20 
(2v + fe) 
aC (ay +98) — Sep 
_ 42 _ "sfr ( le : f z (73) 
Qr = L dy PERA Y1] + |¥2,0 
(2v Eee) 
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For physical point of view and domain of working principle, it is supposed that 
all physical state variables are bounded in domain of interest, where 
Xo min < Xo < Xomaxs Fomin < Yo < Yomaxs and Zo min < Zo < ZoMax: Therefore, Qr 
and Qp can take positives minimum possible values Q p1, min > 0 and Q Fz, min > O. 
This can be done by ensuring the appropriate selection of 4;,, 7, k2, and k4 as 
follows 


4lr:C x tpv,0 (%2,0Max + X34) 

k2 Cp ipo 2, 0Max Cp 2 : a 

kag Cac Di V1,0Max — 2” min 
y, 22,0Max 

z 
x ipv,0 (%2,0max + x39) 
4 Gis ps 2 
. | p (am) 7 
qrı > Vj p 5 E 
lpr Vinin T Cacko 1,0Max * Y2,0Max 
min k4Cpv 


(74) 


Therefore, the boundary-layer system condition for the ` ṣp-subsystem is 
defined as follows 


2 


OVEN coe g E) -2 er 
Wo P£o (Xo, Yo) <- / Qromin T + Qrimin 1,0 So arp (Yo) (75) 


where as <1 and ,( Yo) = y Qramin J3, 0 + Qrimin It, (on 


Now, we consider the composite Lyapunov function candidate of the $- sp - 
subsystem given as follow 


Vs = (1— 11) Vs(Xo) + mV (Yo) (76) 


with 0 <7, <1. The derivative of Vsp presents new terms which represents the 
effect of the interconnection between the slow and fast dynamics. These intercon- 
nections are assumed to satisfy the following conditions 


> (fo (Xo, Yo, ho(Xo)) - fo (Xo. (Xo); ho(%o)) ) < %1d2(¥o) bs (Xo) 


aX 0 
(77) 
(=e) pfo (Zogo (Xo), ho (Xo)) < hh (Yo)¢, (Xo) (78) 
(=) lofo (žo, Fo, ho (Xo)) = rf (Zogo (Xo), ho(%o) | < 113 (Yo) 


(79) 


where the constants 01, 02, and 7, are non-negative. 
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Therefore, from the singular perturbation theory (e.g. Theorem 11.3 in 
Ref. [19]), it follows that the derivative of Vsp is negative-definite for 


62 < e8 for (m <n?) (80) 
with e2 = — “2 — andn? = t= (81) 
2 ay, + 0102 M1 0, + O2 i 


X SFU Stability Analysis: 

The stability of full system is analyzed now. In this step, the results obtained in 
above section will be used in order to prove the asymptotic stability properties of 
the full $ sry-system, which, for convenience, is first rewritten as 


Žo = Fo (Zo Žo) (82) 
eeo = ho (žo Žo) (83) 


i 
where fo = (210 X20 3,0 X40 Jio Foo) , similarly to the $` ṣp general 
asymptotic stability analysis presented in the above section, the X ` şpy-system is 
treated like a two-time-scale singularly perturbed system, where the X sp- 
subsystem is treated as the new slow augmented reduced order. The stability of the 
full system implicates that the previously derived Lyapunov functions for the new 


slow and new fast subsystems, that is, Vst (o) and Vy (Zo), respectively, must 
satisfy certain inequalities. 


* Isolated equilibrium of the origin for the X` ṣpy-subsystem 


The origin (fọ = 0, Žo = 0) is an isolated equilibrium of the 7 .,,j-subsystem 


0 = F,(0, 0) (84) 


0 =ho(0, 0) (85) 


Moreover, Žo = ho (Zo) is the unique root of 0 = ho (Zo Žo)» which vanishes at 
ho(žo)| < 92(||%o||) where 92() is a x function. 


Žo = 0, and | 
* Reduced system condition for the X` spy-subsystem 


Using Lyapunov function candidate Vgp(7) given by expression (76), it is easily 
shown that 


T 

OVsr%o)\' _ (IVs Zo) OVsr(%o) Vso) ƏVse(Zo) Vso) Vsr(Zo) 

OX 0X 1,0 0X 2,0 0X3,0 0X 4,0 dJi o dJ 0 
(86) 


Recalling that Yo = Yo — £ (Xo), and since (61)—(65), one gets 
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~ 7 A g ! À 
1 2C, rC, x i 2r;i 
Cn (2G (FE +) = (20 -+40)) ( (1 i r)a Cpl 2,0, pw ilpv 
OP ry _ 2 2 


‘pv 
ox 2,0 


kay’? (w T (Cack2/2k4Cpv) ie 


z ipv m A 2r;Cpv 
2Cp | 2riCpv T; ! C, (%2,0 + xžo) Tə 
A V 
x (2y + (Cack2/2k4Cpv) ) how" (v aL (Cacka /2k4Cpv)) qr = A192 


f 2 
; x tpv o4 * 1 
app, "CCm (rcp G z) T x0) (y + (Cacks/2ksC)) 
0% 4,0 


Ir = Arq 

QkyT ay? (w’ + (Cack2/2kaCp))? ue FA 
2 riCm \ ~ xy toy Wiipy 

Pm _ CeCmrsl&(1-TB)0+ GE FE) _ 

i — — 3 ’ 
da0 Akay" 42 qr 
2 

OP pa 1eCi 


2 7 Š ipv x I ý 1 
pn Gri (2rcp E | e) (ia sio) ) Qv + (Cackn/2esCyn)) 


0X 4,0 


> 


=A ; 
AkaT ay"? (w’ + (Cack2/2k4Cp))? e A 


g i C, nC, ga: i 2r;i 
C WiC x mw \ fe ¥ 2 (4 ivpv \ ~ LC 2,0 'pv itpv 
ann | elite) C) (2 C E t ea n 
0X2,0 Akaw'? (w’ + (Cack2/2k4Cpv) ) 


2r; Cp i\c 
Ta = ‘dc 


ł =A 
Akay! (y! + (Cack2/2k4Cpv)) 4r2 34ra 


x (2y’ + (Cuck2/2k4Cpv)) 


: 2 
$ 1 
3 (2u 3+ z) x +x; ) 
oP o = Ge P 4 T2 Cp ( 2 žo) 


kT 2’ (w! + (Cack2/2k4Cpy) ) “a 


C, m x. i 2r;i 
A C ipv fi 2,0 pv itpv 
Arq _ Cae (28 (1 T2 ) 20 Cw (E Ton T )) A 
Prz 7 Vr. = 4189 
Fi Fao Akay"? F2 F2 
x ipv 7 . 
5 Cova | 2Cr | a Bo —(%2,0+%3) T 
Di = ‘pv = J1,0 
F3 do aT akay ("+ (Cacka/2aCp)) lr = A94 o> Pry dao 0 
ee One i 
1 Ga(2i +2) - earo) ekl- Eao cie —*8)) 
92,0 T2 Cw po is 
= o 
PP O59 Akay (y + (Cack2/2k4Cpo)) Irn = A109 py 
2 x ipv pa 
0 . Ci 2riCpv T, + C — (204x30) 
p 231.0 = 0 p 220 = 2 pe dm = A 24 
F1 Oa, > E F3 dX4,0 AT akaw? (w'+ (Cacka/2k4Cp)) F2 129 F2? 
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IIo = Ci, =A 
Pio — 4kaTyy? F2 T ANI p 


Substituting the expressions of pFo (to ho o) and by using absolute value 
version of Young’ inequality |ab| < | (a? + b°) 


oV? 
eai pFo(Žosho(Zo)) < — Mix; o — Mox5 9 — M35 o 


-M453 o = M510 z M63 o 


(87) 


2 
where My = (1 - m), M: = (1- m) 2-5" Gao) 4 ( $) , 
2 
M3 = —(1—m) 8, My = (1 ny) = su (Bq we 4 ( +) i 


2 
a ; 1 =n) = 1 2y’ T: 
M; = tp. — M; -6 1 and Me = 1n M; -7 ( Ca Ca oP) 0) 


Max 


A139 + ih, + AsJo9 + (A7 + Aro) )4r2 


| 
ry 
oF 19F2 
Max 


(Fino + Aoy20 +42) (ate R = i) 
(Fino + Aefa o + An) (- a 5 dina) 
iota (Gs Ee) 


x i v0 = 
2 ipo T, A ai y i 
( r Cp ( 7 + =) (2,0 +o) )r 


52,0 + Ag+ Ag 
2Cpv 


M; =m 


Max 


2 
A134 F1) Max + Mar 


Max 


+ ip 


+19 2 
Max 


2 


A2 z 
T NE + Abr +An) ir) 


Max Max 


Vi. 
+19 F2 x I2, F 


2y' (w’ a (Cack2/2k4Cpy )) 71,0 
Max 
49,,+As+ Ao n A3 2 
M. — n, 122 22,0 om > Ae +A 
6 nı 2 oe igs 2 2 Joo + 3+ 9 2 TE 
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Now, it needs to ensure that —MiX7 o — M283 0 — M353 0 — Max o — Ms Ion 
Me Po <0. This can be done by ensuring the appropriate selection of qs1; 452» 453» 
4542 pp and qp» such that qs; and q,, have sufficient large positive values, q,, and 
qp have sufficient small positive values, and q 5 aa iss are limited as follows 


i Aee 14, Aeri 1 fl HVEN 
y T <q < w c ie wa) < as, < E (for m < 1). There- 
F 


T: 


ier 


fore, the reduced yan condition for the X` sry ge is defined as follows 


(ales)! pFo(ZosholZo)) < ahr) (és) 


where a3 < 1 and 


p(o) = Mikio + M253 o + M353 o + M434 o + M51 o it M6950: 


e Boundary-layer system condition for the }>.,,,-system 


Using the Lyapunov function given by expression (58)—(60), and substituting 
pho (Xo, Zo) of expression (57), one obtains 


Vu y Z f a 5 
(EE) whine) = Mg R s ahs) E 
0 
where a4 < 1 and 4(Zo) =\/™@ Zio + tig 0: 


Now, we consider the Piette Lyapunov function candidate of the full system 
given as follow 


Vsru(%o. 20) = (1 — m)Vsr (Žo) + n2Vu (Zo) (90) 


with 0 < m <1. The derivative of Vsru(79,%0) presents new terms, which rep- 
resents the effect of the interconnection between the slow and fast dynamics. These 
interconnections are assumed to satisfy the following conditions: 


2 50 E 

(“se p(Fo(zo»Zo) — Fo (Zoho (o))) < Osba(Zo)s(Zo) OD 
A EE E A ere 

apy. © plo (Zor holo) ) < 0494 (Zo) bso) (92) 


z T 
ƏV q (Žo ma he a Doa 2 
(= p [Fo (ZoŽo0) - Žo (Zoho (£o) )| < 19% (Zo) (93) 
0 
the constants 63, 64, and y, are non-negative. Therefore, from the singular 
perturbation theory (e.g. Theorem 11.3 in Refs. [19, 23, 24]), it follows that the 
derivative of Vsru is negative-definite when 


ey<e® for (n, <n?) (94) 


* @ _ _ aay @__ 
with ey = tam and ny = ay: 
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f : ae T 
Therefore, it can be inferred that the equilibrium Xj = (x¥oxž0xž0xž0) > 


T 
Y> = ( Viodz, o) „and Zò = (zřozž Di of full system, is asymptotically stable for all 


E&i < EP (i = 1,2). Therefore, Part 2 immediately follows from the averaging theory 
(e.g. Theorem 10.4 in Ref. [19]). Theorem is established. 
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Chapter 10 


Battery Internal Fault Monitoring 
Based on Anomaly Detection 
Algorithm 


Nassim Sabri, Abdelhalim Tlemcani and Aissa Chouder 


Abstract 


Battery internal faults are one of the major factors causing safety concern, 
performance degradation, and cost increases. To extend the lifetime of the battery 
and bring more security in the system, internal fault detection of solar battery is 
proposed in this paper using an unsupervised machine learning algorithm based on 
anomaly detection method. The advantages of adopting such a method consist of 
using unlabeled data that meet the battery case in the difficulty of obtaining the 
fault data. In contrast, healthy data can easily be obtained from the battery and 
therefore allows building the anomaly detection algorithm. The effectiveness of the 
proposed method is validated using a simulation platform of a stand-alone photo- 
voltaic system developed in Matlab/Simulink that takes as system input a real 
profile of irradiance and temperature captured from the Centre de Development 
des Energies renewable (CDER), Algeria. The test results in real-time data show the 
ability of the proposed approach to detect the fault occurrence in the battery. 


Keywords: anomaly detection, internal faults, battery, simulation, data 


1. Introduction 


The future development of stand-alone photovoltaic (SAPV) systems will know 
a progression in all countries, especially in remote areas and islands [1]. The battery 
is the main component in SAPV system, since it represents 43% of lifecycle costs for 
this system [2], where this lifecycle could significantly reduce in badly work condi- 
tions. To check continuously the status of battery and aging, monitoring of internal 
resistance is proposed in [3, 4] proposed an estimation of solar irradiation to mon- 
itor the overcharge and internal resistance in Battery. Other studies [5-8] focusing 
to assess the performance of the whole SAPV system using energy parameters. This 
latter provides an analysis of SAPV system, where the battery takes a part in this 
system. Fault diagnosis in SAPV system including the specific fault battery external 
short-circuit is proposed in [9, 10] using Artificial Neural Network, the study is 
experimentally validated an additionally integrated with the system using Matlab 
Graphical User Interface (GUI) [11]. Battery performance evaluation in PV-diesel 
hybrid system is proposed in [12] where the analysis requires many days and 
provides inaccurate results concerning the battery state. 

The failure developed in battery, in particular, the internal faults could bring a 
difficulty to eliminate them by electrical protection elements such as fuse or 
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circuit-breaker by reason of absence an appropriate selection of fuse for this circuit 
[13]. This will result in undetected faults that conduct the battery to work in low 
performance and decrease its lifecycle, furthermore, it increases the risk of fire 
hazard and even an explosion, without account the contamination of equipment’s 
beside [13]. To solve this problem and overcome the limitation of existed battery 
fault analysis, this paper proposes internal faults detection of solar battery using 
anomaly detection technique. This method is commonly utilized for fault detection 
[14] specifically in the photovoltaic system [15, 16], where the proposed anomaly 
detection prediction model has the ability of development with only normal data 
and using the readily available battery voltage and current without requiring extra 
sensor circuit. Furthermore, it could recognize any non-specific failure in battery 
which makes them a promising choice for practical application. 

The training and test for validation have been performed using data from the 
simulation platform of SAPV system, the battery is adequately represented in 
Matlab\SimPowerSystem to allows the creation of fault such as short-circuit and 
ground-fault. In addition, to make the simulation varies along with irradiance level 
and temperature, a real climatic measurements profile captured from CDER is used 
as input to estimate the output. 

The test results of fault detection using anomaly detection method in real-time 
data show the capability of anomaly detection to recognize all the faults in the 
battery. 


2. Modeling of SAPV system 

A typical SAPV system as shown in Figure 1 consists mainly of PV panels, 
Battery, load, and charge regulator [17]. The PV panels are the DC power source, 
and the battery is the storage element. 


2.1 Output PV panel modeling 


The output current and voltage of PV panel is obtained using the common 
one-diode model [18], in which the relationship I-V is given as follow: 


V +R V+R,I 
t ) 1|- = (1) 


I=I -I ! 
etg 


where I,; is the light generated current, Io is the reverse saturation current of 
diode, n is the diode ideality factor, R, and R,, are the series and shunt resistance of 
panel respectively, and V+ is the thermal voltage. 


— Charge regulator | 
+ 


+ + 
3 Battery iis Load 
PV — 
panels æ 


Figure 1. 
Typical SAPV system. 
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2.2 Battery modeling 


The charge and discharge models of lead-acid battery implemented in 
SimPowerSystems simulation environment [19] are summarized in the following: 


2.2.1 Model for charge 
Viatt =Eo —R-i ea koe it + exp (t) (2) 
2.2.2 Model for discharge 
Voatt = Eo —R- i-ke (it +i*) + exp (t) (3) 


where Vpatt is the voltage of battery (V), Eo is the constant voltage of battery 
(V), k is the polarization resistance (ohm), Q is the capacity of battery (Ah), it is the 
present battery charge (Ah), 7 is the current of battery (A), i is the filtered (A) 
current and exp. (t) is the exponential zone voltage (V). 


2.3 Load and charge regulator representation 


The load is represented by a resistance connected in parallel with the PV panels 
and battery. The charge regulator consists of a simple switch on/off placed in one 
side between PV panels and battery in our case. 


3. Proposed fault detection method 


Anomaly detection or known as novelty detection or outlier detection [20] is one 
the most machine learning technique used in fault detection [21], which aim to 
detect abnormalities or unusual operation that can come up, it makes the assump- 
tion that the data are distributed according to Gaussian (or Normal) distribution 
and this latter can be modeled based on two parameters: the mean p and the 
variance 0°. 

Three phases are required to build an anomaly detection model, the first is 
training phase, where the Gaussian distribution is estimated by the parameters u 
and o°, the second is validation phase, in this step, some threshold e is selected as 
a limit of being an anomaly (outlier) compared to the Gaussian probability 
function, the third is testing phase, in which a test is performed to check the 
performance of the model. In the following, Mathematical equations behind 
this approach are given. 


3.1 Gaussian distribution 


Fit a model of Gaussian distribution from data relies on the assumption that 
huge number of data used is normal data. For each feature xj (Gal, non) an 
estimation of Gaussian distribution parameters p; and oj. The Gaussian probability 
density is defined as: 


Pla) =- exp (= =m?) (4) 


189 


Advanced Statistical Modeling, Forecasting, and Fault Detection in Renewable Energy Systems 


3.2 Estimation parameters for Gaussian distribution 


The parameters of Gaussian distribution y; and 0; for the j-th feature are esti- 
mated respectively as follow: 


Hj = S (5) 
oe = 59 (x? — mi (6) 


3.3 Choosing the value of threshold ¢ 


To decide if new examples are anomalies or not, the procedure consists to 
compare the value of Gaussian probability density with a threshold or limit, if the 
probability is lower than a certain threshold then these examples are anomalies. 
However, the selection of the threshold e could be done by the trial-error method 
using the F1 score (Eq. (7)) metric as criteria. As shown in Figure 2, the algorithm 
proceeds to try several values of e, where the chosen value of e corresponds to the 
maximum F1 score [14] defined below: 


_ 2-prec- vec (7) 
prec + rec 


where prec is the precision and rec is the recall, they can be obtained by: 


tp 
rec = (8) 
d tp + fp 
— P 
rec = Pri (9) 


Initialization 
(£, F1_Max) 


F1=F1_Max 


li 
4 


Figure 2. 
Flowchart of the procedure to select the threshold. 
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Where tp is the true positive that means that the algorithm successfully classified 
as positive or anomaly, fp is the false positive that means that the algorithm incor- 
rectly classify it as positive or anomaly, when it is not an anomaly, finally fn is the 
false negative that means that the algorithm wrongly classify it as not anomaly, 
while this sample is an anomaly. 


4, Case study 


A small SAPV system composed of 2 PV panels (212 Wp) connected with a 
battery of 1200 Ah and a load of 50 W is considered in this work. The PV panel used 
is Isofoton 106 W-12 V, where the parameters for this panel are obtained from [22]. 
Industrial lead-acid battery of 12 V contains six cells connected in series, in order to 
investigate internal faults in the battery, each cell is represented by a single battery 
of 2 V connected in series to model a real battery of 12 V, we point out that the cells 
are assumed has identical electrical characteristics. 

The object of this paper is detecting internal faults that occur in the battery, 
where two faults are considered: ground fault and short-circuit (Figure 3), in which 
ground fault situate at three locations: upper, middle and lower, and short-circuit 
between cells consist of 1, 2, 3, 4 or 5 cells circuited. These faults are used to test and 
evaluate the detection approach. 


Figure 3. 
SimPowerSystem model of SAPV system with battery. 


5. Results and discussion 
A real environmental measurement (irradiance and temperature) illustrated in 


Figure 4, taken from CDER, Bouzareah, Algeria, In which nine clear days are used 
as input to simulate the SAPV system under Matlab/Simulink. 
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In the following, some scenarios of faults are performed in order to test the 
capability of anomaly detection to recognize these faults using the last 4 days of data 
constituted of battery voltage and current, The Gaussian probability density is 
calculated and plotted with the threshold e before and after occurring the fault. 

In Figure 5 a test for the Normal operation is realized, it can be seen that the 
probability does not drop below the threshold e apart three false alarms noticed, 
which means there is no fault detected by the detector system. 


0, 
z 
T boi 
= 2” : 
5 ž 
LE ishi 
E £ 
E E 
zor: 
8 a Od dA 72 AA 10 ddd IR 19 15 g au Az 72 oa im i 16 «1m 
Time ni Tine th] 
(a) b 
Figure 4. 


Captured irradiance and temperature. 
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Figure 5. 
Simulation test result: normal operation. 
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Figure 6. 
Simulation test result: ground fault upper cell. 
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As illustrated in Figure 6, when a ground fault in the upper cell occurs, the 
probability decrease under the threshold e after a certain time, and Figure 7 shows a 
ground fault in the middle cell, which takes longer time than the previous fault 
before to be detected. 

In Figure 8 a short-circuit between 2 cells is created and as illustrated this fault 
take much time to descend under the threshold g, while in Figure 9 a similar fault is 


dij ıı 1 1 1 
Calculated prob ability 
the thresheld epsilon 
OE 
0.06 Gee Fault securs 
ati 
0.04 
, Fault detected 
0.02 _ ett? 


g 


$ a4 45 72 56 4128 4144 165 4142 
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Figure 7. 
Simulation test result: ground fault middle cell. 
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Figure 8. 
Simulation test result: short-circuit of 2 cells. 
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Figure 9. 
Simulation test result: short-circuit of 2 cells (Rf x 0). 
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Figure 10. 
Simulation test result: short-circuit of 5 cells. 


created with Rf ~ 0, then the probability drop below the threshold instantly and 
detected. 

Figure 10 indicates the detection of short-circuit of 5 cells after a specific time of 
occurring this fault. From all these figures, it shows that all the faults are detected, 
whereas some faults take much time to detected than other, and this by the fact that 
the faulted voltage and current of Battery is reduced gradually and not immediately, 
and this depend on the external circuitry at moment of fault, in our case this is due 
to the value of Rf. 


6. Conclusions 


In this paper, a detection of battery cell interconnection fault is proposed using 
anomaly detection algorithm. The method has the benefits of using only the steady- 
state of battery and uses the easily available battery voltage and current to predict 
the internal fault in battery. Based on simulation data, the anomaly detection is 
developed and tested for validation, where in the simulation environment, the 
battery unit is viewed as series-connected battery cells. In this way, several typical 
faults such as internal short and ground fault are carried out. The proposed method 
is capable to effectively predict the battery internal faults, where the analysis 
finding reveal that as more battery cell are involved at faults or the fault circuit has 
negligible resistance, the fault detection becoming much faster to indicate the fault 
occurrence. The future work consists to implement this method on real battery as 
well as extend the application of battery fault detection to include battery EV and 
other appliances. Furthermore, another method based on a statistical approach to 
select the threshold e in a better way would be proposed in the future, these 
techniques will perfectly manage to find the outlier of battery operation. 
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Edited by Fouzi Harrou and Ying Sun 


Fault detection, control, and forecasting have a vital role in renewable energy 
systems (Photovoltaics (PV) and wind turbines (WTs)) to improve their productivity, 
efficiency, and safety, and to avoid expensive maintenance. For instance, the main 
crucial and challenging issue in solar and wind energy production is the volatility of 
intermittent power generation due mainly to weather conditions. This fact usually 
limits the integration of PV systems and WTs into the power grid. Hence, accurately 
forecasting power generation in PV and WTs is of great importance for daily/hourly 
efficient management of power grid production, delivery, and storage, as well as for 
decision-making on the energy market. Also, accurate and prompt fault detection and 
diagnosis strategies are required to improve efficiencies of renewable energy systems, 
avoid the high cost of maintenance, and reduce risks of fire hazards, which could 
affect both personnel and installed equipment. This book intends to provide the reader 
with advanced statistical modeling, forecasting, and fault detection techniques in 
renewable energy systems. 
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